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INTRODUCTION 


Research  Initiation  Program  -  1987 

AFOSR  has  provided  funding  for  follow-on  research  efforts  for  the 
participants  in  the  Summer  Faculty  Research  Program.  Initially  this 
program  was  conducted  by  AFOSR  and  popularly  known  as  the  Mini-Grant 
Program.  Since  1983  the  program  has  been  conducted  by  the  Summer  Faculty 
Research  Program  (SFRP)  contractor  and  is  now  called  the  Research 
Initiation  Program  (RIP).  Funding  is  provided  to  establish  RIP  awards  to 
about  half  the  number  of  participants  in  the  SFRP. 

Participants  in  the  1987  SFRP  competed  for  funding  under  the  1987 
RIP.  Participants  submitted  cost  and  technical  proposals  to  the 
contractor  by  1  November  1987,  following  their  participation  in  the  1987 
SF<JP. 

Evaluation  of  these  proposals  was  made  by  the  contractor. 
Evaluation  criteria  consisted  of: 

1.  Technical  Excellence  of  the  proposal 

2.  Continuation  of  the  SFRP  effort 

3.  Cost  sharing  by  the  University 

The  list  of  proposals  selected  for  award  was  forwarded  to  AFOSR  for 
approval  of  funding.  Those  approved  by  AFOSR  were  funded  for  research 

efforts  to  be  completed  by  31  December  1988. 

The  following  summarizes  the  events  for  the  evaluation  of  proposals 
and  award  of  funding  under  the  RIP. 

A.  Rip  proposals  were  submitted  to  the  contractor  by 

1  November  1987.  The  proposals  were  limited  to  $20,000  plus 
cost  sharing  by  the  universities.  The  universities  were 
encouraged  to  cost  share  since  this  is  an  effort  to 

establish  a  long  term  effort  between  the  Air  Force  and  the 
university. 

B.  Proposals  were  evaluated  on  the  criteria  listed  above  and 

the  final  award  approval  was  given  by  AFOSR  after 

consultation  with  the  Air  Force  Laboratories. 

C.  Subcontracts  were  negotiated  with  the  universities.  The 

period  of  performance  of  the  subcontract  was  between 

October  1987  and  December  1988. 

Copies  of  the  Final  Reports  are  presented  in  Volumes  I  through  III 

of  the  1987  Research  Initiation  Program  Report.  There  were  a  total  of  83 

RIP  awards  made  under  the  1987  program. 
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STATISTICS 
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Total  SFRP  Participants  159 
Total  RIP  Proposals  submitted  by  SFRP  117 

Total  RIP  Proposals  submitted  by  GSRP  7 

Total  RIP  Proposals  submitted  124 

Total  RIP'S  funded  to  SFRP  81 

Total  RIP'S  funded  to  GSRP  2 

Total  RIP'S  funded  83 

Total  RIP'S  Proposals  submitted  by  HBCU's  11 

Total  RIP'S  Proposals  funded  to  HBCU's  7 
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LIST  OF  UNIVERSITY  THAT  PARTICIPATED 


Adelphi  University  -  1 
Alabama,  University  of  -  1 
Alaska-Fairbanks,  Univ.  of  -  1 
Alfred  University  -  1 
Arizona  State  University  -  1 
Arkansas  State  University  -  1 
Arkansas,  University  of  -  1 
Auburn  University  -  1 
Bishop  College  -  1 
Capital  University  -  1 
Catholic  Univ.  of  America  -  1 
Cedarville  College  -  1 
Central  State  University  -  1 
Cincinnati,  University  of  -  5 
Colorado,  University  of  -  2 
Dayton,  University  of  -  7 
Dillard  University  -  1 
Drury  College  -  1 
Eastern  Illinois  University  -  1 
Eastern  Kentucky  University  -  1 
Eastern  New  Mexico  University  -  2 
Fairfield  University  -  1 
Florida  A&M  University  -  1 
Florida,  University  of  -  2 
Fort  Lewis  College  -  1 
Gonzaga  University  -  1 
Grambling  State  University  -  1 
Hampton  University  -  1 
Houston,  University  of  -  2 
Howard  University  -  1 
Idaho,  University  of  -  1 
Illinois-Chicago,  Univ.  of  -  2 
Indiana  University  -  1 
Indiana  Univ.  of  Pennsylvania  -  1 
Iowa,  University  of  -  1 
Jackson  State  University  -  1 
Jarvis  Christian  College  -  1 
Jesm  Baromedical  Res.  Inst.  -  1 
John  Hopkins  Evening  College  -  1 
Kansas  State  University  -  1 
Kansas,  University  of  -  1 
Kentucky,  University  of  -  1 


Meharry  Medical  College  -  1 
Memphis  State  University  -  1 
Metropolitan  State  College  -  1 
Michigan  State  University  -  1 
Mississippi  State  University  -  4 
Mississippi,  University  of  -  1 
Missouri-Kansas  City,  Univ.  -  1 
Missouri-Rolla,  Univ.  of  -  3 
Montana,  University  of  -  1 
Montclair  State  College  -  1 
Morehouse  College  -  1 
Nazareth  College  -  1 
Nebraska-Lincoln,  Univ.  of  -  2 
New  Mexico  State  University  -  1 
New  York  State,  Univ.  of  -  3 
N.  Carolina  A&T  State  Univ.  -  1 
N.  Carol ina-Greensboro,  Univ  -  1 
Northwestern  University  -  1 
Ohio  State  University  -  5 
Ohio  University  -  2 
Oklahoma  State  University  -  1 
Oregon  Institute  of  Tech.  -  1 
Oregon  State  University  -  1 
Ouachita  Baptist  University  -  1 
Pace  University  -  1 
Pennsylvania  State  Univ.  -  1 
Point  Loma  College  -  1 
Puerto  Rico-Mayaguez,  Univ.  -  1 
Purdue  University  -  1 
Rochester  Inst,  of  Tech.  -  1 
Rose-Hulman  Inst,  of  Tech.  -  2 
Saint  Paul 's  College  -  1 
San  Francisco  State  Univ.  -  1 
South  Dakota  State  Univ.  -  1 
South  Florida,  University  of  -  2 
Southeastern  Mass.  Univ.  -  2 
Southern  Illinois  University  -  2 
Southern  Mississippi,  Univ.  -  1 
Southern  University  -  2 
St.  Louis  University  -  1 
St.  Mary's  University  -  1 
Talladega  College  -  1 
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Lock  Haven  Univ.  of  Pennsyl-  -  1 
Long  Island  University  -  1 
Louisiana  State  University  -  1 
Louisiana  Tech.  University  -  1 
Lowell,  University  of  -  4 
Texas  Southern  University  -  3 
Texas  Technical  University  -  2 
Texas-Austin,  University  of  -  1 
Tuskegee  University  -  1 
Utah  State  University  ■■  1 
Walla  Walla  College  -  1 
Washington  State  University  -  1 
West  Florida,  University  of  -  1 
Western  Michigan  University  -  3 


Taylor  University  -  1 
Temple  University  -  1 
Tennessee  Technical  Univ.  -  1 
Tennessee,  University  of  -  1 
Texas  A&M  University  -  2 
Wichita  State  University  -  2 
Wilberforce  University  -  1 
Wisconsin-Eau  Claire  Univ.  -  2 
Wisconsin-Madison,  Univ.  of  -  i 
Wisconsin-Whitewater,  Univ.  -  1 
Wittenberg  University  -  1 
Worchester  Polytech.  Inst.  -  2 
Wright  State  University  -  3 
Xavier  University  -  1 
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University  of  Florida 
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University  of  Dayton 
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HARRY  G.  ARMSTRONG  AEROSPACE  MEDICAL  RESEARCH  LABORATORY  / 
(Wright-Patterson  Air  Force  Base) 

(continued) 

Or.  Noel  Nussbaum  Dr.  John  Westerkamp 

Wright  State  University  University  of  Dayton 
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Duke  University 

Specialty:  Biomechanical  Engineering 
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(Arnold  Air  Force  Systems) 
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Eastern  Illinois  University 

Specialty:  Aerospace  Engineering 
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Tennessee  Technical  University 
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Or.  Narayan  C.  Haider 
University  of  South  Florida 
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AVIONICS  LABORATORY 
(Wright-Patterson  Air  Force  Base) 
(continued) 

Prof.  William  K.  Curry 
Rose-Hulman  Inst,  of  Technology 
Specialty;  Computer  Science 
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Wright  State  University 
Specialty;  Computer  Science 

Dr.  George  W.  Zobrist  (1986) 
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Specialty;  Electrical  Engineering 
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(Tyndall  Air  Force  Base) 

Dr.  William  W.  Bannister 
University  of  Lowell 
Specialty;  Organic  Chemistry 

Dr.  William  M.  Bass 

The  University  of  Tennessee 

Specialty:  Physical  Anthropology 

Dr.  ^  ter  Jeffers 
S.U.N.Y. 

Specialty:  Chemistry 

Dr.  William  T.  Cooper  (1986) 
Florida  State  University 
Specialty:  Chemistry 


FLIGHT  DYNAMICS  LABORATORY 
(Wright-Patterson  Air  Force  Base) 

Mr.  Thomas  Enneking  (GSRP) 
University  of  Notre  Dame 
Specialty:  Civil  Engineering 


Dr.  Alastair  McAulay 
Wright  State  University 
Specialty;  Electrical  Engineering 

Dr.  John  Y.  Cheung  (1986) 
University  of  Oklahoma 
Specialty:  Electrical  Engineering 


Dr.  William  Schulz 
Eastern  Kentucky  University 
Specialty:  Chemistry 

Dr.  Joseph  Tedesco 
Auburn  University 
Specialty:  Civil  Engineering 

Dr.  Dennis  Truax 
Mississippi  State  Univeristy 
Specialty:  Civil  Engineering 

Dr.  Yong  S.  Kim  (1986) 

The  Catholic  Univ.  of  America 
Specialty:  Civil  Engineering 


Dr.  Gary  Slater 
University  of  Cincinnati 
Specialty:  Aerospace  Engineering 
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FLIGHT  DYNAHICS  LABORATORY 
(Wright-Patterson  Air  Force  Base) 
(continued) 

Dr.  Oliver  McGee 
Ohio  State  University 
Specialty;  Engineering  Mechanics 

Or.  Shiva  Singh 
Univ.  of  Kentucky 
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University  of  Dayton 
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Dr.  Forrest  Thomas 
University  of  Montana 
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Or.  William  Wolfe 
Ohio  State  University 
Specialty:  Engineering 
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University  of  Portland 
Specialty;  Engineering  Mechanics 
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(continued) 
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Iowa  State  University 
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II.  CORRELATIONS  OF  SPONTANEOUS  IGNITION  TEMPERATURES 
WITH  MOLECULAR  STRUCTURES  OF  FLAMMABLE  COMPOUNDS 
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FINAL  REPORT  TO  UNIVERSAL  ENERGY  SYSTEMS,  INC.,  DAYTON,  OHIO 

William  W.  Bannister,*  Stuart  B.  Clough  and  Sukant  Tripathy 
Department  of  Chemistry,  University  of  Lowell,  Lowell,  MA  01854 

ABSTRACT 

As  noted  in  the  previous  paper  of  this  series,  a  reversed 
relationship  exists  between  molecular  weight  and  spontaneous 
ignition  temperature  (IT) ,  occasioned  by  contact  with  hot 
surfaces;  up  to  a  point,  less  volatile  higher  molecular  weight 
fuels  have  lower  IT's  and  are  more  easily  ignited.  For  higher 
members  of  the  alkane  family  this  trend  may  reverse,  resulting 
in  minimum  IT's  for  the  C5  -  Cg  alkanes.  Branched  chain 
alkanes,  arenas  and  olefins  also  have  higher  IT's  than  analogous 
straight  chain  alkanes.  Important  parameters  in  governing  IT's 
of  fuel  components  include:  (1)  the  effect  of  molecular  weight 
on  the  velocity  of  the  fuel  molecule  near  the  hot  surface; 

(2)  the  effect  of  rigidity  of  the  fuel  molecule  on  its  ability 
to  recoil  from  the  hot  surface;  and  (3)  the  effect  of  specific 
heat  of  the  molecule  with  regard  to  ability  for  dissipation  of 
heat  energy  within  the  molecule.  Molecular  modelling 
considerations  will  be  discussed  to  explain  anomalous  trends. 

As  can  be  seen  in  Figures  1  and  2,  there  is  a  decrease  in  IT  values 
or  alkanes  with  increasing  molecular  weight,  for  the  range  of  alkanes 
rom  methane  (€2^)  through  octane  (C0)  . 


Figure  1.  Spontaneous  Ignition 
Temperatures  for 
Typical  Hydrocarbons 

(References  1  and  2) 
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Precise  data  is  not  available  for  higher  alkanes.  As  shown  in  Fig.  2, 
in  the  lower  projected  track  of  IT  values,  one  possibility  lies  in  an 
asymptotic  sweep  along  a  200*^  isotherm  for  alkanes  beyond  Cg. 

(Thermal  cracking  may  occur  at  this  temperature,  if  the  hydrocarbons 
are  permitted  to  linger  at  the  heated  surface  in  accordance  with 
current  ASTM  procedures. 7- 8  Thus,  larger  alkanes  could  conceivably 
have  higher  ignition  temperatures;  but  by  undergoing  slow  pyrolysis  at 
200°,  some  alkanes  in  the  minimum  IT  region  characteristic  of  Cg  -  Cg 
would  form,  with  these  igniting  at  this  spuriously  low  temperature. 

Figure  2.  Correlations  of  Ignition  Temperatures  for  n-Alkanes 
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Frank  observed  increased  IT  values  for  C^^  through  C20  alkanes, 
attributing  this  to  decreasing  volatilities;’ and  several  other 
anomalously  high  IT  values  have  been  observed  for  decane  and  higher 
alkanes.  Thus,  a  minimum  zone  of  IT's  may  exist  for  Cg  -  Cg  alkanes. 
This  is  shown  in  Figure  2  as  the  upper  projected  track  for  the  alkanes 
beyond  Cg .  Obviously,  either  track  could  actually  pertain;  currently 
available  data  is  insufficient  and  gooc  values  are  probably  poorly 
available  with  existing  equipment  and  technique. 


Fuel  component  characteristics  which  have  been  identified  as  having 
possible  impact  on  the  ignition  temperatures  are  discussed  on  page  3 
and  illustrated  on  page  4  of  this  paper. 
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Molecular  weight  effects  on  speed  and  inertia  for  heavier 
molecules .  These  effects  increase  residence  time  in  the  vicinity 
of  a  hot  surface,  with  increased  time  for  energy  transfer. 

2.  Molecular  rigidity.  Rigid  molecules  have  reduced  rotational 
degrees  of  freedom  within  the  molecular  structures.  This  would 
significantly  decrease  "floppiness"  of  the  molecule  in  its  impact 
with  a  hot  surface,  allowing  the  more  rigid  structure  to  rebound 
readily  and  rapidly  with  correspondingly  reduced  residence  times 
in  the  vicinity  of  the  hot  surface,  and  considerably  reduced 
energy  transfer  from  the  surface  to  the  fuel  molecule. 

3.  Specific  molecular  heat.  For  smaller  molecules  relative  speed  and 
relative  rigidity  are  important.  For  the  lower  alkanes,  each 
addition  of  a  methylene  unit  imparts  a  significant  increase  in 
molecular  weight  and  j.n  internal  degrees  of  rotational  freedom. 

In  going  from  ethane  to  propane,  for  example,  there  is  a  17% 
decrease  in  speed  at  500®C,  whereas  decane  is  only  5%  slower  than 
nonane  at  this  temperature.  Internal  rotational  degrees  of 
freedom  are  even  more  considerably  enhanced  for  propane,  which  is 
very  flexible  compared  with  the  very  rigid  ethane;  but  decane  and 
nonane  have  very  little  difference  in  "mushiness"  due  to  this 
effect.  Thus,  decreased  speed  and  rigidity  arising  from  increased 
molecular  weight  are  much  less  important  for  the  higher  alkanes. 

Moreover,  higher  alkanes  begin  to  enjoy  the  benefit  of  increased 
"molecular  specific  heat".  When  lower  alkanes  such  as  propane  are 
energized  by  radiational  heat,  all  atoms  of  the  small  molecule  are 
fairly  equally  irradiated.  For  larger  molecules,  some  atoms  will 
be  in  the  shadow  of  others.  The  hotter  surface  atoms  can  then 
transfer  some  of  their  increased  energies  to  the  cooler  internal 
atoms.  Thus,  a  higher  ignition  temperature  will  be  required  to 
attain  decomposition  energies  for  a  molecule  large  enough  to 
provide  shadowing  by  some  of  its  atoms  to  other  neighboring  atoms. 
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Figure  3.  Molecular  Characteristics  Affecting  Ignition  Temperatures. 

Typical  JP-8  component  (A) :  Typical  JP-4  component  (B) : 

dodecane  heptane 

CK2-CH2-CH2-CH2-CH2-CH2-CK2-CH2-CH2-CH2-CH2-CK3  CH3-CK2-CK^-CH^-CHo-CH -CH3 


[(1),  (2),  (3)  ...  (8),  (9)  above  indicate  positions  of  molecules  A 
and  B  at  same  time  intervals,  relative  to  hot  surface  area.] 

Molecular  velocity  at  temperature  (T,  °K)  is  inversely  proportional 
to  the  molecular  weight  (M) ;  smaller  molecules  have  greater  speed. 

Average  velocity  =  14,600(T/M)^  cm/sec 

The  slow  moving  heavy  molecule  A  (above)  is  therefore  near  the  hoc 
surface  longer  than  the  lighter  and  faster  B,  and  A  has  more  time 
available  to  absorb  energy  from  the  hot  radiating  surface. 

The  larger  molecule  A  has  more  mobile  C-C  bonds  and  therefore  more 
rotational  degrees  of  freedom  than  the  smaller  B.  Thus,  A  will 
tend  to  be  "mushier"  with  less  elastic  recoil  on  impact  with  the  hot 
surface  than  the  smaller,  more  rigid  B  which  can  bounce  away  more 
readily.  .Again,  molecule  A  will  have  more  time  to  energite  than  B. 

To  offset  this,  and  although  the  larger  molecule  has  greater  surface 
area  enabling  it  to  absorb  more  heat  than  the  smaller  molecule,  the 
large  molecule's  volume  is  increased  even  more  than  is  the  surface 
area  (by  comparison  with  the  smaller  molecule) .  Absorbed  heat  is 
then  dissipated  more  through  the  large  .molecule,  which  accordi.ngly 
is  cooler  than  the  small  molecule.  This  effect  is  more  pronounced 
with  increased  molecular  weights  (and  volumes) .  In  this  case 
A  (dodeca.ne)  needs  more  heat  (with  a  higher  ignition  temperature) 
to  energize  it  to  its  decomposition  point  than  does  B  (heptane)  . 

For  smaller  m.olecules  (e.a.,  comparing  hexane  [Cg]  and  propane 
[Cg]),  this  offsetting  factor  of  increased  molecular  specific  heats 
IS  not  as  i.mportant  as  the  molecular  speed  and  rigiditv  effects:  in 
sucn  cases  tne  heavier  .molecule  has  the  lower  ignition' temcerature'. 
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INTRODUCTION 


The  estimation  of  stature  from  various  skeletal  elements  has 
been  an  area  of  critical  interest  to  physical  anthropologists  for 
many  years.  Past  studies  concerning  stature  estimation  have 
included  those  by  Manouvrier  (1892),  Stevenson  (1929),  Allbrook 
(1961),  Genoves  (1967),  Lundy  (1983)  and  Lundy  and  Feldesman 
(1987).  These  papers  presented  data  and  statistical  formulae  for 
the  estimation  of  stature  for  geographic  and  racial  groups  as 
diverse  as  the  Chinese  and  the  South  African  Black.  The  formulae 
most  often  employed  in  the  United  States  are  those  provided  by 
Trotter  and  Gleser  (1951a;  1951b;  1952;  1958)  and  Trotter  (1970). 
All  of  these  studies  have  shown  a  high  correlation  between  the 
length  of  any  whole,  long  limb  bone  and  stature,  with  the  highest 
single  correlation  being  the  femur. 

Forensic  anthropologists  are  often  confronted  with 
fragmentary  remains  while  the  estimation  of  stature  by 
conventional  formulae  is  dependent  upon  whole  limb  bones.  The 
estimation  of  living  stature  from  long  bones  is  based  upon  the 
principle  that  the  various  long  limb  bones  correlate  positively 
with  stature.  Since  this  is  true,  the  parts  of  each  individual 
long  bone  should  also  be  related  to  stature  even  tnough  they  may 
not  correlate  as  highly  as  the  length  of  the  whole  bone.  Stature 
estimation  from  fragmentary  long  bones  has  been  attempted  by 
Muller  (1935),  Steele  and  McKern  (1969),  and  Steele  (1970).  In 
practice,  these  methods  are  difficult  to  apply  because  of 
difficulty  in  identifing  the  exact  location  of  necessary 
anatomical  landmarks.  Steele  (1970:87)  himself  recognized  this 
problem  when  he  stated: 
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Another  landmark  not  used  in  the  present  study,  but  used 
previously  by  the  author  and  Dr.  McKern,  involves  the 
humerus  alone.  The  landmark  in  question  is  the  point  of 
greatest  narrowing  at  the  middle  of  the  diaphysis.  Although 
this  narrowing  can  be  related  to  a  point,  posterolaterally 
in  the  radial  groove  where  the  lateral  supracondylar  ridge 
becomes  indistinct  (on  a  level  with  the  distal  termination 
of  the  surface  for  m.  deltoideus) ,  in  practice  the  point 
proved  impossible  to  locate  with  sufficient  constancy  to 
give  useful  results. 

Most  practicing  physical  anthropologists  have  also  experienced 
difficulties  in  locating  many  of  the  other  anatomical  landmarks 
used  to  predict  stature  from  fragmentary  long  bones.  Indeed  many 
do  not  use  either  the  Steele  and  McKern  (1969)  or  the  Steele 
(1970)  formulae,  and  hence  do  not  attempt  to  estimate  stature 
when  no  whole  long  limb  bones  are  present  because  of  this  very 
difficulty.  However,  since  skeletal  remains  are  often 
fragmentary,  it  would  be  useful  to  devise  a  more  reliable 
technique  for  stature  estimation  from  smaller  segments  of  long 
bones . 


OBJECTIVE 

The  objective  of  this  study  is  to  assess  the  feasibility  of 
stature  estimation  using  small  fragments  of  the  femur.  As 
described  above,  the  Steele  and  McKern  (1969)  and  Steele  (1970) 
methods  are  plagued  by  the  difficulty  of  identifying  the 
anatomical  landmarks  which  define  the  various  segments  of  the 
limb  bones  from  which  estimations  are  calculated.  The  present 
approach  seeks  to  avoid  this  pitfall  by  using  standard,  clearly 
defined  measurements  taken  on  the  distal,  proximal,  and  midshaft 
regions  of  the  bone.  Many  of  these  measurements  are  already 
familiar  to  physical  anthropologists,  who  regularly  use  them  for 
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comparative  population  studies  or  in  some  cases,  the  estimation 
of  gender.  It  is  also  hoped  that,  by  presenting  data  and 
providing  equations  for  the  estimation  of  stature  from  small 
segments  of  all  areas  of  the  bone,  any  recovered  section  might  be 
used  to  estimate  stature.  It  will  normally  be  the  case  that 
smaller  segments  of  bone  will  exhibit  weaker  relationships  to 
stature  than  larger  segments.  Because  the  emphasis  of  this  study 
is  on  smaller  segments,  the  equations  generated  are  therefore 
expected  to  offer  less  accuracy  for  both  maximum  femur  length 
and  stature.  However,  when  small  segments  alone  are  available, 
this  approach  will  offer  the  best  opportunity  for  estimation  of 
stature.  The  relationship  of  these  dimensions  to  maximum  femur 
length  and  stature  has  never  been  systematically  examined. 

The  femur  was  selected  as  the  Initial  skeletal  element  to  be 
analysed  because  in  addition  to  its  high  correlation  with 
stature,  it  is  (along  with  the  feet)  one  of  the  bones  most 
frequently  recovered  from  Air  Force  crash  sites.  It  is  a  large, 
durable  bone  protected  by  both  large  amounts  of  soft  tissue  and 
the  seat  and  harness  mechanisms  of  the  aircraft. 

METHODOLOGY 

The  sample  measured  was  obtained  from  the  Terry  anatomical 
collection  housed  at  the  Smithsonian  Institution's  National 
Museum  of  Natural  History  in  Washington,  D.  C. .  The  skeletons  in 
the  Terry  collection  are  from  cadavers  dissected  at  the 
Washington  University  School  of  Medicine,  St.  Louis,  Mo. 
Accurate  data  for  age,  sex,  race,  and  cadaver  stature  are 
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available.  A  seunple  of  approximately  200  males  and  200  females 
from  Blacks  and  incites,  yielding  a  total  sample  of  800 
individuals  was  obtained.  To  be  included  in  the  sample  an 
individual  must  have  race  and  sex  recorded  on  the  morgue  record. 
Only  femora  exhibiting  no  gross  pathological  changes  of  the 
distal,  proximal,  or  midshaft  regions  were  included  in  the  study. 
Left  femora  were  measured  whenever  possible,  though  right  femora 
were  substituted  if  breakage  or  pathological  conditions  precluded 
the  use  of  the  left.  Preference  for  inclusion  in  the  sample  was 
given  to  individuals  for  whom  age  and  cadaver  stature  were  also 
recorded.  With  few  exceptions,  photographs  of  the  cadavers  were 
examined  in  order  to  be  certain  that  the  stature  was  recorded 
when  the  body  was  placed  vertically  with  the  feet  in  natural 
plantar  position. 

The  measurements  employed  are  defined  in  Table  1  and 
illustrated  in  Figure  1.  All  seven  measurements  are  standard  and 
defined  in  Martin  (1957) .  The  definition  of  number  7  was 
modified  from  a  midshaft  diameter  to  a  minimum  transverse  shaft 
diameter  in  order  to  avoid  the  necessity  of  locating  a  midpoint 
on  a  fragment.  All  the  measurements  can  be  recorded  with  ease. 
Maximum  femur  length  was  measured  on  an  osteometric  board?  all 
other  measurements  were  obtained  with  sliding  calipers.  All 
measurements  were  taken  by  Tal  Simmons,  a  Ph.D.  student  in 
physical  anthropology  at  the  University  of  Tennessee. 


(BCB) 


FIG.l 


STANDARD  MEASUREMENTS  OF  THE  FEMUR  USED  IN  THE 
REGRESSION  EQUATIONS  AGAINST  MAXIMUM  FEMUR 
LENGTH  (FML)., 


Table  1  :  Measurements  (abbreviations  in  parentheses)  used  in 
estimating  stature  from  femoral  fragments. 


Martin  No. 
1 

18 

15 

13 

7 

21 

25 

26 


Measurement 

Maximvim  Femoral  Length  (FML) 

Vertical  diameter  of  the  femoral  head  (VHD) 
Vertical  diameter  of  the  femoral  neck  (VND) 
Proximal  Breadth  of  the  Femur  (VHA) 

Trcmsverse  dicuneter  of  midshaft  (minimum)  (WSD) 
Bicondylar  Breadth  (BOB) 

Epicondylar  Breadth  (FDL) 

Lateral  Condyle  Height  (LCH) 

Medial  Condyle  Height  (MCH) 


Data  were  entered  on  an  IBM  compatible  PC  using  dBase 
Ill-plus.  Analysis  of  the  data  was  carried  out  using  the  SAS 
package  on  the  University  of  Tennessee's  Vax  cluster.  Equations 
for  the  prediction  of  stature  were  generated  in  two  ways.  First, 
stature  estimation  was  performed  by  the  usual  method  of 
regressing  stature  onto  bone  measurements.  Second,  because  the 
Terry  collection  records  do  not  always  list  the  cadaver  stature 
of  the  individual,  maximum  femur  length  was  also  regressed  onto 
bone  measurements.  By  determining  the  maximum  femur  length  of  an 
individual  from  a  fragment,  the  stature  may  then  be  estimated 
also,  though  the  standard  error  of  the  estimate  will  be 
increased. 


RESULTS 

Age,  Stature,  and  Maximum  Femoral  Lencrth 

Stature,  femur  length  and  age  of  the  present  sample  is 
compared  to  other  samples  which  have  been  used  in  stature 
estimation  in  Table  2.  The  means  and  standard  deviations  for 
age,  stature,  and  maximum  femur  length  are  quite  close  to  those 
presented  by  Trotter  and  Gleser  (1952)  for  the  Terry  Collection. 
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They  differ,  however,  from  those  presented  by  Steele  (1970) 
presumably  because  of  that  author's  small  sample  size  and  his 
elimination  of  older  individuals  (those  over  age  70)  from  his 
study.  Jantz  and  Moore-Jansen  (1988)  have  shown  that  the  Terry 
Collection  differs  from  the  contemporary  White  and  Black 
population,  present  day  populations  exhibiting  significantly 
longer  femora.  Table  2  also  shows  that  Terry  Collection  heights 
are  considerably  below  Trotter  and  Gleser's  (1958)  Korean  war 
dead  and  Modem  Forensic  cases.  Consequently,  regression  formulae 
developed  from  the  Terry  collection  and  applied  to  modern  people 
may  underestimate  stature. 

Sunmarv  of  Fragment  Measurements 

Table  3  presents  the  means  and  standard  deviations  for  each 
of  the  segments  measured.  Correlations  for  these  measurements 
with  femur  length  and  cadaver  stature  are  given  in  Table  4.  The 
correlations  show  the  strength  of  association  between  the  segment 
measurements,  stature  and  femur  length.  At  the  outset  we  had 
expected  that  VHA  (proximal  femur  breadth)  would  be  the  best 
predictor  and  that  is  borne  out  by  the  correlations  in  males,  but 
in  females  several  other  measurements  are  more  highly  correlated. 
In  general,  correlations  rarely  exceed  0.65. 

Comparison  with  Steele's  Results 

Our  goal  at  the  outset  was  to  attempt  stature  estimation 
from  fragements  which  would  equal  or  exceed  Steele's  in  accuracy 
but  would  use  standard  measurements.  We  can  examine  this  question 
by  comparing  our  segment  correlation  with  those  from  Steele  that 
are  more  or  less  comparable.  Steele's  segment  1  is  measured  from 
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Table  2.  Means  and  standard  deviations  of  ageCyrs) , 
stature* (cm.)  and  femur  length (cm.)  of  various  samples 
compared  to  present  sample. 


White  Males 


Terry (1) 

200 

58.73 

13.88 

167.75 

7.59 

45.64 

2.68 

Terry (2) 

255 

61.66 

12.55 

167.89 

7.34 

45.66 

2.45 

Terry (3) 

61 

52.97 

4.98 

168.44 

4.98 

46.24 

3.06 

War  dead (4) 

545 

23.14 

4.31 

173.90 

6.63 

47.29 

2.36 

War  dead (5) 

1265 

— • 

174.39 

6.55 

47.17 

2.30 

Modem  (6) 

113 

37.26 

17.08 

176.07 

8.06 

47.18 

2.54 

Black  Males 

Terry (1) 

203 

47.33 

16.77 

170.18 

7.97 

47.47 

2.93 

Terry ( 2 ) 

360 

49.46 

15.51 

170.23 

7.81 

47.42 

2.97 

Terry (3) 

42 

43.25 

13.21 

172.02 

7.84 

47.92 

3.28 

War  dead (4) 

54 

25.07 

4.98 

172.11 

6.14 

48.34 

2.26 

War  dead (5) 

191 

~ 

173.86 

6.65 

48.41 

2.48 

Modem  (6) 

29 

37.45 

21.52 

179.15 

6.68 

48.24 

2.81 

White  Females 

Terry (1) 

201 

63.19 

15.59 

157.07 

7.76 

42.71 

2.36 

Terry ( 2 ) 

63 

63.93 

16.07 

158.18 

7.51 

42.96 

2.53 

Terry ( 3 ) 

52 

63.35 

17.02 

157.62 

7.96 

42.69 

2.71 

Modern (6) 

89 

37.14 

21.52 

165.54 

5.94 

43.72 

2.02 

Black  Females 

Terry (1) 

199 

48.11 

18.44 

159.04 

6.63 

43.99 

2.41 

Terry ( 2 ) 

177 

47.21 

17.65 

158.39 

6.53 

43.71 

2.39 

Terry (3) 

57 

39.58 

15.52 

159.88 

6.22 

43.96 

2.30 

Modem  (6) 

24 

34.37 

19.73 

165.00 

10.80 

45.52 

3.23 

♦Terry  collection  stature*cadaver  stature-2 . 5cm.  All  others 
are  stature  obtained  during  life. 

1.  Present  study 

2.  Trotter  and  Gleser  (1952) 

3.  Steele  (1970) 

4.  Trotter  and  Gleser  (1952) 

5.  Trotter  and  Gleser  (1958) 

6.  Jantz  and  Moore-Jansen  (1988) 


the  proximal  point  on  the  head  of  the  femur  to  the  midpoint  of 
the  greater  trochanter.  It  is  comparable  to  our  VHA  (proximal 
femur  breadth)  in  the  size  and  location  of  the  fragment. 


Table  3.  Means  and  standard  deviations  (mm.)  of  the  Terry 
collection  segements,  by  race  and  sex. 


W  Males  B  Males  W  Females  B  Females 


Variable 

X 

SD 

X 

SD 

X 

SD 

X 

SD 

VHA 

99.10 

5.87 

98.99 

5.77 

88.24 

5.18 

88.98 

5.24 

VHD 

48.27 

3.17 

47.65 

2.69 

42.54 

2.50 

41.95 

2.35 

VND 

33.09 

3.05 

31.16 

2.60 

29.27 

2.62 

27.17 

2.48 

WSD 

28.26 

2.16 

28.05 

2.21 

25.10 

2.13 

25.18 

2.03 

LCH 

41.35 

2.91 

42.33 

3.00 

36.30 

2.53 

37.05 

2.34 

MCH 

40.92 

3.01 

41.94 

3.22 

37.22 

2.90 

37.75 

2.70 

BCB 

77.81 

4.30 

77.71 

4.57 

68.44 

3.75 

67.74 

3.94 

FDL 

83.42 

4.47 

83.00 

4.22 

74.57 

3.54 

73.94 

3.95 

Table  5  shows  the  correlations  of  our  VHA  (proximal  femur 
breadth)  and  Steele's  segment  1  with  stature  and  maximum  femur 
length.  For  stature,  our  correlations  exceed  Steele's  in  three  of 
the  four  groups;  only  in  White  males  is  it  slightly  lower.  For 
femur  length,  Steele's  correlations  are  higher  in  three  of  the 
four  groups.  All  16  of  the  correlation  coefficients  are  similar, 
indicating  approximately  equal  predictive  efficiency. 

Steele's  segment  4  is  comparable  to  our  MCH  (medial  condyle 
height) .  His  correlations  with  height  and  femur  length  exceed 
ours  in  four  cases  while  ours  exceed  his  in  four  cases.  However, 
our  LCH  (lateral  condyle  height)  yields  consistently  higher 
correlations  than  MCH  in  our  data  (see  table  4) .  Comparing  LCH  to 
Steele's  segment  4  shows  our  correlations  to  exceed  his  in  7  of  8 
comparisons.  It  is  evident  that  lateral  condyle  height  has  the 
higher  correlation  with  stature  and  femur  length  and  is  to  be 
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preferred  as  the  predictor  variable  when  both  condyles  are 
present. 

Table  4.  Correlations  of  femur  segments  with  maximim  femur 
length  and  stature  by  race  and  sex. 


W  Males  B  Males  W  Females  B  Females 


Var. 

Ht 

FML 

Ht 

FML 

Ht 

FML 

HT 

FML 

VHA 

0.587 

0.606 

0.564 

0.592 

0.526 

0.632 

0.432 

0.513 

VHD 

0.462 

0.526 

0.499 

0.454 

0.406 

0.596 

0.540 

0.585 

VND 

0.312 

0.384 

0.393 

0.315 

0.409 

0.409 

0.461 

0.422 

WSD 

0.386 

0.281 

0.251 

0.276 

0.428 

0.295 

0.367 

0.277 

LCH 

0.557 

0.571 

0.503 

0.452 

0.677 

0.665 

0.571 

0.585 

MCH 

0.417 

0.459 

0.391 

0.404 

J.595 

0.518 

0.470 

0.403 

BCB 

0.512 

0.541 

0.509 

0.440 

0.294 

0.445 

0.220 

0.345 

FDL 

0.493 

0.521 

0.560 

0.465 

0.428 

0.537 

0.329 

0.415 

Table  5.  Correlations  of  upper  femur  breadth  (VHA)  and 
Steele's  segment  1  with  height  and  femur  length 


Steele's 
Segment  1 

VHA  (Proxima 
Breadth 

Height 

Femur  Igth 

Height 

Femur  Igth 

White 

males 

0.602 

0.651 

0.588 

0.607 

White 

females 

0.498 

0.623 

0.526 

0.632- 

Black 

males 

0.472 

0.606 

0.564 

0.593 

Black 

females 

0.431 

0.543 

0.432 

0.514 
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Formulae  for  the  Estimation  of  Stature  and  Maximum  Femur  Length 

In  this  section  we  present  the  statistical  constants  for 
estimating  height  and  femur  length  from  various  fragments.  We 
will  present  only  those  equations  which  offer  the  best  predictive 
capability  and  seem  useful  in  forensic  work. 

Table  6  presents  the  slopes,  intercepts  and  standard  errors 
of  estimates  for  each  race/sex  group  for  stature  and  femur  length 
using  the  variable  VHA  (proximal  femur  breadth)  VHD  (femur  head 
diameter)  LCH  (lateral  condyle  height)  MCH  (medial  condyle 
height)  as  predictor  variables.  While  not  presented  as  an 
equation,  these  constants  are  used  in  the  same  manner  as 
traditional  regression  equations.  Thus  to  estimate  stature  for  a 
white  male,  using  VHA  one  takes  the  slope  times  VHA  and  adds  the 
intercept.  The  standard  error  of  estimate  then  provides  an 
indication  of  the  range  within  which  the  true  stature  is  likely 
to  fall.  For  example,  suppose  femur  fragments  from  a  White  male 
yield  the  following  measurements;  VHA=95mm. ;  LCH='40mm. ;  VHD=45mm. 
These  measurements  would  be  used  to  estimate  stature  as  follows: 

Stature=  0.78  x  95mm.  +  89.64  *  163.74cm.  +  6.10cm. 

Stature*  1.47  x  40mm.  +  107.09  *  165.89cm.  +6. 24cm. 

Stature*  1.11  x  45mm.  +  113.89  *  163.84cm.  +  6.77cm. 

The  average  ages  for  all  groups  in  this  study  are  similar  to 
those  presented  by  Trotter  and  Gleser  (1952)  in  their  work  with 
the  Terry  Collection.  In  addition  to  concern  regarding  the 
secular  trend  in  stature  increase  discussed  above,  Trotter  and 
Gleser  (1951a)  addressed  the  issue  of  the  effects  of  aging  on 
stature.  The  correction  of  stature  for  age  related  effects  is 


iq-i  1 


Table  6.  Regression  constants  for  estimating  stature*  and  femur 
length**  from  VHA,  LCH  and  VHD.  These  are  the  three  best 
predictor  variables.  See  text  for  instructions. 


Height 

Femur  length 

slope  intercept  S.E. 

slope  intercept  S.E. 

Predictor  variedsle:  ' 

VHA  (proximal  femur  breadth) 

White  males 

0.78 

89.64 

6.10 

0.25i 

14.81 

2.10 

White  females 

0.73 

91.54 

6.67 

0.21 

21.50 

2.10 

Black  males 

0.79 

91.70 

6.60 

0.32 

13.64 

2.30 

Black  females 

0.59 

107.10 

6.00 

0.25 

21.90 

2.04 

Predictor  variable: 

LCH  (lateral 

condyle  height) 

White  males 

1.47 

107.09 

6.24 

0.54 

20.86 

2.18 

Black  males 

1.34 

113.23 

6.91 

0.46 

25.44 

2.59 

White  females 

1.94 

86.10 

5.77 

0.42 

24.96 

2.09 

Black  females 

1.59 

100.07 

5.47 

0.56 

20.80 

2.03 

Predictor  variable:  VHD  (femur  head  diamter) 


White 

males 

1.11 

113.89 

6.77 

0.43 

23.57 

2.32 

Black 

Males 

1.51 

97.82 

6.92 

0.54 

19.45 

2.56 

White 

females 

1.35 

99.22 

7.16 

0.47 

20.22 

2.06 

Black 

females 

1.59 

92.43 

5.59 

0.58 

17.12 

1.99 

*  Stature*  predicror  variable  x  slope  +  Intercept  +  S.E. 

**  Femur  length*  predictor  variable  x  slope  +  Intercept  +  S.E. 


given  by  Trotter  and  Gleser  (1952);  they  suggest  that  this 
formula  be  applied  to  all  individuals  over  the  age  of  35  for  whom 
stature  is  to  be  estimated. 
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DISCUSSION 


The  method  of  femur  length  and  stattire  estimation 
presented  here  has  several  advantages  over  previous  attempts. 
First,  the  anatomical  landmarks  by  which  these  measurements  are 
defined  are  standard,  well  defined,  and  easy  to  locate.  Many  of 
them  are  already  in  standard  usage  for  the  estimation  of  sex  or 
are  standard  measurements  taken  for  comparative  purposes  on  post 
cranial  material.  Second,  the  standard  errors  of  estimates  were 
equal  to  or  lower  than  those  presented  by  Steele  (1970)  for  his 
smaller  segments.  And  third,  the  sample  size  used  in  this  study 
is  approximately  four  times  that  used  by  Steele  (1970) ,  and  thus 
the  regression  equations  are  more  accurately  estimated. 

Perhaps  the  most  significant  comparison  of  the  technique 
presented  in  this  study  with  that  of  Steele  (1970)  concerns  the 
standard  errors  of  estimate.  For  comparable  segments,  such  as 
those  involving  the  proximal  end  (VHA)  or  the  distal  end  (LCH) 
our  standard  errors  of  estimate  are  uniformly  lower  than 
Steele's.  That  in  turn  indicates  that  our  stature  estimates  are 
more  accurate.  To  some  extent  this  results  from  higher  standard 
deviations  for  femur  length  and  height  in  Steele's  sample.  This 
is  apparently  a  consequence  of  Steele's  smaller  sample  size. 

Our  results  have  demonstrated  the  feasibility  of  estimating 
stature  and  femur  length  from  small  fragments  using  standard 
measurements.  It  is  becoming  increasingly  clear  that  the  Terry 
collection  skeletons  were  drawn  from  a  population  considerably 
smaller  than  contemporary  Americans.  Therefore,  our  results 
should  be  viewed  more  as  indicative  of  the  feasibility  of  the 

technique  than  as  providing  formula  applicable  for  forensic  work 
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on  contemporary  people 


CONCLUSIONS 

The  technique  for  the  estimation  of  stature  from  fragments 
of  the  femur  presented  in  this  report  represents  an  improvement 
over  methods  which  are  currently  in  use.  It  is  desirable  to 
develop  formulae  using  skeletons  from  contemporary  populations, 
if  possible.  It  is  also  possible  that  certain  other  bones  likely 
to  survive  air  crashes  because  of  protection,  such  as  the  bones 
of  the  foot,  calcaneous  and  metacarpals  in  particular,  would  lend 
themselves  to  stature  estimation.  Future  research  should  be 
directed  along  these  lines. 
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In  order  to  more  fully  understand  the  role  of  soil  and 
sedimentary  organic  matter  (SOM)  in  the  migration  and  fate  of 
organic  pollutants  in  aquatic  systems,  a  fast  and  accurate  liquid 
chromatography  method  has  been  developed  for  the  determination  of 
sorption  isotherms.  The  method  is  a  modified  Frontal  Analysis  by 
Characteristic  Point  (FACP)  technique  which  accounts  for 
hydrodynamic  dispersion  and  which  uses  conventional  KPLC 
equipment  coupled  with  commercially  available  chromatography 
software.  We  have  used  this  technique  to  generate  sorption 
isotherms  for  various  fractions  of  the  total  SOM.  Isotherms  are 
then  fit  to  the  Freundlich  equation,  from  which  water  -  organic 
carbon  distribution  coefficients  (Kq^  values)  can  be  calculated. 
In  this  presentation  we  report  on  our  studies  of  the  roles  played 
by  various  fractions  of  soil  organic  matter  (lipids,  humic/ fulvic 
acids  and  humin)  in  the  overall  sorption  process.  By  using  solute 
probes  which  interact  with  solid  phases  in  well-defined  ways,  we 
have  been  able  to  determine  the  contributions  of  acidic,  basic 
and  dipolar  sites  in  the  various  SOM  fractions,  as  well  as  the 
contribution  of  nonspecific  hydrophobic  partitioning,  in  the 
sorption  of  nonionic,  polar  organic  solutes.  It  appears  that 
these  solutes  interact  with  soil  organic  matter  primarily  through 
specific  interactions  with  active  sites  in  the  humin  fraction  of 
the  organic  matrix,  although  there  is  some  hydrophobic 
partitioning  into  the  lipid  fraction  as  well. 


INTRODUCTION 


It  is  now  generally  recognized  that  the  sorption  of 
sparingly  soluble  organic  solutes  onto  soils  and  aquifer  media  is 
dominated  hy  the  natural  organic  carbon  fraction  of  these 
materials.  Lambert  (196'o,  1968, 1978)  and  Lambert  et  al.  (1965) 
first  demonstrated  this  principle  with  studies  of  neutral  organic 
pesticide  upt?ke  onto  a  variety  of  soils.  Lambert  and  co-workers 
(1965)  also  suggested  that  soil  organic  matter  might  act  as  an 
organic  solvent  in  a  solvent  extraction  process,  and  thus  the 
movement  of  neutral  organics  in  soils  could  be  considered 
analogous  to  a  chromatographic  process,  with  dissolved  solutet 
partitioning  between  a  polar  aqueous  mobile  phase  and  a  non-polar 
stationary  phase  of  soil  organic  matter  coating  an  inorganic 
matrix.  The  partitioning  of  a  neutral  organic  solute  in  such  a 
system  should  then  be  correlated  with  the  solute's  partitioning 
between  water  and  an  immiscible  organic  solvent,  and  Briggs 
(1973)  and  Chiou  et  al.  (1979)  developed  empirical  equations 
which  related  sorption  in  soil  systems  to  a  solute's  octanol- 
water  partition  coefficient. 

The  importance  of  sorption,  including  adsorption  and/or 
partitioning  on/ into  soil  organic  matter,  in  the  transport  of 
nonionic  organic  solutes  can  be  addressed  by  first  considering 
the  relative  residence  time,  or  retardation  factor,  (tj-)j^,  of  a 
solute  as  it  travels  from  source  to  observation  point  in  an 
aquifer.  This  retardation  factor  is  defined  by  equation  1: 


where  t* and  represent  the  average  travel  times  of  solute  i 
and  water  molecules,  respectively.  The  effect  of  sorption  on  this 
factor  can  be  quantitatively  estimated  through  equation  2: 

(tr)i  «  1  +  (Kd)i*  (2) 

aq 

where  (K^j)  is  the  distribution  coefficient  for  solute  i  between 
water  and  the  aquifer  media  and  and  gg  represent  the  average 
•bulk  density  and  effective  porosity  of  the  media,  respectively. 

Distribution  coefficients  for  the  sorption  of  nonionic 
solutes  on  complex,  heterogeneous  soils  and  aquifer  materials 
actually  represent  the  weighted  sum  of  all  individual 
contributions: 

^^d^i  “  ^j(^d)i,j  “  ^inorg*^inorg  ^org*^org 

where  fj  represents  the  fraction  of  phase  j  and  j  the 
distribution  coefficient  for  solute  i  on  that  phase.  However, 
Karickhoff  et  al.  (1979)  evaluated  the  sorption  of  many 
hydrophobic  organic  solutes  on  a  variety  of  organic  rich  (>l% 
organic  matter)  natural  sediments  and  found  that  sorptive  effects 
cf  the  inorganic  matrix  were  negligible.  Equation  3  could 
therefore  be  written 


(Kd)i 


foc*(K 


OC>  1 


:4) 


where  (K^^,)  is  the  distribution  coefficient  for  the  sorption  of 


i  onto  the  organic  fraction  alone  and  is  the  mass  percentage 
of  this  fraction  in  the  mineral  assemblage.  Karickhoff  et  al. 
(1979)  also  evaluated  the  relationship  between  and  a  variety 
of  solute  properties  and  found  first  that,  for  sorption  onto  the 
fine  fractions  (<50  um  particles)  of  the  sediments,  a  solute 
tended  to  have  a  reasonably  constant  regardless'  of  the 
origin  of  the  organic  matter.  These  workers  also  observed,  in 
agreement  with  previous  studies,  a  highly  significant  correlation 
between  and  a  solute's  octanol/water  partition  coefficient, 

^ow* 


log  Kqc  =  a  *  log(KQ^,)  +  b  (5) 

Sorption  onto  the  sand  fraction  (particles  >50  um)  could  also  be 
described  by  equation  5,  but  the  fitting  coefficients  were 
different  than  those  for  the  fine  fraction. 

Substituting  equations  4  and  5  into  equation  2  yields 

(tj.)i  =  1  +  C*foc*(Kov)i*  +  C'  (6) 

aq 

The  relative  residence  time  of  organic  solutes  in  aquifers  can 
therefore  be  predicted  from  equation  6  if  accurate  measurements 
or  estimates  of  fQ^,  (KQy)j^,  gg  and  gg  are  available. 

Numerous  studies  have  confirmed  the  conclusions  of  Lambert 
and  co-workers,  Briggs,  Chiou  et  al.  and  Karickhoff  et  al. 
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(Hamaker  and  Thompson,  1972;  Schwartzenbach  and  Westall,  1981; 
Chiou  et  al.,  1983;  Chiou  et  al,  1985;  Karickhoff,  1984).  Tnese 
studies  invariably  involve  hydrophobic  organic  solutes ,  such  as 
halogenated  hydrocarbons,  and  organic  rich  soils  and  sediments. 
For  such  systems  it  is  generally  accepted  that  interactions 
between  solutes  and  soil  organic  matter  result  from  hydrophobic 
effects  (Tanford,  1980) .  The  driving  force  for  such  'hydrophobic 
interactions'  is  the  increase  in  entropy  associated  with  the 
removal  of  the  structured  hydration  region  that  surrounds  a 
solute  in  water.  This  hydrophobic  model  also  has  been  shown  to  be 
valid  when  organic  cosolvents  are  present,  and  therefore  is  more 
generally  termed  a  'solvophobic'  model  (Banerjee  et  al.,  1980; 
Rao  et  al.,  1985;  Nkedi-Kizza  et  al.,  1985). 

In  spite  of  these  advances  in  understanding  the 
thermodynamic  driving  forces  that  control  the  adsorption  and/or 
partitioning  of  organic  solutes  between  water  and  soils,  there 
are  still  many  questions  which  have  not  been  answered.  Even  for 
hydrophobic  solutes  which  should  interact  only  through 
hydrophobic  forces,  values  have  been  observed  to  vary  by 
as  much  as  factors  of  3  to  5  for  different  roineral/organic 
assemblages  (Curtis  et  al.,  1986).  Mingelgrin  et  al.  (1983)  found 
literature  values  of  Kqjjj  (K^j  expressed  on  an  organic  matter 
rather  than  organic  carbon  basis)  that  varied  by  more  than  an 
order  of  magnitude.  This  is  an  indication  that  the  fitting 
coefficients  of  equation  5  could  be  a  function  of  the  nature  of 
the  soil  organic  matter  involved.  It  is  well  known  in  the  field 
of  organic  geochemistry  that  the  chemical  character  of  organic 
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matter  in  sediments  depends  both  on  the  paleoecological 
conditions  at  the  time  of  deposition  and  the  extent  of  diagenesis 
within  the  sediment  (Tissot  and  Welte,  1984).  It  should  not  be 
surprising,  therefore,  to  find  that  soil  and  sedimentary  organic 
matter  might  be  highly  variable  in  its  ability  to  interact  with 
organic  solutes.  Indeed,  Lambert  et  al.  (1965)  suggested  this  in 
their  original  work.  Hartley  (1960)  argued  that  the  fats-waxes- 
resins  fraction  of  soil  organic  matter  was  actually  responsible 
for  the  sorption  of  nonionic  organic  solutes.  In  probably  the 
most  definitive  study  to  date  on  this  subject,  Garbarini  and  Lion 
(1986)  showed  that  different  components  of  soil  organic  matter 
had  widely  varying  affinities  for  toluene  and  trichloroethylene 
(TCE) ,  both  nonionic,  hydrophobic  organic  solutes. 

As  noted  previously,  virtually  all  the  studies  cited 
involved  nonpolar,  hydrophobic  organic  solutes.  If  soil  organic 
matter  does  vary  in  its  affinity  for  organic  solutes,  this 
variability  would  become  more  pronounced  as  the  water  solubility 
of  the  organic  solute  increases.  With  increasing  water  solubility 
the  non-specific  hydrophobic  effect  becomes  less  important  in  the 
sorption  process,  while  specific  solute-organic  matter 
interactions  correspondingly  increase  in  importance.  In  such 
cases  it  would  therefore  be  necessary  to  consider  the  specific 
adsorptive  characteristics  of  soil  and  aquifer  organic  matter. 

In  this  and  the  following  paper  we  report  the  results  of  our 
initial  studies  of  the  sorptive  properties  and  chemical 
characteristics  of  several  fractions  of  soil  organic  matter  using 
the  'Inverse  Chromatography'  technique.  In  this  approach  the 
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column  material  itself  is  characterized  by  observing  the 
retention  behavior  of  solutes  which  are  chosen  because  they  are 
known  to  interact  in  chemically  well-defined  ways.  Tne  explosive 
growth  in  more  conventional  chromatographic  separations  thar  has 
occurred  over  the  past  15  years  has  overshadowed  the  rather 
unique  ability  of  chromatographic  methods  to  measure  meaningful 
physiochemical  parameters. 

In  order  to  ensure  that  the  organic  matter  under 
investigation  is  indeed  variable  in  both  composition  and 
adsorptive  character,  we  have  chosen  to  study  fractions  of  a 
single  organic  matter  sample  rather  than  a  variety  of  different 
samples.  An  organic  rich  sediment  (peat)  was  first  fractionated 
into  standard  organic  geochemical  classes:  lipids,  humic  and 
fulvic  acids,  and  non-extractable  humin.  By  comparing  the 
sorption  of  the  chosen  solutes  on  these  various  organic  matter 
fractions  with  the  sorption  observed  with  the  total  organic 
matter,  we  have  been  able  to  define  the  relative  importance  of 
each  of  these  fractions  in  determining  the  overall,  specific 
chemical  character  of  soil  organic  matter. 

In  this  first  paper  we  report  on  our  studies  of  the 
interactions  of  polar  organic  solutes  with,  soil  organic  matter 
fractions  using  a  dispersion  corrected  High  Performance  Liquid 
Chromatography  (HPLC)/  Frontal  Analysis  by  Characteristic  Point 
(FACP)  technique.  We  will  also  present  the  results  of  structural 
and  compositional  studies  carried  out  primarily  by  solid  srate 
C  NMR  spectroscopy.  This  is  a  powerful  technique  for  studying 


insoluble  organic  material,  but  a  technique  which  has  not  yet 
found  significant  applications  in  environmental  chemistry. 
Analyses  of  the  solid  state  NMR  spectra  have  allowed  us  to 
correlate  soil  organic  matter  structure  with  sorption  properties. 
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EXPERIMENTAL 


Substrate  Preparation 

Commercially  available  peat  was  prepared  for  fractionation 
by  lyophilizing  for  24  hours  to  obtain  what  will  be  termed  the 
total  soil  organic  matter  (SOM)  fraction.  Non-bound  lipid 
material  was  then  removed  by  repeated  methylene  chloride-methanol 
(2:1)  extractions,  followed  by  vacuum  filtration  until  a  clear 
organic  phase  was  achieved.  The  resulting  lipid  material  was 
rotary-evaporated  to  dryness.  The  residue  (SOM-Lipid)  was  split 
and  one-half  subjected  to  a  48  hour  saponification  in  IN  NaOK, 
followed  by  repeated  washings  with  O.OIN  NaOH  until  the 
supernatant  became  clear.  This  procedure  removed  the  humic  and 
fulvic  acids.  The  non-extractable  residue  was  subsequently  washed 
with  methylene  chloride  to  remove  bound  lipids  freed  during  the 
saponification.  The  remaining  residue  is  termed  the  Humin 
fraction.  The  entire  fractionation  scheme  is  summarized  in  Figure 


Column  Preparation 

The  three  insoluble  residue  fractions  were  each  dried  by 
lyophilization,  frozen  in  liquid  nitrogen  and  ground  in  a  mortar 
and  pestle.  The  material  was  then  dry-sieved  to  isolate  a  38-53 
micron  size  fraction.  To  ensure  reasonable  chromatographv 
retention  times  an  appropriate  amount  of  organic  matter  was 
weighed  and  mechanically  mixed  with  37-44  micron  pellicular 
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silica  gel  to  yield  a  5%  w/w  organic  matter  mixture.  Lipid 
material  was  solublized  in  methylene  chloride  and  rotary- 
evaporated  onto  the  silica  gel  support  to  yield  a  1%  w/w  coated 
packing  material.  The  prepared  materials  were  dry  packed  into 
stainless  steel  HPLC  analytical  columns  (4.1  mm  I.D.)  fitted  with 
5  micron  column  frits.  Prior  to  analysis  each  column  was  allowed 
to  equilibrate  for  4-6  hours  under  continuous  flow  of  deionized 
water.  After  analysis  packing  material  was  removed  from  the 
columns,  dried,  and  percent  organic  carbon  determined  using  a 
Horiba  model  PIR  2000  General  Purpose  Infrared  Gas  Analyzer. 

Materials  and  Reagents 

Solute  probes  were  purchased  from  Mallinkrodt  (St.  Louis, 
Mo.)  and  were  of  analytical  reagent  grade.  Aniline  was  purified 
by  distillation  prior  to  use.  Vydac  pellicular  silica  gel  (37-44 
micron)  was  purchased  from  Alltech  Associates  Inc.  (Deerfield, 
Ill.).  Water  used  in  the  analysis  was  purified  by  distillation 
and  deionized  using  a  Millipore  Milli-Q  system. 

NMR  Spectroscopy 

CP-MAS  imR  spectra  were  obtained  with  a  Eruker  WP-200  SY 
spectrometer  equipped  with  an  IBM  solids  accessory  operating  at 
50  MHz  for  carbon.  The  spectrometer  was  tuned  with  t-butyl 
benzene  and  the  chemical  shifts  are  reported  as  parts  per  million 
downfield  from  tetramethylsilane  .  Lyophilized  samples  were  spun 
at  the  magic  angle  at  4  KHz.  Use  and  optimization  of  this 
instrument  for  studying  sedimentary  organic  carbon  has  been 
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described  previously  (Cooper  et  al. ,1985;1986) .  Dipolar  dephasing 
delay  spectra  were  obtained  using  a  contact  time  of  i  millisecond 
and  a  50  microsecond  delay  period  prior  to  data  acquisition 
during  which  the  proton  decoupler  is  turned  off.  A  180° 
refocusing  pulse  was  inserted  at  the  mid-point  of  this  delay 
period. 

Frontal  chromatograms  used  in  the  dispersion  corrected  FACP 
technique  were  generated  using  a  conventional  HPLC  syster 
consisting  of  a  Valeo  six-port  valve  outfitted  with  a  13  ml 
sample  loop,  a  Waters  Model  510  Solvent  Delivery  System,  and  a 
Waters  Model  440  Absorbance  Detector  operating  at  254  nm 
wavelength.  All  data  was  collected  by  and  processed  using  a 
Hewlett-Packard  HP-1000  Lab  Automation  System. 

Isotherm  Generation 

Isotherms  were  calculated  from  the  rear  diffuse  boundary 
(i.e.  the  desorption  branch)  of  a  single  frontal  chromatographic 
peak.  This  technique  is  termed  Elution  by  Characteristic  Point 
(ECP) ,  or  Frontal  Analysis  by  Characteristic  Point  (FACP) .  It  was 
first  introduced  by  Cremer  (1961)  and  Huber  and  Kuelemans  (1962), 
and  has  been  thoroughly  discussed  by  Conder  and  Young  (1979). 
With  this  technique  each  point  on  the  isotherm  is  calculated 
according  to  equation  7. 
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c 


(7) 


q(C)  =  *  (V-Vo)*dC 

m  Q 

where  q(C)  is  the  uptake,  or  the  concentration  of  adsorbed  solute 
in  equilibrium  with  the  aqueous  mobile  phase  concentration  C,  m 
is  the  mass  of  organic  matter  present  in  the  column,  V  is  the 
volume  eluted  at  concentration  C,  and  Vq  is  the  void  volume 
contribution.  Void  volumes  were  determined  from  frontal 
chromatograms  of  non-retained  solutes. 

Ti.e  FACP  technique,  in  contrast  to  the  more  familiar  Frontal 
Analysis  (FA)  method,  does  not  take  into  account  dispersional 
broadening  that  occurs  in  the  column.  In  order  to  refine  the  FACP 
technique  we  have  introduced  a  dispersion  correction  which 
employs  the  diffuse  rear  boundary  of  a  non-retained  solute 
breakthrough  curve.  The  advantage  of  this  approach  is  that  the 
entire,  dispersion  corrected  isotherm  can  be  determined  from  the 
diffuse  rear  boundary  of  a  single  chromatogram,  whereas  the  FA 
method  requires  a  separate  chromatogram  for  each  isotherm  point. 
Solutes  were  determined  to  be  non-retained  by  comparing  their 
retention  volumes  to  those  of  calcium  chloride,  as  well  as  column 
void  volumes  determined  gravimetrically  (Jacobson  et  al.,  1984). 
These  comparisons  agreed  within  3%. 
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RESULTS 


Reliability  q1  tM  HP.LC/fMg  M§thod 

The  dispersion  corrected  HPLC/FACP  method  was  compared  with 
both  a  frontal  chromatography  technique  and  the  more  conventional 
batch  isotherm  method.  Kq^,  values  calculated  by  the  uncorrected 
FACP  method  differed  from  frontal  Kqj,  values  by  up  to  14%,  but 
when  the  dispersion  correction  was  applied  the  difference  was 
reduced  to  within  2%.  When  compared  to  K^^'s  obtained  by  the 
batch  isotherm  method,  dispersion  corrected  HPLC/FACP  values 
differed  by  as  little  as  1%  to  a  maximum  of  11%.  Comparisons  of 
Kqj,  values  generated  by  these  various  methods  are  summarized  in 
Table  1,  where  the  worst  cases  observed  are  included. 

The  dispersion  corrected  HPLC/FACP  method  is  tnerefore  seen 
to  produce  results  which  are  in  substantial  agreement  with  those 
obtained  by  other,  more  conventional  techniques.  It  is  important 
to  note,  however,  that  an  entire  isotherm  can  be  generated  from  a 
single  chromatogram  by  this  method.  The  frontal  chromatography 
technique  requires  separate  chromatograms  for  each  data  point  on 
the  isotherm,  while  the  batch  method  involves  an  individual 
experiment  for  each  data  point  (extraction,  analysis,  etc.).  Once 
a  column  has  been  prepared,  the  HPLC/FACP  method  can  generate  an 
isotherm  approximately  every  30  minutes.  The  speed  of  this 
technique  thus  allows  a  wide  spectrum  of  solute  probes  to  be 
investigated  in  a  relatively  short  time  and  with  a  minimum  of 
sample  required. 


Soil  Constituent  Trends 


Isochenn  data  for  the  four  isolated  soil  organic  matter 
fractions  were  generated  using  four  solute  'probes'.  These  solute 
probes  were  chosen  because  they  are  known  to  interact  with 
surfaces  in  well  defined  ways  (Snyder,  1974,  1978).  Cresol  ic  a 
'soft'  acid  and  interacts  strongly  with  basic  sites;  aniline  a 
'soft'  base  that  interacts  with  acidic  sites;  nitrobenzene 
interacts  strongly  with  dipolar  sites;  and  benzene  is  non- 
functionalized  and  therefore  interacts  via  non-specific 
dispersion  forces.  While  both  cresol  and  aniline  are  weakly 
ionizable,  at  neutral  pH  their  dissociation  constants  are  such 
that  they  exist  in  essentially  unionized  form  (pKgj[ cresol]  = 
10.2,  pKjj[aniline]  =  9.3). 

Sorption  data  for  the  four  probes  and  four  organic  matter 
fractions  are  summarized  in  Figure  2,  where  experimental 
values  are  plotted  vs.  calculated  K^^'s.  The  fitting  coefficients 
(a  and  b  of  equation  5)  are  summarized  in  Table  2.  Interestingly, 
all  fractions  exhibit  a  reasonably  linear  relationship  between 
Kqj,  and  Kq^,  as  predicted  by  equation  5,  although  in  some  cases 
this  linear  relationship  holds  for  only  three  of  the  four 
solutes.  Surprisingly,  however,  it  is  only  the  lipid  fraction 
that  exhibits  the  positive  slope  observed  by  others  when  using 
nonpolar,  hydrophobic  solutes.  Soil  Organic  Matter  (SOM),  SOM- 
Lipids,  and  Humin  all  show  a  negative  relationship  between 
and  Kq^.  We  interpret  these  results  as  indicating  that  only 
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sorption  onto  the  lipid  fraction  of  SOM  is  driven  by  non¬ 
specific  hydrophobic  (or  solvophobic)  partitioning.  Sorption  onto 
the  other  fractions  appears  to  result  from  specific  interactions 
between  polar  sites  in  the  organic  matter  and  the  polar  solutes. 
Since  polarity  generally  increases  with  increasing  water 
solubility,  and  therefore  decreases  with  increasing  octanol-water 
partition  coefficients,  such  specific  interactions  would  then  be 
expected  to  decrease  as  values  increase.  This,  of  course,  is 
precisely  the  behavior  we  have  observed  with  the  insoluble  soil 
organic  matter  fractions. 

Consideration  of  the  chemical  nature  of  the  solute  probes 
used  in  these  studies  gives  insight  into  the  observed  behavior. 
Aniline  (K^y  *  0.98),  the  soft  base  probe,  exhibits  the  highest 
Kqp  values  on  all  three  insoluble  SOM  fractions,  indicating  a 
relatively  large  population  of  acidic  sites  in  these  fractions. 
These  sites  are  probably  due  to  carboxyl  and  phenolic  groups,  a 
conclusion  which  is  confirmed  by  the  NMR  spectra  (Figure  3). 
Nitrobenzene  (Kgy  »  1.85)  shows  the  next  highest  values, 

indicating  numerous  polar  but  non  acidic  sites.  Again,  NMR 

is  helpful  in  interpreting  the  nitrobenzene  data,  as  the  spectra 
show  substantial  methoxy  (carbohydrate)  structure.  Benzene  (K.,. ,  = 
2.13),  the  non-specific  probe,  shows  the  smallest  affinity  for 
the  insoluble  fractions,  again  suggesting  that  these  materials 
are  highly  polar  and  functionalized.  The  soft  acid  probe  cresol 
(^ow  *  1*54)  exhibits  erratic  behavior  which  we  attribute  to 
either  strong  hydrogen  bonding  with  the  aqueous  solvent 
exclusion  from  the  largely  acidic  organic  matrix. 


or 


It  is  interesting  to  note  from  the  data  of  Table  2  that 
removal  of  the  lipid  material  from  the  intact  SOM  increases  the 
slope  and  intercept  of  the  log  vs.  log  plot.  This 

indicates  that  binding  of  polar  solutes,  when  normalized  for  the 
fraction  of  orgonrc  carbon  comprising  the  sorbent,  increases  once 
the  lipid  fraction  is  removed.  This  increased  sorptivity  may  be  a 
result  of  the  lipid  components  hindering  access  to  polar  sites  in 
the  remaining  organic  matter  matrix,  or  a  concentration  of  polar 
sites  in  the  residue.  A  similar  effect  is  observed  vnen  the  humic 
and  fulvic  acids  are  removed,  with  the  humin  fraction  exhibiting 
the  largest  sorptive  capacity  for  polar  solutes  per  mass  of 
organic  carbon  than  any  of  the  fractions  studied.  We  interpret 
this  trend  in  increasing  sorptivity  with  removal  of  extractable 
organic  matter  as  an  indication  that  the  polar,  active  sites  in 
soil  organic  matter  are  concentrated  in  the  non-extractable  humin 
fraction. 

Relationship  Between  Sorptive  Characteristics  and  Structure  of 
Soil  Organic  Matter 

With  the  development  of  high  power  decoupling,  cross 
polarization  and  magic  angle  spinning  techniques  (Pines  et  al., 
1973;  Schaefer  and  Stejskal,  1976),  it  is  now  possible  to  obtain 
NMR  spectra  of  solid  ,  insoluble  materials  that,  in  favorable 
situations,  approach  the  quality  of  solution  spectra.  However, 
the  inherent  heterogeneity  of  soil  and  sedimentary  organic  matter 
results  in  peaks  which  are  sufficiently  broad  so  as  to  preclude 
any  realistic  interpretation  of  the  'fine  structure'  that  may  be 
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observed.  For  that  reason  solid  state  NMR  spectra  of  SOM  are 
normally  analyzed  by  comparing  relative  intensities  in  four 
spectral  'windows':  0-50  ppm  (aliphatic  carbon),  50-100  ppm 
(methoxy  [carbohydrate]  carbon),  100-160  ppm  (olefinic  and 
aromatic  carbon) ,  and  160-190  ppm  (carboxyl  and  amide  carbon) . 

NMR  spectra  of  SOM,  SOM-Lipid,  and  Humin  are  displayed  in 
Figure  3.  One  of  the  most  striking  features  of  these  spectra  is 
the  pronounced  relative  decrease  in  the  aliphatic  peak  (ca.  30 
ppm)  when  the  lipids  are  removed  (i.e.  Figure  3(b)  compared  to 
3(a)).  This  indicates  that  non-functionalized,  hydrophobic 
aliphatic  carbon  is  the  chief  component  of  the  lipid  fraction. 
This  explains  the  sorption  results  for  both  the  isolated  lipid 
fraction  (hydrophobic  partitioning)  ,  as  well  as  the  differences 
between  the  SOM  and  SOM-Lipid  fractions.  While  both  of  these 
fractions  appear  to  interact  with  the  polar  solute  probes  through 
specific  interactions  at  polar  sites  in  the  organic  matter 
matrix,  the  relative  number  of  sites  in  the  SOM-Lipid  fraction 
appears  to  be  greater  than  the  relative  number  of  sites  in  the 
unextracted  SOM.  In  other  words,  the  lipid  fraction  appears  to 
dilute  the  concentration  of  active  polar  sites  in  the  total  SOM 
fraction.  The  lipid  fraction  might  also  physically  and/or 
chemically  block  some  active  sites  from  interacting  with  the 
solute  probes. 

Removal  of  humic  and  fulvic  acids  (Humin,  Figure  3(c))  does 
little  to  the  overall  structure  of  the  organi'^  matter,  except  to 
reduce  the  relative  intensity  of  the  carboxyl  peak  (ca.  170-180 


ppm).  This  is  an  entirely  expected  result,  of  course,  since 
humics  and  fulvics  are  thought  to  be  primarily  polycarboxylate 
type  materials.  The  important  feature  of  the  Humin  spectrum  is 
that  peaks  arising  from  the  important  polar  carbon  functional 
groups  persist  after  extraction  of  lipids  and  humic/fulvic  acids, 
and  the  relative  intensities  of  these  peaks  are  approximately  the 
same  as  those  observed  in  the  spectrum  of  unextracted  SOM  (4(a)). 
These  persistent  polar  carbon  functional  groups  include  methoxy 
(70-80  ppm),  olefinic  (100-120  ppm),  aromatic  (120-140  ppm), 
substituted  aromatic  (primarily  phenolic,  140-160  ppm) ,  and  to  a 
lesser  extent  carboxyl/amide  (170-190  ppm) .  Also  of  interest  in 
the  Humin  spectrum  is  the  intense  aliphatic  peak  at  ca.  40  ppm, 
which  is  slightly  downfield  from  the  aliphatic  peaks  in  the  other 
two  spectra.  We  attribute  this  peak  to  -CH2-  groups  adjacent  to 
methoxy  carbon,  suggesting  that  a  large  percentage  of  aliphatic 
carbon  in  the  Humin  fraction  is  tied  up  in  polar  carbohydrate 
material,  thus  reducing  the  hydrophobic  character  of  this 
aliphatic  carbon.  Interpretation  of  the  NMR  spectra  thus  explains 
nicely  the  increased  polar  sorptivity  of  the  Humin  fraction 
relative  to  SOM  and  SOM-Lipid. 

Figure  3  also  includes  the  50  usee  dipolar  dephased  spectra 
of  each  SOM  fraction.  Dipolar  dephasing  (DD)  exploits  the  time 
dependence  of  .magnetic  relaxation  to  further  enhance  spectral 
resolution,  and  this  technique  represents  what  has  been  termed 
the  'second  generation'  of  solid  state  NMR  technology  (Maciel, 
1984).  DD-CP-MAS  NMR  is  finding  utilization  in  a  variety  of 
studies  of  organic  matter  in  soils  (Wilson  et  al.,  1983),  resins 
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(Wilson  et  al.,  1984),  anoxic  marine  sediments  (Cooper  et  al., 
1985*  1986),  and  peats  (Orem  and  Hatcher,  1987). 

To  a  first  approximation,  relaxation  of  C  nuclei  in  the 
dipolar  dephasing  experiment  is  directly  reT.ated  to  the  number  of 
bonded  protons.  Thus,  intensities  arising  from  -•CH2-  groups  decay 
faster  than  those  from  -CH-  groups,  'which  will  themselves  decay 
faster  than  quaternary  aliphatic  carbon.  Similarly,  unsubstituted 
aromatic  carbon  nuclei  decay  faster  than  substituted  aromatic 
carbon.  The  exception  to  this  general  .trend  are  “CH3  groups, 
which  have  long  relaxation  times  due  to  their  rotational 
mobility.  Our  previous  work  (Cooper  et  al.,  1985;1986)  and  that 
of  others  indicate  that  spectra  obtained  with  a  50  usee  delay 
show  resonances  arising  only  from  non-protonated  ^^C  nuclei,  as 
well  as  from  "CH3  groups. 

Comparisons  of  the  50  usee  DD  spectra  of  the  various  SOM 
fractions  support  previous  conclusions  regarding  the  sorptivity 
and  structure  of  these  materials.  Removal  of  the  lipids  from  the 
SOM  has  little  effect  on  the  DD  spectra  (Figure  3(b)  compared  to 
Figure  3(a)).  Distinct  resonances  in  the  DD  spectrum  of  humin 
(Figure  3(c)  are  for  the  most  part  obscured  due  to  the  low  signal 
to  noise  ratio  that  results  from  the  lower  carbon  content  of  this 
fraction.  However,  the  broad  peak  between  100  and  150  ppm,  and 
the  carboxyl  peak  at  175  ppm,  indicate  substantial  non- 
protonated,  functionalized  olefenic,  aromatic  and  carboxyl/amide 
carbon  in  the  humin.  The  dipolar  dephasing  NMR  experiment  thus 
supports  the  conclusion  that  active,  polar  sorption  sites  in  SOM 
are  located  in  the  non-extractable  humin  fraction. 
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CONCLUSIONS 


We  have  shown  that  polar,  nonionic  organic  solutes  interact 
with  soil  organic  matter  primarily  through  specific  interactions 
with  active  sites  in  the  humin  fraction  of  the  organic  matrix, 
and  not  as  a  result  of  the  hydrophobic  forces  that  are 
responsible  for  the  sorption  of  nonpolar  organic  solutes.  It  is 
only  sorption  onto  the  lipid  component  of  soil  organic  matter 
that  appears  to  be  a  true  hydrophobic  phenomenon.  Our  results 
therefore  substantiate  the  arguments  of  Hartley  (i960),  who 
described  the  fats-waxes-resins  fraction  of  soils  (i.e.  soil 
lipids)  as  a  solvent  for  organic  pesticides.  Our  results  are  also 
in  substantial  agreement  with  those  of  Garbarini  and  Lion  (1986), 
who  observed  wide  variations  in  values  for  nonpolar, 

hydrophobic  solutes  with  different  soil  organic  matter  fractions. 
Interestingly,  Garbarini  and  Lion  also  observed  the  highest 
values  with  nonextractable  humin. 

We  believe  that  these  results  and  the  results  of  others 
cited  have  important  implications  for  the  modeling  of  organic 
solute  migration  in  ground  water.  For  example,  with  regard  to 
hydrophobic  solutes,  it  may  be  necessary  to  recast  equation  6  in 
terms  of  the  weight  fraction  of  lipids  in  the  soil  matrix,  not 
the  total  organic  matter  fraction.  On  a  broader  scale,  however, 
it  is  probably  unreasonable  to  expect  a  single  empirical 
parameter  (e.g.  K^^,)  to  adequately  predict  sorption  in  soils, 
sediments  and  aquifers,  given  the  variations  in  organic  matter 
composition  which  are  to  be  expected,  as  well  as  our  data  and 
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those  of  others  which  indicate  that  the  various  organic  matter 
fractions  differ  substantially  in  their  sorptive  cnaracteristics. 
Qualitative  interpretation  of  the  conventional  and  dipolar 
dephased  NMR  spectra  indicate  that  polar  sorption  sites  in 

the  humin  fraction  arise  from  aromatic,  carboxyl/amide  and 
carbohydrate  carbon  nuclei,  many  of  which  are  non-protonated. 
These  qualitative  observations  are  in  agreement  with  the 
statistical  analyses  of  Garbarini  and  Lion  (1986;  ,  who  concluded 
that  elemental  C,K  and  0  data  are  necessary  parameters  for 
accurately  describing  solute-soil  organic  matter  interactions. 

We  also  believe  that  much  more  work  is  necessary  before 
accurate,  a  priori  estimates  of  polar  solute  sorption  can  be 
made.  Current  efforts  underway  in  our  laboratory  are  directed 
toward  the  development  of  a  more  general  model  of  sorption  of 
polar  solutes  based  on  Snyder's  concept  of  liquid  chromatographic 
selectivity  (1974,1978). 
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Table  1.  Comparison  of  Kj^  Values  Obtained  by  Different 
Techniques. 


System 


Aniline/Humin  o-Cresol/Humin 


Uncorrected  HPLC/FACP 

Corrected  HPLC/Facp  - 

Frontal  Analysis  - 

Corrected  HPLC/FACP  (dynamic)  40 

Conventional  Batch  (static)  39 


10. 

9. 

9.3 


Table  2.  Linear  Fitting  Coefficients  fro  Equation  5. 


S-vmfas.e 

1 

h 

correlation  coefficient 

Lipid* 

0.31 

1.24 

0.95 

SOM 

-0.36 

1.97 

0.80 

Som-Lipid 

-0.44 

2.53 

0.92 

Humin 

-0.61 

3.16 

0.92 

*  o-cresol 

was  excluded 
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FIGURE  CAPTIONS 


Figure  1.  Extraction  scheme  used  to  isolate  soil  organic  matter 
fractions. 

Figure  2.  Plot  of  log  versus  log  for  the  isolated  organic 
matter  fractions. 

Figure  3.  CP-MAS  NMR  spectra  of  (a)  total  soil  organic  matter 

/ 

(SOM)  (b)  SOM-Lipids  (c)  Humin.  Top  spectra  represent  0 
usee  DD  delay,  while  bottom  represent  50  usee  delay. 
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Gas  chromatography  has  been  used  to  measure  the  molal 
activity  coefficients  of  a  series  of  nonpolar  organic  solutes  in 
sedimentary  organic  matter  and  its  various  extractable 
components.  Organic  matter  from  a  peat  deposit  of  higher  plant 
origin  in  the  Florida  Everglades  was  isolated  at  three  depths  in 
the  peat  core.  The  organic  matter  was  then  fractionated  into  its 
lipid,  humic  acid  and  humin  components,  and  the  activity 
coefficients  measured  for  solutes  ranging  in  electronic 
polarizability  from  pentane  to  dichloroethane.  Activity 
coefficients,  which  are  a  thermodynamic  measure  of  the  affinity 
of  a  solute  for  a  polymeric  solution,  show  a  reasonable 
correlation  with  electronic  polarizability  for  the  total  organic 
matter  and  its  lipid  and  humin  fractions.  There  is  also  a  marked 
increase  in  the  lipophilic  character  of  the  humin  fraction  with 
increasing  depth  of  burial.  In  spite  of  this  change  in  humin 
solution  capabilities,  however,  all  of  the  f^olutes  exhibit 
decreased  affinities  for  the  deeper  total  organic  matter,  which 
we  believe  is  due  to  its  lower  lipid  content,  as  well  as  a 
somewhat  lower  affinity  for  the  more  mature  lipid  fraction.  Solid 
state  I'll®  spectra  of  the  unfractionated  organic  matter  and 
the  humic  acid  and  humin  fractions  give  some  insight  into  the 
relationship  between  their  structure  and  sorption 
characteristics,  although  the  chemical  basis  for  the  evolution  of 
the  nonpolar,  lipophilic  humin  at  depth  is  not  readily  apparent. 
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INTRODDCTIOK 


In  the  preceding  manuscript  (Antworth  et  al.,  1987)  we 
demonstrated  that  the  sorption  of  polar,  nonionic  organic  solutes 
onto  fractions  of  soil  organic  matter  depended  on  the  composition 
of  the  organic  matter  fraction.  The  results  presented  in  that 
paper  suggest  that  hydrophobic  (solvophobic)  partitioning  of 
these  solutes  occurs  only  with  the  lipid  fraction,  and  that 
specific  polar  interactions  between  the  solutes  ano  active  sites 
in  the  nonextractable  humin  fraction  actually  dominate  the 
sorption  of  polar  solutes.  It  would  thus  appear  that  the 
empirical  relationship  between  octanol-water  partition 
coefficients  and  soil  and  sedimentary  organic  matter  partitioning 
(Briggs,  1973;  Chiou  et  al.,  1979;  Karickhoff  et  al.,  1979) 
cannot  be  extended  to  polar  organic  solutes  with  water 
solubilities  greater  than  that  of  benzene  (-1000  ppm) .  It  would 
also  appear  that  the  composition  and  structure  of  the  organic 
matter  involved  must  also  be  considered.  Indeed,  the  work  of 
Garbarini  and  Lion  (1986),  who  correlated  the  sorption  of  toluene 
and  trichloroethylene  into  a  variety  of  organic  matter  samples 
with  C  and  0  content  of  the  organic  matter,  suggest  that  soil 
organic  matter  composition  must  be  considered  for  accurately 
predicting  sorption  of  even  hydrophobic  solutes  (solubilitites 
less  than  -100  ppm) .  Garbarini  and  Lion  also  observed  that,  when 
normalized  for  organic  carbon  content,  humin  showed  the  highest 
affinity  for  hydrophobic  solutes  of  all  the  organic  matter 
studied.  Interestingly,  we  observed  the  same  phenomenon  with 
polar,  more  soluble  solutes. 
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Accurate  predictions  of  solute  -  organic  matter  partitioning 
would  thus  appear  to  require  some  knowledge  of  the  organic  matter 
involved.  However,  for  such  predictions  to  be  included  in 
realistic  models  of  the  movement  and  behavior  of  organic  solutes 
in  hydrologic  systems,  they  must  rely  on  parameters  which  can  be 
reasonably  acquired  in  a  laboratory.  We  suggest  that  the  relative 
proportion**  rf  lipids,  humic/ful  vie  acids,  and  humin  comprising 
soil  and  sedimentary  organic  matter,  along  with  the  elemental 
composition  (C,H,N,0}  of  each  fraction  which  was  first  suggested 
by  Garbarini  and  Lion,  represent  a  reasonable  first  approximation 
of  the  character  of  organic  matter. 

In  order  to  more  fully  understand  the  chemical  nature  of 
these  organic  mattet  fractions,  as  well  as  to  provide  a  more 
thorough  scientific  rationale  for  theses  suggestions,  we  have 
embarked  on  a  detailed  study  of  the  solution  properties  of  soil 
lipids,  humic/fulvic  acids  and  humin.  We  have  again  chosen 
organic-rich  peat  as  our  model  soil  organic  matter  since  this 
material,  low  in  inorganic  content,  can  be  considered  essentially 
pure  organic  matter.  This  greatly  simplifies  the  experimental 
aspects  of  the  work.  We  do  not  mean  to  imply,  however,  thac  all 
soil  and  sedimentary  organic  carbon  will  behave  as  does  peat.  At 
this  point  our  primary  interest  is  more  in  understanding  the 
differences  in  chemical  properties  of  organic  matter  fractions, 
rather  than  the  differences  between  peat  and  low-carbon  soils  and 
sediments.  The  similarities  between  the  various  organic  fractions 
of  peat  and  analogous  fractions  from  other  environments  will  be 
the  subject  of  future  investigations. 
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We  have  chosen  to  characterize  these  organic  matter 
fractions,  which  can  be  considered  mixtures  of  discrete 
geopolymers,  by  one  form  of  the  Inverse  Chromatography  Technique. 
Specifically,  we  have  measured  activity  coefficients  of  organic 
solutes  in  these  polymeric  mixtvires  by  gas  chromatography.  It  has 
been  experimentally  verified  that,  under  the  proper  conditions, 
gas  chromatography  can  provide  reliable  thermodynamic  values 
which  can  be  used  to  study  the  interactions  of  molecules  with 
different  solvents  (Conder  et  al.,  1977;  Conder  et  al.,  1986). 
Particularly  attractive  features  of  the  gas  chromatographic 
approach  are  its  essential  simplicity,  non-destructive  nature, 
the  speed  with  which  large  amounts  of  data  can  be  obtained,  and 
the  ability  to  work  at  solute  concentrations  close  to  what  is 
considered  infinite  dilution. 

In  this  work  we  have  again  used  different  fractions  of  total 
organic  matter  as  stationary  phases.  Solutes  are  injected  in  the 
gaseous  state,  and  their  retention  volumes  then  related  to  molal 
activity  coefficients  by  equation  1. 

RTWj  Pi°(Bn  -  '^i)  Po(2Bi2  “  Vi)J3^ 

In  =  Ln - - - + -  (1) 

P^VyjlO^  RT  RT 

where  is  the  activity  coefficient,  T  the  absolute  temperature,  R 
the  gas  constant,  W2  the  weight  of  stationary  phase  in  the 
column,  Pj^°  the  saturation  pressure  of  the  pure  solute,  Pq  the 
pressure  at  the  column  outlet,  Vj^  the  net  retention  volume,  and 
and  82^2  second  virial  coefficients  of  the  pure  solute  and 
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the  solute  in  the  carrier  gas,  respectively.  This  equation  was 
first  derived  by  Martire  (1974)  and  it  has  been  shown  that  if  the 
column  pressure  differential  is  less  than  two  atmospheres,  the 
carrier  gas  can  be  treated  as  an  ideal  gas,  and  the  contributions 
of  the  last  two  terms  are  insignificant.  Activity  coefficients 
less  than  one  are  interpreted  as  an  indication  of  high  affinity 
for  the  stationary  phase,  while  values  greater  than  one  indicate 
low  affinity,  or  unfavorable  mixing  of  the  solute  with  the 
stationary  phase. 

The  activity  coefficient  measurements  just  described  have 
been  supported  by  both  conventional  (CP-MAS)  and  dipolar  dephased 
(DD-CP-MAS)  Cross  Polarization-Magic  Angle  Spinning  NMR 
spectroscopy.  The  dipolar  dephasing  experiment  takes  advantage  of 
the  faster  relaxation  of  protonated  carbons.  By  allowing 
protonated  carbons  to  dephase  during  a  dipolar  coupling  period 
prior  to  data  acquisition,  the  resonances  of  C-H  and  C-H2  groups 
can  be  suppressed  so  that  the  remaining  signal  corresponds  to 
non-protonated  carbons,  as  well  as  methyl  carbons  that  couple 
weakly  because  of  their  rapid  rotational  mobility.  By  using  an 
optimum  delay  period  (40-50  usee)  the  protonated  carbon  signal 
can  be  removed  from  the  total  signal,  leaving  resonances 
primarily  from  non-protonated  carbon  atoms. 
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EXPERZMEMTAL 


Gas  Chromatography . 

All  chromatographic  measurements  were  performed  with  a 
Varian  3700  gas  chromatograph  equipped  with  a  flame  ionization 
detector  that  was  maintained  at  200°C.  For  infinite  dilution 
retention  volume  measurements  the  oven  was  keep  at  25°C  and  the 
temperature  of  the  column  was  monitored  to  +  0.1°C  with  a  copper- 
constantan  thermocouple.  The  inlet  and  outlet  pressures  were 
monitored  with  an  amplified  pressuie  transducer  (105  series, 
Onega; . 

ME  Spectroscorv 

NMR  spectra  were  obtained  with  a  Bruker  WP-200  SY 
spectrometer  equipped  with  an  IBM  solids  control  accessory 
operating  at  50  MHz  for  carbon.  The  spectrometer  was  tuned  with 
t-butyl  benzene  and  chemical  shift  scales  are  reported  as  parts 
per  million  downfield  from  tetramethylsilane.  Samples  were  spun 
at  the  magic  angle  at  approximately  4.5  KHz.  A  contact  time  of  l 
msec  was  used  throughout,  and  for  dipolar  dephasing  experiments  a 
40  usee  delay  was  inserted  prior  to  data  acquisition.  This  delay 
consisted  of  two  20  usee  intervals  separated  by  a  10  usee,  180° 
refocusing  pulse.  The  repetition  rate  for  all  NMR  experiments  was 
3  sec.  The  data  was  acquired  as  2K  sets,  zero-filled  to  8K  and 
exponentially  multiplied  prior  to  Fourier  transformation. 

Extractions 

The  sedimentary  organic  matter  used  in  this  study  was  from  a 
2  meter  peat  core  obtair.ed  from  the  Everglades  National  Park  (FL) 
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during  March  of  1986.  The  core  was  taken  from  a  Mahogany  Hammock 
adjacent  to  the  PaHayOkee  visitor  site.  The  core  was  refrigerated 
and  stored  in  a  sealed  plastic  container  under  c  nitroger, 
atmosphere  in  the  field,  then  frozen  upon  return  to  the  lab  until 
it  was  ready  for  use.  Three  horizons  of  8  cm  were  cut  from  the 
core  at  12,  35  and  151  cm  and  labeled  as  horizons  A,  B  and  C, 
respectively.  Each  sample  horizon  was  lyophilized,  mechanically 
homogenized  and  stored  in  sealed  vials  before  analysis. 

The  total  organic  matter  (TOM)  from  each  horizon  was 
continuously  extracted  with  a  2:1  mixture  of  dichloromethane  and 
methanol  until  a  clear  supernatent  was  achieved.  The  material 
extracted  in  this  operation  was  defined  as  the  lipid  fraction. 
The  residue  left  after  the  lipid  extraction  was  then  extracted 
with  0.5  N  NaOH  under  a  nitrogen  atmosphere  to  remove  humic  and 
fulvic  acids.  The  insoluble  material  was  separated  from  the 
supernatent  by  centrifugation  at  12000  RPM.  This  operation  was 
repeated  until  no  further  material  could  be  extracted  from  the 
residue  and  the  supernatent  was  clear.  The  remaining  residue  was 
then  washed  with  distilled  water  and  lyophilized.  This  is  the 
humin  fraction.  The  NaOH  solution  was  acidified  with  HCl  to  a  pK 
of  1.0  and  the  humic  acids  that  precipitated  were  separated  by 
centrifugation  at  12000  RPM,  resuspended  in  an  aqueous  solution 
of  pH  3  to  protonate  phenolic  and  carboxylic  sites,  and  then 
lyophilized.  The  isolated  fulvic  acids  were  not  used  in  this 
study. 

After  extraction  of  the  lipids,  the  residue  was  dried  and 
re-weighed.  Lipid  contents  in  horizons  A,B  end  C  were  found  to  be 
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4.3%,  3.8%  and  1.3%  by  weight,  respectively. 


Activity  Coefficient  Measurements 

Activity  coefficients  were  determinated  according  to  the 
procedures  described  by  Conder  and  Young  (1977)  .  The  carrier  gas 
flow  rate  was  adjusted  until  the  retention  volumes  of  each  probe 
showed  no  dependence  on  flow  rate.  This  procedure  was  done 
separately  for  each  of  the  fractions.  Columns  were  made  with  1/16 
inch  ID  glass  capillary  tubes.  Lipids  were  coated  onto  carefully 
silanized  100/200  mesh  silica  gel.  Humic  acids,  humin  and  TOM 
were  lyophilized  and  size  fractionated  to  100/200  mesh,  then 
packed  into  columns  without  any  additional  silica  support 
particles.  The  probes  were  injected  as  head  space  vapor  and  the 
amount  injected  was  adjusted  individ'*:ally  for  each  sample  so  that 
the  infinite  dilution  condition  was  achieved  (i.e.  no  dependence 
of  sample  size  on  retention  volume) . 


results  and  discussion 


Activity  go?ff.lgi.ent-s 

Activity  coefficients  for  the  organic  matter  fractions  from 
the  deeper,  more  altered  horizons  B  (35-42  cm)  and  C  (151-158  cm) 
are  summarized  in  Table  1.  The  solutes  used  in  these  studies  are 
arranged  from  top  to  bottom  in  Table  1  according  to  increasing 
electronic  polarizability.  Electronic  polarizability  determines 
the  contribution  of  nonspecific  interactions  to  the  total 
interaction  energy  involved  in  the  chrcmatograpn^c  process 
(Barrer,  1966) ,  and  a  correlation  between  sorption  affinity  and 
polarizability  is  an  indication  that  nonspecific  dispersion 
forces  are  responsible  for  the  observed  sorption,  thus  implying  a 
partition  process. 

Solute  activity  coefficients  are  lowest  for  the  lipid 
fractions,  as  expected,  indicating  that  the  probes  have  the 
highest  affinity  for  this  material.  Chromatographic  peaks  were 
very  symmetrical,  and  retention  volumes  showed  no  dependence  on 
flow  rate  or  concentration.  Except  for  the  relatively  low 
coefficient  for  carbon  tetrachloride  in  the  lipids  from  horizon 
C,  activity  coefficients  exhibit  a  negative  relationship  with 
polarizability  (i.e.  increasing  sorption  with  increasing 
polarizability)  Interestingly,  all  probes  except  carbon 
tetrachloride  have  higher  activity  coefficients  in  (lower 
affinities  for)  the  lipid  material  from  the  deeper  horizon, 
suggesting  that  the  qualitative  character  of  lipids  in  this 
sediment  changes  with  depth,  along  with  quantitative  changes 
(note  that  lipids  comprise  3.8%  of  the  organic  matter  in  horizon 
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B,  but  only  1.3%  in  horizon  C) . 


The  huroic  acid  fractions  show  very  little  affinity  for  the 
probes,  as  indicated  by  the  large  activity  coefficients  obtained. 
We  also  observed  that  peaks  were  less  symnetrical  than  those 
obtained  with  the  lipids,  and  the  concentrations  needed  to 
achieve  infinite  dilution  were  significantly  lower.  These  data 
and  observations  indicate  that  the  polar  and  ionic  character  of 
humic  acids  results  in  very  little  tendency  for  nonpolar  solutes 
to  partition  into  them.  A  weak  surface  adsorption  process  is 
probably  responsible  for  what  little  sorptio;,  is  observed.  There 
does  not  appear  to  be  any  significant  change  in  the  sorption 
properties  of  humic  acids  with  depth. 

Humin  showed  intermediate  affinities  for  the  probes,  with 
sorption  of  the  chlorinated  probes  being  particularly  strong. 
Peak  shapes  were  significantly  more  symmetrical  than  those  for 
the  humic  acids,  although  some  asymmetry  was  still  apparent.  As 
with  the  lipids,  increasing  affinity  correlates  with  increasing 
polarizability,  again  excepting  carbon  tetrachloride.  In  contrast 
to  the  lipids,  however,  there  is  a  very  significant  increase  in 
sorption  of  the  probes  by  the  lower  horizon.  This  data  suggest 
substantive  diagenetic  alteration  of  the  humin  fraction, 
resulting  in  a  much  less  polar,  more  lipophilic  material.  Other 
researchers  have  observed  similar  phenomena.  Shin  et  al .  (1970) 
observed  that  DDT  adsorption  was  greater  in  more  humified  soil 
organic  matter.  Garbarini  and  Lion  (1986)  demonstrated  a 
significant  correlation  between  atomic  C/0  ratios  of  various 


organic  matter  samples  and  the  affinity  of  the  organic  matter  for 
the  hydrophobic  solutes  toluene  and  trichloroethylene. 

Another  explanation  for  this  change  may  be  the  inclusion  of 
non-solvent  extractable  lipids  into  the  humin  fraction  during 
diagenesis.  In  contrast  to  the  studies  reported  in  the  previous 
manuscript  (Antworth  et  al.,  1988),  only  solvent  extractable 
lipids  were  removed  from  the  organic  matrix  here.  ^ Bound'  lipids 
must  first  be  freed  by  saponification,  ther  removed  by  solvent 
extraction.  Several  workers  have  noted  increases  in  the  amount  of 
bound  lipids  with  depth  of  burial  (Mackenzie  et  al.,  1982  and 
references  cited  therein, ,  and  it  is  presumed  that  diagenetic 
alteration  of  sedimentary  organic  matter  results  in  certain 
lipids  forming  ester,  amide,  ether,  etc.  linkages  with  the 
organic,  and  possibly  inorganic,  matrices. 

Table  1  also  includes  activity  coefficients  for  the 
original,  unfractionated  organic  matter.  As  with  the  humin  and 
lipids,  sorption  increases  with  increasing  solute  polarizability. 
Unlike  the  humin  fraction,  however,  sorption  into  the  total, 
unfractionated  sedimentary  organic  matter  decreases  with 
increasing  horizon  depth,  in  spite  of  a  significantly  more 
lipophilic  humin  component.  This  is  a  reasonable  result  if  we 
presume  that  the  non-bound,  solvent  extractable  lipid  fraction 
dominates  the  sorption  of  nonpolar  organic  solutes  by  the  intact, 
unfractionated  organic  matter,  since  this  lipid  fraction 
comprises  only  1.3%  of  the  TOM  of  horizon  C,  but  3.8%  of  the  TOM 
of  horizon  B. 
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In  order  to  further  illustrate  the  dominance  of  the  lipid 
fraction  in  the  sorption  process,  we  have  calculated  the  ratio  of 
the  solute  activity  coefficients  of  the  humic  acid,  humin  and  TOM 
fractions  to  those  in  the  lipid  fractions  for  the  same  compound 
and  depth.  The  results,  summarized  in  Table  2,  provide  an 
estimate  of  the  relative  affinities  of  these  .nonpolar  compounds 
for  the  various  fractions.  Again,  smaller  values  in  Table  2 
reflect  affinities  more  similar  to  those  for  the  lipid  fraction; 
that  is,  a  more  *  lipid-1  iJce ' ,  or  lipophilic,  material.  The  data 
of  Table  2  reinforces  the  fact  that,  with  increasing  depth  of 
burial,  the  total  organic  matter  contained  in  the  peat  sediment 
becomes  slightly  less  lipophilic,  in  spite  of  a  substantial 
increase  in  the  lipophilic  character  of  the  humin  fraction.  This 
again  suggests  that  the  free,  solvent-extractable  portion  of  soil 
and  sedimentary  lipids  dominates  the  sorption  of  nonpolar  organic 
solutes. 

Spectroscopy 

The  zero-delay  CP-MAS  NMR  spectra  of  the  total 
sedimentary  organic  matter  and  the  humic  acid  and  humin  fractions 
are  displayed  in  Figure  1.  Included  are  spectra  from  horizon  A 
(12-20  cm)  in  order  to  give  a  broader  perspective  of  the 
diagenetic  processes  within  this  core. 


The  primary  changes  observed  in  the  unfractionated  organic 
matter  are  the  relative  loss  of  intensities  in  resonances  arising 
from  alcoholic  (54  ppm) ,  carbohydrate  (74  ppm'  and  carboxyl  (175 
ppm)  carbon,  with  corresponding  increases  in  aliphatic  (30  ppm) 
and  aromatic  (130,154  ppm)  carbon.  These  changes  probably  reflect 
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removal  of  the  more  labile  carbohydrates,  leaving  a  material 
enriched  in  aliphatic  and  aromatic  structures.  Orem  and  Hatcher 
(1987)  observed  the  same  basic  trends  in  whole  peats  from  a 
sawgrass  plain  in  the  Everglades,  although  that  peat  appears  to 
contain  a  significantly  greater  carbohydrate  fraction.  Indeed, 
the  spectrum  of  Orem  and  Hatcher's  whole  peat  from  the  70-90  cm 
horizon  is  almost  identical  to  the  spectrum  we  obtained  on 
unfractionated  peat  from  the  151-158  cm  horizon. 

There  appears  to  be  little  significant  change  in  tl.e  humic 
acid  fraction  with  depth.  Slight  decreases  in  the  aliphatic  and 
carboxyl  peaks  relative  to  the  aromatic  peaks  are  noted,  however. 

Interpretation  of  the  humin  spectra  is  much  less 
straightforward.  First,  all  of  the  humin  fractions,  even  that 
from  horizon  A,  have  a  significant  aliphatic  component.  These 
spectra  thus  resemble  those  of  deeper  humin  samples  (10  cm  and 
below)  obtained  by  Orem  and  Hatcher.  The  evolution  of  our  humin, 
however,  appears  quite  different  from  that  observed  by  them, 
where  there  was  a  steady  loss  of  carbohydrate,  alcoholic  and 
carboxyl  carbon  in  humin  from  the  surface  (0-5  cm)  through  humin 
in  the  55-70  cm  horizon.  While  we  observe  this  same  trend  in 
humin  evolution  on  going  from  the  12-20  cm  horizon  to  the  35-43 
cm  horizon,  both  '  carboxyl  and  alcoholic  peaks  reappear  in  humin 
from  the  151-158  cm  horizon.  We  do  not  as  yet  fully  understand 
these  observations,  but  the  humin  in  horizon  C  may  include  some 
incorporated  humic  acids  which  are  no  longer  extractable. 
Alternatively,  there  may  have  been  some  change  in  the  sedimentary 
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environment  which  accelerated  the  production  of  extractable  numic 
acids  in  horizon  B  (or  hindered  the  humic  acid  evolution  in 
horizon  C) .  These  last  arguments  presume  that  the  hypothesis  of 
Orem  and  Hatcher  that  humic  acids  originate  as  decomposition 
products  of  humin  is  correct. 

As  a  final  word  on  the  NMR  spectra  of  these  humin  samples, 
it  should  be  noted  that  a  second  set  of  spectra  were  obtained  in 
order  to  ensure  that  we  were  not  observing  some  spectra] 
artifact.  These  spectra  were  identical  to  the  originals,  and  we 
thus  believe  that  we  are  observing  a  real  phenomenon. 

With  regard  to  the  present  study  on  the  solution  properties 
of  organic  matter  fractions,  the  most  important  and  consistent 
difference  in  humin  from  horizon  B  and  horizon  C  appears  to  be  an 
increase  in  the  relative  proportion  of  aromatic  carbon. 
Furthermore,  dipolar  dephased  NMR  spectra  (Figure  2)  show 
persistent  aromatic  resonances  after  40  usee,  indicating  that 
this  aromatic  carbon  is  significantly  nonprotonated ;  that  is, 
highly  cross-linked  to  other  carbon  moeties.  Current  efforts  are 
directed  at  a  better  understanding  of  degraded  humin  structure, 
and  why  this  structure  results  in  a  lipophilic  polymeric  material 
when  compared  to  its  less  altered  precursor. 
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SUMMARY 


The  sorption  of  nonpolar,  hydrophobic  solutes  into 
sedimentary  organic  matter,  and  the  lipid  and  humin  fractions  of 
the  organic  matter,  is  seen  to  correlate  with  electronic 
polarizability  of  the  solutes,  implying  a  non-specific 
dispersion,  or  partition,  process.  As  expected,  the  lipid 
fraction  shows  by  far  the  highest  affinity  for  these  solutes. 
Sorption  by  the  total,  unfractionated  organic  matter  seems  to  be 
dominated  by  this  extractable,  non-bound  lipid  component,  even  at 
depth  where  humin  has  acquired  a  significantly  lipophilic 
character. 

The  results  of  this  physiochemical  investigation  of 
sedimentary  organic  matter  and  its  various  components  seem  to 
support  previous  results  of  Hartley  (1960) ,  who  described  the 
fats-waxes-resin  fraction  (i.e.  lipids)  as  the  fraction  of  soil 
organic  matter  into  which  hydrophobic  solutes  partition,  as  well 
as  those  of  Shin  et  al.  (1970)  and  Garbarini  and  Lion  (1986), 
both  of  whom  observed  increasing  hydrophobic  solute  sorption  with 
increasing  maturation  of  soil  and  sedimentary  organic  matter. 
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Table  1.  Molal  Activity  Coefficients  Obtained  at  25°C. 


Lipids 

Humic 

Acids 

Humin 

TOM 

Horizon 

B 

C 

B 

C 

B 

C 

B 

C 

Depth ( cm) 

35-43  150-158 

35-43 

150-158 

35-43  150 

1-158 

35-43  150 

1-158 

^5^12 

0.036 

0.043 

>10 

>10 

5.2 

1.5 

4.5 

5.7 

C6Hi4 

0.013 

0.015 

9.0 

>10 

1.3 

0.5*. 

2.4 

2.3 

CCI4 

0.013 

0.004 

7.9 

7.0 

1.6 

0.60 

1.6 

2.0 

C2HCI3 

0.008 

0.009 

>10 

>10 

0.45 

0.20 

0.49 

0.54 

C2H4CI2 

0.003 

0.009 

8.5 

9.8 

0.36 

0.12 

0.55 

0.66 

Table  2.  Relative  Activity  Coefficients,  Normalized  to 

Corresponding  Activity  Coefficient  for  Lipid  Fraction. 


Humic 

Acid 

Humin 

TOM 

Horizon 

B 

C 

B  C 

B 

C 

Depth (cm) 

35-43 

150-158 

35-43  150-158 

35-43  150 

1-158 

Solute 

C5H12 

278 

232 

144  34 

125 

132 

^6^14 

692 

666 

100  36 

184 

153 

CCI4 

607 

1750 

123  150 

184 

500 

C2HCI3 

1250 

1111 

56  22 

61 

60 

C2K^Cl2 

2833 

3266 

120  40 

183 

220 
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FIGURE  CAPTIOKS 


Figure  1.  CP-MAS  NMR  spectra  of  total  organic  matter,  humic 

acids  and  humin  as  a  function  of  depth. 

Figure  2.  Dipolar  dephased  CP-MAS  NMR  spectra  (0-4C  usee)  of 
humin  from  35-43  and  151-158  cm  horizons. 
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ABSTRACT 


A  reinforced  concrete  arch  structure  is  modeled  under  a 
sloping  embankment  consisting  of  three  different  protective 
layers;  compacted  soil,  concrete  slab,  and  rock  cover.  The 
variables  taken  into  account  are:  bedding  condition, 
embankment  slope,  embankment  soil  type,  and  type  of  rock 
cover . 

Structural  responses  are  predicted  through  finite 
element  code,  CAKDE  (Culvert  ANalysls  and  DEslgn).  The 
results  of  CANDE  analyses  Include  nonlinear  constitutive 
models  for  characterizing  arch  structure  and  soil. 
Incremental  construction  and  free  vertical  movement  with  a 
symmetric  mesh,  and  allowed  the  comparison  of  structural 
responses  for  different  structural  systems  modeled. 

Also,  a  preliminary  lnv«*3tigatlon  to  study  the  behavior 
of  a  buried  arch  structure  using  the  centrifuge  model 
technique  Is  performed  and  described  in  the  report. 


SECTION  I 


INTRODUCTION 


A.  PURPOSE 

During  the  recent  years,  several  research  programs  were 
undertaken  by  the  United  States  Air  Force  (References  1,2,3,25) 
to  develop  the  protective  antipenetration  system  for  military 
facilities  (e.g.,  arch-,  circular-,  and  box-type  structures) 
constructed  on  the  ground  surface  (Figure  l.l).  Protective 
concepts  consisting  of  rock  rubble/boulder  and  burster  slab  over 
soils  were  Introduced,  and  full-scale  models  were  tested  to 
determine  the  levels  of  protection  against  conventional  weapons 
(Figure  1.2).  A  parametric  study  on  the  thickness  of  rock 
rubble/boulder  layer,  ratio  of  size  of  rock  rubble/boulder  to 
size  of  the  weapons,  thickness  of  burster  slab,  and  density  and 
thickness  of  sand,  layer  was  performed  and  qualitative  results 
have  been  obtained  from  the  field  model  tests.  Based  on 
Information  observed,  the  system  has  proven  effective  to  resist 
the  penetration  of  conventional  weapons  and  It  is  now  being 
considered  for  application  over  military  facilities  as  a 
hardening  protective  structure.  However,  the  designers  are 
facing  one  critical  question  before  the  application  of  the 
system  (Reference  25);  Will  these  military  structures  be  safe 
under  static  loadings  due  to  the  weight  of  embankment  soil,  rock 
rubble/boulder  layers  and/or  burster  slab?  The  purpose  of  the 
research,  therefore,  Is  to  study  the  behavior  of  the  buried  arch 
structuies  subjected  to  static  loadings  by  using  centrifuge  and 
numerical  modeling  techniques. 
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b)  Rock/Rubble  Protective  Layer 
Figure  1.2 
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B.  OBJECTIVE 

The  primary  objective  of  this  study  is  to  Investigate  the 
soil-structure  interaction  problems  of  a  buried  concrete  arch 
structure  subjected  to  static  landings.  More  specifically,  the 
objectives  of  this  study  are  to: 

1.  Develop  an  experimental  technique  to  simulate  the 
behavior  of  a  prototype  arch  structure  in  a  centrifuge. 

2.  Conduct  a  preliminary  series  of  centrifuge  model 
experiments  on  a  small  capacity  centrifuge. 

3.  Predict  the  behavior  of  the  arch  using  finite  element 
analysis . 

4.  Gain  experience  necessary  for  model  studies  In  a  large 
capacity  centrifuge. 

C.  SCOPE 

With  these  objectives  in  mind,  the  research  is  reported  in 
the  following  sections.  They  contain  a  description  of  the 
centrifuge  facility  used  for  model  study;  design  and  fabrication 
of  an  accompanying  testing  apparatus;  preliminary  model  test 
results,  and  recommendation  for  further  centrifuge  model  study; 
and  investigation  of  the  influence  of  soil  types  for  in-situ 
bedding  and  backfill,  and  slope  of  the  embankment  using  finite 
element  analysis.  Finally,  some  potential  topics  for  future 


research  in  this  area  are  mentioned. 


SECTION  II 


CENTRIFUGE  MODEL  STUDY 

A.  FUNDAMENTALS  OF  CENTRIFUGE  MODELING 

One  of  the  most  ideal  approaches  for  obtaining  information 
on  the  behavior  of  prototype  structures  is  full-scale  model 
testing.  A  full  scale  model  with  the  necessary  instrumentation 
(i.e.,  soil  stress  meters,  pore  water  pressure  transducers, 
settlement  gauges  and  strain  gauges,  etc.)  could  give  the  best 
results  for  estimating  prototype  behavior. 

Unfortunately,  full-scale  model  testing  has  serious  major 
drawbacks:  mainly,  cost  and  time  of  construction  and  operation. 

Because  of  these  reasons,  small  scale  model  testing  is  becoming 
a  favorite  testing  method  in  geotechnical  engineering. 

However,  use  of  small  scale  model  tests  in  the  laboratory 
is  severely  limited  when  the  gravity  body  force  of  the  structure 
Itself  is  the  principal  load  on  the  system,  such  as  in  dams  and 
embankments.  This  limitation  is  due  to  two  major  factors.  One 
is  that  soil  characteristics  are  nonlinear  and  ove''burden 
dependent,  and  the  other  Involves  stress  magnitudes.  The 
stresses  in  a  small  scale  model  due  to  its  own  weight  are  much 
smaller  in  magnitude  than  those  in  the  corresponding  prototype 
system.  To  eliminate  these  deficiencies,  the  centrifuge 
modeling  technique  has  been  used  since  1931.  The  increase  in 
unit  weight  is  achieved  by  placing  a  model  in  a  centrifuge  and 
spinning  it  to  produce  an  acceleration  field  that  is  equivalent 
to  an  artificial  gravitational  field.  In  this  way  the  state  of 
stress  at  every  point  in  the  model  under  an  artificial 
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gravitational  field  is  equal  to  that  at  the  homologous  point  In 
the  prototype.  This  technique  has  been  applied  to  a  variety  of 
geotechnical  problems  (i.e.,  slope  stability,  reinforced  earth, 
foundation  design,  offshore  gravity  platforms,  rockflll  dams, 
tunnels,  and  burled  circular  pipes)  reported  elsewhere 
(References  5,10,13,14,23,24,28,29  and  30).  Table  2.1 
summarizes  the  scaling  rules  that  have  been  developed  for 
relating  the  behavior  of  the  centrifuge  model  and  an  equivalent 
full  scale  geotechnical  structure. 

B.  THE  AFESC  CENTRIFUG'i  FACILITY 

1.  The  Centrifuge 

The  centrifuge  Installed  at  Tyndall  Air  Force  Base  (AFB) 
was  originally  built  for  testing  of  electronic  and  mechanical 
packages  under  G-loadings  that  simulate  those  experienced  in 
actual  flight.  It  was  initially  installed  at  Kirtland  AFB,  New 
Mexico.  The  machine  was  modified  for  geotechnical  studies  in 
1981  and  moved  to  the  Air  Force  Engineering  and  Services  Center 
(AFESC)  at  Tyndall  AFB  in  1986.  The  centrifuge  was  assembled 
and  housed  in  a  7-ft.  high,  16-ft.  diameter,  9-10.  thick 
reinforced  concrete  retaining  structure  (Figure  2.1). 

The  Center's  centrifuge  is  a  medium  size  rotary  accelerator 
designed  by  Genlsco,  Inc.,  of  California  (Model  E-185).  The 
machine  is  capable  of  applying  ccntrolled  centrifugal 
accelerations  up  to  100  g's  and  a  limit  of  30,000  g-lb.  at  a 
nominal  radius  of  6  ft.  The  payload  of  the  test  package  can  be 
as  heavy  as  500  lb.  The  machine  includes  a  variable-speed 
hydraulic  drive  system,  timing  belt  drive,  RPM  pick-up,  rotating 
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Table  2.1  Scaling  Relations  (Reference  29) 


Quantity 

Full  Scale 
(Prototype) 

Centrifugal 
Model  at  n  g’s 

Linear  Dimension 

1 

1  /n 

Area 

1 

Volume 

1 

l/n^ 

Time 

In  Dynamic  Terms 

1 

1/n 

In  Diffusion  Cases 

1 

l/n^ 

In  Viscous  Flow  Cases 

1 

1 

Velocity  (Dlstance/Tlme) 

1 

1 

Acceleration  (Distance/Time) 

1 

n 

Mass 

1 

l/n^ 

Force 

1 

l/n^ 

Energy 

1 

1/n^ 

Stress  (Force/Area) 

1 

1 

Strain  (Displacement/Unit  Length  1 

1 

Density 

1 

1 

Energy  Density 

1 

1 

Frequency 

In  Dynamic  Terms 

1 

n 
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arm  asseably,  ternlnal  box  for  control  and  test  connections. 

The  hydraulic  drive  system  consists  of  an  electric  drive 
motor  coupled  to  a  variable-displacement  pump,  a 
constant-displacement  fluid  motor  couple  to  the  rotor  drive 
shaft,  and  two  pressure-relief  valves.  A  A-way 
solenoid-operated  spring-loaded  control  valve  provides 
rotational  stability,  maximum  hunting,  and  constant  torque 
characteristics.  The  hydraulic  power  unit  is  enclosed  in  the 
console,  and  the  fluid  motor  is  Installed  at  the  bottom  of  the 
rotor  assembly  (Figure  2.2).  Maximum  pressure  of  the  hydraulic 
system  has  been  set  at  approximately  2,200  psl. 

The  boom  consists  of  two  symmetrical  cantilever  arms,  and 
an  adjustable  30-ln.  square  cradle-type  mounting  platform 
carried  on  aircraft-type  roller  bearings.  The  platform  is  held 
by  two  arms  attached  to  the  spokes  by  two  pivots.  It  can  be 
locked  into  horizontal,  45,  90,  135  or  180  degree  positions. 

This  permits  the  soil  surface  to  remain  perpendicular  to  the 
vector  sum  of  the  centrifuge  acceleration  and  the  acceleration 
due  to  gravity. 

One  of  the  unique  features  of  the  machine  is  the  automatic 

dynamic  balancing  system.  Dynamic  balance  is  accomplished 

automatically  by  vertical  displacement  of  the  two  arms  of  the 

boom  until  the  center  of  gravity  of  the  opposing  mass  lies  in 

the  same  horizontal  plane.  It  requires  approximately  20  seconds 

for  full  compensation.  If  compensation  cannot  be  accomplished 

within  the  12-in.  travel  limitations,  the  balance  motor  is 

automatically  shut  down  by  limit  switches.  Corrective  measures 

are  indicated  by  the  relative  attitudes  of  the  mounting 
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Figure  2.2  Rotary  Accelerator  (Model  E-185)  -  Outline  Drawing 
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2.  Video  Recording  System 

Video  camera  and  recording  devices  are  also  available  as 
part  of  the  monitoring  system.  This  video  camera  Is  mounted 
near  the  hub  of  the  centrifuge.  This  system  provides  a 
continuous  and  Instantaneous  monitor  of  the  test  while  In 
progress,  and  a  permanent,  replayable  record  of  model  tests. 

3.  Data  Acquisition  System 

Digital  oscilloscopes  (Nlcolet  Model  4094)  have  been  used 
for  collecting  sensor  signals  from  model  structures  for 
centrifuge  model  study.  The  Model  4094  digital  oscilloscope 
with  multi-channel  plugs  Is  capable  of  performing  acquisition 
with  12  bit  resolution  at  a  rate  of  2  million  sample/sec  (16 
bits  samples).  The  system  allows  zero  time  to  be  set  anywhere 
In  the  displaying  area  and  enables  two-channel  measurements  with 
each  channel  set  at  an  Independent  pre-  or  post-trigger  delay. 

It  displays  these  two  functions  as  well  as  a  combination  of 
them.  Also,  the  system,  a  combination  of  them.  A  HP  plotter 
Interfaced  with  the  oscilloscope  has  been  used  to  create  hard 
copy  records. 

C.  FABRICATION  OF  MODEL  BOX 

A  requirement  of  centrifuge  modeling  Is  the  fabrication  of 

a  model  box  chat  contains  the  small-scale  model  of  the  structure 

together  with  the  soils  sample.  The  dimensions  of  the  model  box 

were  made  small  enough  to  fit  the  centrifuge  swing  bucket 

platform,  and  large  enough  to  house  the  structure  model  and  a 

sufficient  amount  of  soil  sample.  The  model  box  has  a  platform 
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base  of  2A  in.  x  24  in.  and  a  total  height  of  17.875  in.  It  was 
totally  made  of  2024-T-35  aluminum  plates  welded  at  the  sides 
and  the  bottom  at  almost  100%  penetration.  The  model  box  has  a 
top  lid  as  well  a  front  lid  for  easy  access  to  the  model  and 
cleaning  purposes.  Both  lids  are  secured  by  steel  screws,  the 
front  lid  is  designed  to  withstand  the  internal  lateral  pressure 
caused  by  the  soil  at  increased  gravity  conditions  (100  g's), 
while  the  top  lid  is  designed  to  withstand  the  pressure 
generated  by  the  surcharged  system  (50  psi). 

The  base  of  the  model  box  extends  2  in.  from  every  side  and 
has  a  hole  drilled  into  each  corner  so  as  to  enable  the 
researchers  to  bolt  down  the  box  onto  the  centrifuge  platform  to 
prevent  it  from  slipping.  Although  theory  states  that  the  box 
would  not  slip  off  the  swing  bucket  while  the  centrifuge  is 
running,  this  is  merely  a  safety  precaution  taken  by  the 
researchers.  For  easy  and  safe  handling  of  the  model  box,  the 
top  edges  of  the  two  opposing  sides  are  extended  and  holed-out 
to  accomodate  crane  hooks.  A  detailed  diagram  of  the  model  box 
is  Illustrated  in  Figure  2.3 

D.  ARCH  STRUCTJRE  MODEL  AND  STRAIN  GAUGES 

The  model  semi-circular  concrete  arch  structure  was  made  of 
a  microconcrete  mix  composed  of  gypsum  molding  plaster  and 
silica  sand,  with  gypsum: aggregate twater  ratio  of  1:2:1  by 
weight.  The  structure  has  an  inner  diameter  of  9  in.  and  a 
shell  thickness  of  0.3  in.  The  total  length  was  16  in.  The 
structure  model  was  complete  with  miniature  reinforcement  set  on 
both  circumferential  and  longitudinal  directions.  Reinforcement 
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Figure  2.3  The  Centrifuge  Model  Box 
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bar  dlaneter  was  22'-gauge.  The  model  arch  structure  also 
Included  a  reinforced  floor  slab  made  of  the  same  material  (see 
Figure  2.4). 

Five  pairs  of  strain  gauges  (gauge  type  CEA-13-125UW-120) 
manufactured  by  Measurement  Group,  Inc.  were  used  to  measure 
tcraln  along  the  circumference  of  the  structure  at  midspan. 
Gauges  were  mounted  at  5,  45,  90,  135,  175  degrees  and  hooked  up 
In  a  full-bridge  fashion.  The  arch  structure  with  attached 
strain  gauges  placed  in  the  model  container  is  shown  In  Figure 
2.5. 


E.  MODEL  PREPARATION 

1.  Miniaturized  Concrete  Arch  Structure 
The  structures  used  In  the  tests  are  l/60th  size  of  the 
prototype  structure.  It  Is  made  of  a  microconcrete  mix  composed 
of  gypsum  plaster  as  cement  and  silica  sand  as  aggregate. 

Scaling  down  the  prototype  to  l/60th  size  resulted  with  a  model 
that  requires  a  membrane  thickness  of  0.3  Inch.  Due  to  this, 
the  consistency  of  the  microconcrete  mix  had  to  be  very  slurry 
for  It  to  go  down  the  mold.  The  final  mix  had  a 
gypsum : aggregate  :water  ratio  of  1:2:1. 

The  mold  is  made  up  of  two  aluminum  tubings,  one  with  an 
outer  diameter  equal  to  the  desired  innner  diameter  of  the 
model,  the  other  with  an  inner  diameter  equal  to  the  desired 
outer  diameter  of  the  model.  The  frame  work  is  made  of 
plexiglass  which  also  serves  to  hold  the  reinforcement  cage  in 
place  while  letting  the  mix  to  cure. 

The  miniaturized  reinforcement  used  was  22-gauge  wires  in 
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Figure  2.5  The  Instrumented  Arch  Structure  in  the  Model  Box 


22-19 


both  the  longitudinal  and  circumferential  directions.  Prior  to 
the  installation  of  the  wires.  It  was  necessary  to  apply  a  small 
amount  of  pretension  to  remove  all  bends  and  colls.  Also,  It 
was  necessary  to  roughen  up  the  surface  of  the  wires  with  very 
coarse  sand  paper  to  assure  better  bonding  to  the  microconcrete. 

The  models  were  left  to  cure  for  two  days  and  then  were 
coated  with  concrete  shellac  to  preserve  the  water  content. 
Allowing  the  models^to  dry  for  more  than  the  set  amount  of  time 
produced  a  very  chalky  and  brittle  surface.  After  appllcaton  of 
shellac,  the  structure  was  left  to  cure  for  at  least  another  two 
days  before  testing. 

2.  Soil  Fluvlator 

To  attain  reliable  data  from  the  tests,  most  of  the 
variables  must  be  held  constant  for  every  test.  This  Includes 
the  soli  sample  properties.  Although  grain  site  and  sample 
height  of  the  soil  can  easily  be  produced.  It  Is  the 
reproduction  of  the  soil  sample  density  that  Is  Important.  For 
this  purpose,  the  researchers  use  the  sand  raining  method  or, 
pluvlatlon.  This  method  Is  commonly  used  due  to  its  simplicity 
and  for  the  simple  reason  that  it  Is  the  best  method  of 
reproducing  soli  samples.  The  purpose  Involves  forming  a  sand 
mass  to  a  homogeneous  condition  with  a  specified  density,  using 
a  process  that  can  be  repeated  as  many  times  without  the 
variation  of  the  sand  properties  other  than  as  chosen  by  the 
researchers . 

A  sand  pluvlator  was  designed  and  built  at  AFESC/RDCO 
facility  at  Tyndall  Air  Force  Base  (Figure  2.6).  It  has  a  total 
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Figure  2.6  The  Soil  Pluviator 


height  of  approxlnately  90  Inches  and  covers  a  tnlnlfflum  area  of 
12  square  feet*  The  pluvlator  is  capable  of  producing  sand 
samples  of  extreme  densities.  By  selecting  the  proper 
perforation  board,  the  researcher  can  control  the  desired 
discharge  rate  resulting  to  a  specific  sample  density. 

The  soil  strata  used  as  a  model  consisted  of  three 
pluvlation  steps.  The  first  step  was  pluvlatlng  a  shallow  sand 
layer  on  the  bottom  of  the  model  box  to  act  as  the  bedding 
layer.  At  the  completion  of  this  atep,  the  model  box  was  pulled 
out  from  underneath  the  pluvlator  and  the  model  was  then  put  in 
place.  The  second  step  was  to  bury  the  structure  to  the  desired 
embankment  depth  using  the  same  process.  Again,  the  model  box 
is  pulled  out  from  underneath  the  pluvlator  to  set  the  burster 
slab  model  In  place.  The  final  layer  of  sand  is  then  pluvlated 
and  the  model  box  is  ready  to  be  mounted  on  the  centrifuge  swing 
bucket . 


3.  Loading  Model  Box 

The  first  step  taken  before  loading  the  model  box  onto  the 
centrifuge  swing  bucket  was  to  weigh  the  model  box  complete  with 
the  lid  with  its  contents  to  determine  the  amount  of 
counterweight  to  be  used.  Since  the  model  box  was  set  on  tne 
center  of  the  platform,  the  counterweight  needed  was  equal  to 
box  weight  Itself  (+  5  lbs).  As  long  as  the  red  indicator  light 
remains  off  after  setting  the  model  box  and  the  counterweight, 
the  system  is  considered  balanced.  Note  that  the  weight  on 
either  side  should  not  exceed  500  lbs.  The  model  box  was  then 

secured  onto  the  centrifuge  swing  bucket  by  bolting  it  down  on 
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the  platform  (Figure  2.7). 

As  soon  as  the  model  box  Is  secured  onto  the  swing  bucket, 
all  Instrumentation  is  wired  up  to  the  slip  rings  taking  note  of 
all  channel  designations.  At  the  centrifuge  control  panel,  all 
wires  were  hooked  up  to  designated  ports  on  the  channel  switch 
and  onto  the  strain  indicator. 
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SECTION  III 


FINITE  ELEMENT  ANALYSIS 

A.  DESCRIPTION  OF  THE  COMPUTER  CODE 

The  computer  code,  CANDE  (Culvert  ANalysls  and  DEalgn), 
uaed  In  this  study  was  developed  by  Katona  et  al.  (References 
11  and  12)*  The  basic  assumptions  of  the  program  are:  plane 
strain  geometry  and  loading,  small  displacement  theory,  and 
quasistatic  response.  The  following  description  summarizes 
salient  features  of  CANDE  used  in  the  analysis. 

1.  Modeling  of  Structure 

The  box  culvert  was  modeled  as  a  series  of  straight  line 
segments  with  a  one^'dimensional  beamocolumn  bending  element 
comprising  each  segment*  This  bending  element  employs  a 
plane-strain  formulation  which  neglects  shearing  deformations. 
Two  nodes  with  three  degrees  of  freedome  per  node  (horizontal 
and  vertical  displacement  and  a  rotation)  define  each  element. 

At  each  node,  cubic  and  linear  displacement  approximations  were 
used- in  calculations  of  displacements  in  the  transverse  and 
axial  directions,  respectively.  Figure  3.1  shows  the  coordinate 
system  of  the  beam  element. 

2.  Modeling  of  Soil 

Nonconforming  quadrilateral  and  triangular  elements  were 
used  to  represent  the  soil  (see  Figure  3.2).  The  quadrilateral. 
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defined  by  four  nodes  with  two  degrees  of  freedom  (horizontal 
and  vertical  displacements)  per  node  comprise  two  triangles  with 
complete,  quadratic,  Interpolation  functions  specified  within 
each  triangle.  Application  of  appropriate  constraints  and 
static  condensation  procedures  result  in  a  four-node 
quadrilateral  element  (Reference  11). 

There  are  four  material  characterizations  available  for  the 
soil  properties:  linear  elastic;  Incremental  elastic 
(overburden  dependent);  extended-Hardln;  and  Duncan’s  hyperbolic 
stress-strain  relationship,  which  employs  tangent  Young's 
modulus  and  tangent  bulk  modulus  formulations  (R'cif erences  7,  12, 
33). 


3.  Modeling  of  Slip 

A  constrained  finite  element  formulation,  which  Is  based  on 
a  generalized  principle  of  virtual  work,  was  used  to  model  the 
relative  movement  of  the  soil  with  respect  to  the  pipe  at  the 
soil-structure  Interface,  and  relative  movement  of  the  soil  with 
respect  to  the  soil-soil  Interface  (l.e.,  trench).  An  Interface 
element  Is  defined  by  a  set  of  paired  nodes  joining  two  elements 
as  shown  In  Figure  3.3.  The  paired  nodes  which  have  two  de.gree 
of  freedom  per  node  (vertical  and  horizontal  displacement), 
initially  take  the  same  position  In  space  before  any  loading, 
but  are  assigned  to  separate  elements.  Therefore,  each  node 
responc‘s  individually  under  any  applied  loading.  In  addition,  a 
third  node,  which  is  assigned  between  the  paired  nodes,  provides 
normal  and  shear  forces  existing  between  the  oaired  nodes.  It 

should  be  observed  that  these  Interface  forces  arise  only  from 
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Figure  3.3  Nodes  for  Constraint  Interface  Element  in  CANDE 
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the  adjoining  element  interaction  and  applied  forces  at  the 
paired  nodes. 

A.  Modeling  of  Incremental  Construction 
CANDE  uses  Incremental  Iterative  solution  techniques  to 
represent  the  placing  of  the  culvert  and  embankment  In  a  series 
of  soil  layers.  The  basis  of  the  techniques  Is  superposition  of 
solutions  from  successive  soil  layers  and  a  newly  added  layer. 

At  each  increment,  Iterative  calculations  determine  soil  and 
structure  moduli  until  equilibrium  Is  approximated  within  an 
allowable  error. 

B.  FINITE  ELEMENT  PARAMETERS 

1.  Finite  Element  Mesh  and  Boundary  Conditions 
The  finite  element  grid  with  boundard  conditions  used  In 
this  study  Is  shown  In  Figure.  3. A.  Since  the  model  and  Its 
loading  is  symmetric,  only  half  of  the  mode  was  analyzed. 
Boundary  conditions  on  the  culvert  used  In  the  analysis  are  show 
In  the  figure:  a  fixed  movement  conditions  In  the  horizontal 
direction  and  free  movement  condition  In  the  vertical  direction. 
The  arch  structure  was  represented  by  12  beam-column  elements. 


2.  Soil-Culvert  Interface  and  Incremental  Solution 
Procedure 

A  fixed  condition  was  used  at  the  soil-culvert  interface. 
Although  it  is  possible  to  represent  the  slip  between  the 
continuum  (soil)  and  beam  elements,  the  fixed  condition  was  used 
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Figure  3.4a  Mesh  No. 
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in  this  preliminary  study  due  to  large  amounts  of  CPU  (Central 
Processing  Unit)  time,  and  substantial  computational  costs  when 
slip  was  allowed. 

The  incremental  solution  procedure  for  embankment  and 
trench  simulated  the  actual  Installation  process  of  placing  soil 
layers  in  a  series  of  lifts.  Figure  3.5  shows  the  construction 
increment  numbers  of  element  groups  entering  sequentially  into  a 
system.  The  first  construction  increment  Included  placing  all 
bedding  pad,  in  situ  soil,  and  the  arch  structural  elements. 
Subsequent  Increments,  numbers  2  through  9,  were  gravity  loaded 
elements  of  fill  soil. 

3.  Material  Properties 

The  soil  model  employed  in  the  study  is  a  characterization 
proposed  by  Duncan  et  al.  (References  7  and  12)  which  has  had  a 
substantial  history  of  development  and  application  over  the  last 
decade.  This  soil  model  characterizes  soil  behavior  with  a 
variable  tangent  Young's  modulus  and  tangent  bulk  modulus. 

Eight  parameters  are  needed  to  define  a  particular  soil  in 
loading:  R,  n,  R^,  c,  0 ,  andA0  for  the  tangent  Young's  modulus; 

and  and  m  .for  the  tangent  bulk  modulus.  Table  3.1  shows  the 
material  properties  used  for  soil,  concrete  arch,  rock,  and 
burster  slab  in  the  analysis. 
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Figure  3.5  Incremental  Construction 
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Table  3.1  Duncan's  Material  Characterization  Used  for  Soil  Model 


CONCRETE 


Coirpressive  strength,  -Pc 

5,000  loi 

Young's  modulus  in  elastic  range,  E 

4,070,000  pel 

Poisson's  ratio,  t) 

0.25 

Unit  weight,  TT 

150  pcf 

Strain  at  which  tensile  cracking  occurs 

0. 0001 

Strain  at  elastic  limit  in  compression 

0.0013 

Strain  at  initial  compressive  strength. 

0. 0040 

REINFORCING  STEEL 

Yeild  stress,  -fij 

60,000  psi 

Young's  modulus,  C 

29,000,000  pei 

Poisson's  ratio,  V 

0.3 

Spacing  of  reinforcement 

1.5  in. 

Area  of  inner  cage  reinforcement 

0.31  sq.  in./in. 

Area  of  outer  cage  reinforcement 

0.31  sq.  in./in. 

Thickness  of  concrete  cover  to  center  of 
inner  cage 

1.25  in. 

Thickness  of  concrete  cover  to  center  of 
outer  cage 

1.25  in. 

Table  3.2  Concrete  and  Steel  Properties  Used  for  Arch  Structure 
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Type  Young's  modulus  Poisson's  ratio  Unit  weight 

E  V  le 

Rock  1,500,000  pei  0.2  86.94  pcf 

Burster  slab  4,070,000  psi  0.25  145.0  prf 


Table  3.3  Material  Properties  used  for  Bock  and  Burster  Slab 
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SECTION  IV 


RESULTS  AND  DISCUSSIONS 

A.  CENTRIFUGE  MODEL  STUDY 

Since  the  authors  are  -the  first  time-users  of  the 
centrifuge  facility  at  Tyndall  Air  Force  Base  with  a 
small-scale  structure  and  instrumer tat  Ion ,  a  series  of 
centrifuge  model  tests  was  first  performed  as  a  pilot  study  to 
check  the  proper  operations  of  the  centrifuge  machine  and 
Instruments.  During  this  preliminary  test,  the  Instruments 
faced  two  major  difficult  problems:  (1)  a  limited  number  of 
slip  rings  available  for  data  transmission,  and  (2)  brush  noise 
caused  by  dynamic  vibration  occurred  at  the  centrifuge  slip 
rings.  The  voltage  fluctuations  generated  by  the  noise  at  the 
slip  rings  were  found  to  be  much  greater  than  the  Instrument 
signals  themselves.  The  major  centrifuge  model  tests, 
therefore,  were  not  able  to  be  performed  during  the  stay  at 
Tyndall  Air  Force  Base  due  to  the  limitation  and  unavailability 
of  testing  equipments.  Rather,  this  section  discusses  the 
various  computer-based  data  acquisition  systems  required  to 
overcome  these  problems  for  a  future  centrlifuge  model  study. 
The  description  given  below  is  a  brief  one;  readers  Interested 
In  more  detailed  information  should  refer  elsewhere  (References 
8,9,17-22,26,32). 

Transmission  of  Signals; 

Direct  transmission  of  signals  is  used  the  most  in  model 
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studies  mainly  because  of  Its  simplicity  and  relatively  low 
cost.  Also  there  Is  no  distortion  of  signals.  When  used  In 
centrifuge  modeling,  direct  transmission  consists  of  the 
physical  linkage  of  system  components  by  wire  cables. 
Brush-and-sllp  ring  riggings  are  a  commonly  used  method  .of 
providing  a  signal  path  from  the  centrifuge  shaft  to  the 
outside.  Brush  noise  Is  caused  by  dynamic  vibration  In  the  area 
In  which  the  brush  Is  in  contact  with  the  slip  ring.  This 
vibration  In  contact  area  gives  a  variation  In  circuit 
resistance.  An  Increase  In  signal  level  before  reaching  slip 
rings,  through  amplication,  can  give  a  desirable  Increase  In  the 
signal-to-nolse  ratio.  Another  beneficial  technique  is  to 
decrease  the  total  brush  resistance  by  providing  multiple 
parallel  current  paths;  noise  would  be  reduced  because 
variations  In  brush  contact  for  Individual  brushes  would  not  be 
as  signlflcnt  when  other  brushes  are  in  the  path.  If  a  lirge 
number  of  slip  rings  dedicated  to  signal  transmission  are 
available,  the  direct  transmission  system  should  be  considered. 
It  has  the  advantage  of  simplicity  and  uninterrupted 
transmission  of  signals. 

However,  for  most  of  the  smaller  size  centrifuges,  there 
are  only  a  limited  number  of  slip  rings  for  data  transmission.  A 
multiplexer-demultiplexer  unit  (Indirect  method)  is  therefore 
used  to  facilitate  the  transmission  of  signals  from  a  relatively 
large  number  of  sensors.  The  two  most  commonly  used 
multiplexer-demultiplexer  units  are  the  time  division  unit  and 
the  frequency  modulation  (FM)  frequency  division  unit.  The  time 
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divisioTi  unit  works  as  a  synchronized  channel  selector  as  shown 
in  Figure  4.1  where  the  frequency  division  unit  works  as  a 
dedicated  frequency  carrier  along  with  FM  detectors  as  shown  in 
Figure  4.2.  Of  the  two  indirect  methods  of  signal  transmission, 
the  time  division  technique  is  relatively  simple  and  less  costly 
to  Install;  however,  additional  control  signals  are  needed  to 
synchronize  the  functions  of  the  multiplexer  and  demultiplexer. 
This  additional  control  and  operation  time  impedes  the  sampling 
rate  of  othe  system.  The  frequency  division  technique,  on  the 
other  hand,  is  costly  but  capable  of  transmitting  signals  at 
much  higher  rates.  Since  the  FM  system  is  sensitive  only  to 
frequency  variations  and  not  amplitudes,  it  is  immune  to  most 
noises  induced  along  the  transmission  line. 

It  should  also  be  mentioned  that  wireless  sensors  can  be 
adopted  to  transmit  radio  frequency  signals  thus  bypassing  the 
use  of  connecting  wires  between  sensors  and  the  signal 
conditioning  unit.  However ,  these  sensors  are  very  expensive 
and  usually  too  big  in  size  and  too  heavy,  not  readily  suitable 
for  geotechnical  use.  The  author  does  not  foresee  the  need  to 
use  wireless  sensors  in  centrifuge  modeling  in  the  future. 


B.  FINITE  ELEMENT  ANALYSIS 

1.  The  Influence  of  Installation  Type 
The  first  part  of  this  investigation  deals  with  the 
influence  of  the  type  of  installation  on  the  relative 
displacement  of  the  crown  of  the  structure  (refer  to  Figures  4.3 
thru  4.6).  The  comparisons  have  been  made  between  a  system  with 
a  3-layer  rock  cover  versus  one  with  a  5-layer  rock  cover.  All 
other  variables  were  considered  in  the  comparison  to  determine 
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CL85  BEDDING,  38  DEG  SLOPE 


Figure  M.3  Relative  Crov/n  Deflection  vs 
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Figure  Relative  Crown  Deflection 
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Figure  U.5  Relative  Crown  Deflection  vs.  Fill  Height  for  Case  I2A 
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one  with  a  5~layer  rock  cover.  All  other  variables  were 
considered  In  the  comparison  to  determine  the  degree  of 
Influence  of  Installaton  type  on  the  displacement  of  the  crown. 

For  Case  IIA  where  the  structure  Is  modeled  on  CL85  (clay) 
soil  type  bedding  and  30-degree  slope  embankment,  the  finite 
element  analysis  predicted  a  steel  safety  factor  equal  to  1.00 
at  the  end  of  eight  construction  Increments  when  a  3-layer  rock 
cover  Is  Installed,  meaning  that  the  stell  reinforcement  of  the 
structure  had  failed  under  these  conditions.  A  5-layer  rock 
cover  proved  to  be  the  best  alternative  In  this  case. 

As  for  Case  IIB  where  the  structure  Is  modeled  on  SM85 
(sand)  soil  type  bedding  and  30-degree  slope  embankment,  the 
Influence  of  the  type  of  rock  layer  embankment  can  be  seen  at 
fill  heights  greater  than  112  ft.  It  can  clearly  be  seen  that 
the  system  under  a  5-layer  rock  embankment  caused  less  crown 
deflection  than  the  3-layer  type.  Minimum  relative  crown 
deflection  at  fill  height  of  182  ft.  Is  4.76  Inches  and  occurs 
when  the  embankment  Is  of  SC90  (sand-clay)  type  soil. 

Case  I2A  deals  with  the  structure  modeled  on  CL85  soil 
bedding  and  45-degree  slope  embankment.  With  the  case  where  the 
5-layer  rock  embankment  system  was  Installed,  the  finite  element 
analysis  predicted  a  stell  safety  factor  and  concrete 
compression  safety  factor  both  equal  to  1.00  at  the  end  of  four 
construction  Increments,  meaning  that  the  structure  had  failed 
under  these  conditions.  In  this  case,  therefore,  a  3-layer  rock 
cover  would  be  the  best  and  only  solution. 

The  last  installation  case  examined  is  Case  I2B  where  the 
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structure  is  modeled  on  SM85  soil  type  bedding  and  a  45-degree 
slope  embankment*  In  this  case  the  Influence  of  the  type  of 
rock-layer  Installed  Is  much  more  clearly  defined.  Case  I2B 
shows  that  the  system  with  a  3-layer  rock  cover  causes  less 
deflection  of  the  crown  of  the  structure.  Differences  In 
deflection  of  the  crown  Is  close  to  the  1.25-lnch  range.  Even 
at  a  higher  fill  height,  the  system  with  a  3-layer  rock  cover 
causes  less  deflection  than  the  system  with  the  5-layer  rock 
cover.  Minimum  relative  crown  deflection  at  a  fill  height  of 
182  ft  is  5.31  inches  and  occurs  when  the  embankment  soli  Is  of 
SC90  type  and  the  rock  cover  Is  of  the  3-layer  type. 

The  second  part  of  the  investigation  deals  with  the 
influence  of  installation  type  on  the  bending  moment  developed 
at  the  crown  of  the  structure  depending  upon  the  type  of  rock 
cover  installed  on  the  system  (refer  to  Figures  4.7  thru  4.10). 
Similar  to  the  discussion  above,  the  comparison  has  been  made 
between  a  system  with  a  3-layer  rock  cover  versus  one  with  a 
5-layer  rock  cover.  In  Case  IIA,  the  numerical  analysis  of  the 
system  predicted  the  failure  of  the  structure  under  a  3-layer 
rock  cover  at  the  end  of  eight  construction  increments.  And  as 
for  Case  I2A,  the  same  situation  arised  with  a  structure  under  a 
5-layer  rock  cover  but  after  the  application  of  only  four  load 
increments.  Therefore,  the  comparisons  were  limited  to  two 
cases . 

Case  IIB  showed  that  at  a  fill  height  greater  than  182 
feet,  the  system  with  the  5-layer  rock  cover  produced  less 
bending  moment  at  the  crown  compared  to  the  one  with  the  3-layer 
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Figure  4.9  Crown  Moment  vs.  Fill  Height  for  Case  I?A 
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rock  cover.  The  difference  between  the  moments  occurring  at  the 
crown  due  to  the  two  types  of  rock  cover  averaged  7988.667 
In-lb/in.  Minimum  moment  developed  at  the  crown  of  the 
structure  has  a  magnitude  of  78640  in-lb/in  and  occurs  when  the 
embankment  is  of  SM90  type  soil. 

Case  I2B  showed  that  the  structure  modeled  with  a  3“layer 
rock  cover  produced  much  less  moment  at  the  crown  *  compared  to 
that  with  a  5-layer  rock  cover.  The  average  difference  of  the 
moment  developed  at  the  crown  is  45594.0  in-lb/ln  and  mlniniutn 
moment  with  a  magnitude  of  110336.10  In-lb/in.  with  the  use  of 
CL90  soil  for  the  embankment.  Although  the  soil  type  of  the 
embankment  has  considerable  contribution  to  this  diference,  one 
can  still  see  the  dramatic  rate  of  Increase  in  the  moment 
developed  at  the  crown  beginning  at  load  increment  7  when  the 
first  section  of  the  rock  cover  is  Installed.  Although  the 
3-layer  type  rock  cover  produces  a  much  higher  fill  height,  the 
moment  developed  at  the  crown  of  the  structure  is  significantly 
lower  than  that  generated  by  the  5-layer  type. 

Overall,  the  results  discussed  in  this  section  are  quite 
consistent.  Observing  that  the  moment  developed  at  the 
structural  crown  is  less  when  the  3-layer  type  rock  cover  is 
installed  at  a  fill  height  of  182  ft.  But  since  layer  type  and 
embankment  height  are  coupled,  the  5-layer  type  rock  cover 
proves  to  produce  the  least  amount  of  bending  moment  at  the 
crown  of  the  structure  at  the  completion  of  rock  cover 
installation.  The  least  amount  of  moment  is  developed  then  the 
5-layer  rock  cover  is  installed  in  combination  with  SM90  soil 
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type  for  the  embankment. 

The  thrust  developed  at  the  crown  of  the  structure  Is  also 
dependent  upon  the  type  of  Installation  that  the  structure  Is 
built  under  (refer  to  Figure  4.11  thru  4.14).  This  section 
describes  thrust  In  relation  to  the  type  of  rock  cover  Installed 
on  the  system.  The  comparisons  have  been  made  between  only  two 
cases  wherein  the  structure  withstood  the  entire  loading  scheme, 
these  Include  Case  IIB  and  I2B.  The  loading  scheme  for  Case  IIA 
with  a  3-layer  rock  embankment  and  Case  I2A  cused  structural 
failure  after  eight  and  four  step  load  Increments,  respectively. 

Thrust  developed  at  the  crown  of  the  structure  with  a 
30-degree  slope  embankment  and  SM85  soil  type  bedding  (Case  IIB) 
was  minimal  with  the  employment  of  the  3-layer  rock  cover. 

The  thrust  ranged  from  651  Ib/ln  when  the  embankment  Is  of  CL90 
type  soil  to  as  low  as  -531  Ib/ln  when  the  embankment  Is  of  SC90 
type  soil.  As  for  the  5-layer  type  rock  cover,  the  values  of 
thrust  ranged  from  -1664  Ib/ln  when  the  embankment  is  of  SC90 
type  soil  to  as  high  as  -1810  Ib/in  when  the  embankment  is  of 
CL90  type  solid.  Again,  even  though  the  fill  height  generated 
by  the  3-layer  rock  embankment  Is  much  higher  than  the  5-layer 
type,  it  Is  observed  that  the  final  thrust  magnitude  developed 
at  the  crown  is  less  with  the  3-layer  type  rock  cover. 

Case  I2B  shows  a  totally  contradicting  set  of  values.  At 
the  completion  of  both  types  of  rock  covers,  final  thrust 
magnitudes  developed  at  the  crown  ranged  from  at  leaslt  -1413 
Ib/in  when  the  rock  cover  type  is  5-layer  and  the  embankment  is 
SM90  type  soil  to  as  high  as  1927  Ib/in  when  the  rock  cover  type 
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Figure  M.ll  Crov/n  Thrust  vs.  Fill  !leif»ht  for  Case  IIA 


sri85  BEDDING,  30  DEG  SLOPE 


008003  00'009l  00'00ZI  00'008  00‘00t’  00'0 


Ni/ai  'mo<^3  iu  isndHi 

22-56 


Fi^'ure  M.12  Cj'own  Thrust  vs 


CL85  BEDDING,  45  DEG  SLOPE 


Ni/an  'NMOdo  lb  isndHi 


22-57 


Figure  4.13  Crown  Thrust  vs.  Fill  Height  for  Case  I2A 


SM85  BEDDING,  45  DEG  SLOPE 


\  \ 


\  \  \ 

\  \ 

H  \\X  \ 

W  '-^A 
\\^\\\ 


^\v\\ 


-4' 


w 


'4 


'^WV 


^  4J  ^  ^ 

c  c  c  c  c  c 

(U  (U  (U  <1^  O  (1> 

C  E  £  £  c  s 

u.  u  u  :£  u  is. 

c  s:  c  £  =  c 

R3  (0  <Q  (Q  lO  fQ 

js  jQ  i:  X} 

E  E  E  E  E  E 

(1)  (U  Q)  <U  Q> 

c  o  o  o  o  o 

O'  O'  O'  O'  O'  O' 

r  o  J  z  o  o 

w  c/5  u  o)  (n  o 


(A  (/}  lA  U)  cn  U) 

U  U  Li  u  U  U 

(U  C/  CU  0/  o 

<C  rt3  <C  <0  itJ  rtj 


r>  n  1/^  lO  in 

I  I  I  I  I  I 

n  o  +  X  0 


n'V'on 


v\xA 

\'v^x 


00-0083 


00-0091 


00'00ZI 


00-008 


r 

'V.v  ' 


00-00^ 


NI/81  'NM083  lU  ISndHl 


22-53 


rip>ure  M.JM  Crown  Thrust  vs 


>  I 


»  4 

I  ' 
1  • 


1 


i. 


is  3-layer  and  the  embankment  is  SM90  type  soil.  The 
interesting  part  of  this  case  that  when  the  rock  cover  is  of  the 
3-layer  type,  and  the  embankment  is  of  the  SC90  type  soil,  the 
final  thrust  magnitude  developed  at  the  crown  is  1669  Ib/in,  the 
second  lowest  value  for  thrust  found  in  this  case.  It  is, 
therefore,  difficult  to  make  any  conclusions  about  the  Influence 
of  the  rock  cover  type  on  the  thrust  for  this  particular  case. 

2.  The  Influence  of  the  Slope  of  Embankment 

To  determine  the  effect  of  the  embankment  slope  on  the 
deflection  of  the  crown  of  the  structure,  certain  variables  had 
to  be  held  constant  to  be  able  to  easily  focus  on  the  effect  of 
the  slope.  The  main  variable  that  was  held  constant  in  all 
cases  of  the  analysis  is  the  number  of  layers  used  for  the  rock 
cover,  in  this  case  the  3-layer  type,  is  selected.  All  other 
variables  such  as  the  bedding  and  soil  embankment  conditions 
were  used,  but  values  for  crown  displacement  were  compared 
individually  for  every  case. 

The  first  three  cases  studied  in  this  discussion  is  a 
system  with  CL85  bedding  and  combination  with  SM90,  SC90,  and 
CL90  soil  type  embankments  designated  as  Case  IIIA  (refer  to 
Figure  4.15).  Unfortunately,  the  structural  model  failed  after 
eight  load  increments  for  all  cases  with  a  30-degree  slope 
embankment.  Comparison  between  the  effect  of  the  embankment 
slope  was  therefore  limited  to  the  values  of  the  crown 
displacement  after  installation  of  step  load  increment  6,  the 
completion  of  the  concrete  slab.  In  the  first  case,  where  the 
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Figure  'i.JSa  Relative  Crown  Deflection  vs.  Fill  Height  for  Case  IIIA 
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Figure  'i.iSb  Relative  Crown  Deflection  vs.  Fill  Height  for  Case  IIIA 
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embankment  Is  made  up  of  SM90  type  soil,  maximum  difference 
between  crown  deflection  Is  only  0.13  Inch  after  six  load 
Increments.  The  cases  where  the  embankment  Is  made  up  of  SC90 
and  CL90  type  soil,  the  Influence  of  the  embankment  slope  Is 
much  more  apparent.  The  difference  between  crown  deflection 
after  six  load  increments  was  found  to  be  0.33  Inch  In  both 
cases  < 

The  next  three  cases  examined  In  this  section  Is  a 
structure  modeled  on  a  bedding  condition  using  SM85  type  soil  In 
combination  with  SM90,SC90,  and  CL90  type  soil  for  the 
embankment  designated  as  Case  II2B  (refer  to  Figure  4.16).  All 
structural  models  proved  to  withstand  the  entire  load  scheme  for 
these  cases,  therefore,  producing  a  complete  picture  of  crown 
deflection  with  both  the  30  and  45  degree  slope  embankment.  For 
the  model  with 

the  SM90  soil  embankment,  the  system  with  the  30-degree 
embankment  produced  greater  crown  deflectln  than  Its 
counterpart.  The  difference  between  the  two  was  0.33  inch  at 
the  end  of  nine  load  step  Increments.  For  the  model  with  SC90 
soil  embankment,  a  better  distinction  between  the  crown 
deflection  for  both  slope  conditions.  In  this  case,  the 
30-degree  slope  embankment  causes  a  crown  defelction  of  1.01 
inches  more  than  the  45-degree  slope  embankment  at  the  end  of 
nine  step  load  ""•ements.  Finally,  for  the  model  with  CL90 
soil  embankmen>.,  .  embankment  slope  of  45-degrees  produces  less 
crown  deflection  than  the  30-degree  embankment  slope,  the 
difference  being  at  least  0.86  inch. 
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In  less  crown  deflection  for  systems  with  SM85  soil  for  beddln^ 
conditions,  and  minimal  when  combined  with  SC90  soil  for  the 
embankment  for  all  cases  with  a  3**layer  rock  cover. 

A  procedure  similar  to  the  previous  one  is  used  in  the 
analysis  of  the  influence  of  the  embankment  slope  on  the  bending 
moment  developed  at  the  crown  of  the  structure.  That  is,  the 
installation  type  is  held  constant,  here  taken  to  be  the  3-layer 
type  rock  cover,  to  easily  focus  on  the  effect  of  the  slope  and 
only  the  slope  on  the  moment  developed  at  the  crown.  All  other 
variables  are  taken  into  account,  though,  to  assure  consisten'^y 
of  the  findings.  Two  bedding  conditions,  CL85  and  SM85,  plus 
three  embankment  conditions:  SM90,  SC90,  CL90;  are  included  in 
the  analysis. 

The  first  three  cases  examined  is  a  system  w^th  CL85  soil 
bedding  in  combination  with  SM90,  SC90,  and  CL90  soil 
embankments  (refer  to  Figure  A. 17).  The  structural  models  failed 
after  the  completion  eight  load  increments  for  all  cases  with 
the  30-degree  slope  embankment.  Then  again,  the  comparison 
between  the  two  slope  conditions  for  the  moment  developed  at  the 
crown  is  limited  to  the  values  of  crown  moment  developed  after 
six  load  increments.'  The  first  case,  where  the  embankment  is 
made  up  of  SM90  type  soil,  shows  that  by  installing  a  30-degree 
slope  embankment,  the  moment  at  the  crown  can  be  reduced  by  as 
much  as  10647  in-lb/in,  less  than  that  developed  with  a 
45-degree  slope  embankment  at  the  end  of  six  load  increments. 

The  same  is  true  with  the  next  case,  where  the  embankment  is 
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Figure  M.17b  Crown  Moment  vs.  Fill  Height  for  Case  IIIB 
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made  up  of  SC90  type  soil.  The  noaent  developed  at  the  crown  Is 
4528  In-lb/ln  less  with  a  30-degree  slope  embankment  Installed 
than  with  a  45-degree  slope  embankment.  Finally  with  an 
embankment  type  soil  of  CL90  used,  a  30-degree  slope  embankment 
produces’  5939  in-lb/in  less  moment  than  a  45-degree  slope 
embankment . 

The  next  three  cases  examined  Is  a  system  with  SM85  soil 
bedding  in  combination  with  SM90,  and  CL90  soil  embankments 
(refer  to  Figure  4.18).  For  the  model  with  SM90  type  soil  for 
the  embankment  material,  the  bending  moment  developed  at  the 
crown  of  the  structure  was  found  to  be  greater  when  a  45-degree 
slope  embankment  Is  Installed.  The  difference  In  magnitude 
being  7705  In-lb/in  at  the  end  of  four  load  increments  up  to 
27364  In-lb/ln  at  the  completion  of  nine  load  Increments. 

Results  from  the  tests  using  the  structure  model  with  SC90  type 
embankment  soil  are  consistent  with  the  previous  one  but  less 
the  differences  In  the  moment  are  less  dramatic.  Results  show 
that  again,  by  using  a  30-degree  slope  embankment,  the  moment 
developed  at  the  crown  Is  less  than  that  with  a  45-degree 
embankment.  The  differences  being  4919  in-lb/ln  at  the  end  of 
four  load  Increments  up  to  13344  in-lb/in  at  the  completion  of 
nine  load  Increments.  Finally,  with  the  case  of  the  system 
consisting  of  an  embankment  made  up  of  CL90  type  soil.  It  was 
also  found  that  the  embankment  with  a  30-degree  slope  produces 
less  moment  at  the  crown.  This  is  the  case  wherein  "he 
influence  of  the  slope  on  the  moment  is  most  apparent.  At  the 
end  of  nine  load  increments,  the  crown  moment  with  a  45-degree 
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slope  embankoent  was  found  to  be  156650  in-lb/ln,  nearly 

665  in-lb/ln  greater  than  the  moment  developed  at  the  crown  with 

a  30-degree  slope  embankment. 

The  thrust  developed  at  the  crown  of  the  structure  Is  also 
highly  influenced  by  the  slope  of  the  embankment.  Again,  to 
make  the  analysis  simpler  and  to  be  able  to  focus  on  the 
Influence  of  only  the  slope,  the  installation  type  is  held 
constant  for  all  cases  and  taken  to  be  the  3-layer  type  rock 
cover.  The  bedding  types  Included  in  the  analysis  are  SM85  and 
CL85,  and  emoankment  types  used  are  SM90,  SC90,  and  CL90. 

The  first  three  cases  analyzed  are  systems  in  CL85  bedding 
(refer  to  Figure  4.19).  This  condition  though,  structural 
models  failed  after  eight  load  increments  for  systems  with 
30-degree  slope  embankments.  Therefore,  like  the  former 
analyses  made  in  this  section  the  observations  will  be  limited 
to  the  crown  thrust  after  six  load  Increments  for  each  case. 

Even  though  the  45-degree  condition  is  not  feasible  in  practice, 
analyses  is  still  performed  to  determine  the  consistency  of  the 
results.  The  first  case  analyzed  is  the  system  constructed  in  a 
SM90  soil  embankment.  Here,  the  results  show  that  an  embankment 
with  a  45-degree  slope  produces  less  thrust  at  the  crown  than  an 
embankment  with  a  30-degree  slop-J.  It  was  found  that  at  the  end 
of  six  construction  Increments,  a  45-degree  slope  produces  less 
thrust  at  the  crown  than  the  30-degree  slope.  The  difference 
being  269  Ib/in,  nearly  42%  smaller  with  a  steeper  slope. 
Similarly,  the  case  where  the  embankment  is  made  of  SC90  type 
soil,  the  same  results  occur.  An  embankment  with  a  45-degree 
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Figure  M,19b  Crown  Thrust  vs.  Fill  Height  for  Case  I I 1C 
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slope  produces  nearly  322  less  thrust  than  an  embankment  with  a 
30-degree  slope.  Finally,  system  with  an  embankment  of  CL90 
type  soil,  unsurprisingly  yields  the  same  results.  An 
embankment  slope  of  AS-degrees  produces  nearly  302  less  thrust 
than  the  embankment  with  a  30-degree  slope.  In  general,  for  a 
syst'm  with  CL85  soil  bedding,  an  embankment  of  SM90  type  soil 
Installed  In  a  AS-degree  slope  produces  the  least  thrust  at  the 
crown  after  the  completion  of  all  step  load  Increments. 

The  next  three  cases  examined  were  systems  wherein  the 
structure  Is  modeled  on  SM85  type  soil  bedding  In  combination 
with  SM90,  SC90,  and  CL90  type  soil  for  the  embankment  (refer  to 
Figure  A. 20).  In  all  three  cases  It  was  found  that  at  the  end 
of  six  step  load  Increments,  the  system  with  the  30-degree  slope 

J 

embankment  produced  more  thrust  at  the  crown  than  the  A5-degree 
option,  the  average  difference  of  the  thrust  developed  at  the 
crown  being  2A2  Ib/ln.  But  at  the  completion  of  all  nine  step 
load  Increments,  It  was  found  that  an  embankment  with  a 
30-degree  slope  shall  produce  a  final  crown  thrust  magnitude  up 
to  at  least  1862  to  at  most  2182.  With  SM90  soil 
embankment,  the  difference  of  the  final  thrust  magnitude 
produced  at  the  crown  between  the  30-degree  and  A5-degree  slope 
Is  1320  Ib/ln  favoring  the  former.  With  SC90  soil  embankment, 
the  difference  was  found  to  be  1212  Ib/ln.  The  system  producing 
the  least  thrust  at  the  end  of  all  step  load  Increments  is  the 
one  with  the  SC90  soil  for  the  embankment,  with  a  final  crown 
thrust  magnitude  of  531  Ib/ln. 
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3.  The  Influence  of  Soil  Type  of  Embankment 

Crown  deflection  Is  greatly  influenced  by  soil  type, 
meaning  that  the  physical  properties  of  the  soil,  not  only  the 
density,  plays  an  Important  part  in  the  determination  of  load 
transfer  from  soil  to  structure.  This  section  discusses  the 
results  found  on  the  Influence  of  soil  type  on  the  relative 
deflection  of  the  crown.  To  minimize  the  number  of  cases  to  be 
examined  and  focus  the  attention  on  the  Influence  of  only  the 
embankment  soil  on  the  structure,  the  system  models  were  set  to 
have  a  constant  embankment  slope,  taken  to  be  30’ degrees 
designated  as  Case  IIIl.  All  other  variables  such  as  bedding 
type  and  the  number  of  layers  for  the  rock  cover  are  taken  into 
account  (refer  to  Figure  4.21). 

The  first  case  studied  is  a  system  model  on  CL85  soil 
bedding  with  a  3->layer  rock  cover.  Although  it  was  found  that 
at  the  end  of  six  load  Increments,  crown  relative  settlement  had 
a  standard  deviation  of  0.07  in.  with  the  minimum  of  8.38  in. 
occurlng  when  the  embankment  is  of  SM90  type  to  maximum  of  8.55 
in.  occurlng  when  the  embankment  is  of  SC90  type,  the  structural 
models  failed.  Meaning  that  these  types  of  embankment  soil 
Installed  within  these  conditions  is  not  feasible  for  -  practical 
purposes . 

The  second  case  studied  is  a  system  similar  to  the  one 
described  above  only  that  this  is  modeled  with  a  5-layer  rock 
cover.  Relative  crown  deflection  at  the  completion  of  all  nine 
load  increments  had  a  standard  deviation  of  0.048  in.  with  the 
minimum  being  10.58  in.  occurlng  when  the  embankment  is  of  SM90 
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Figure  u.21d  Relative  Crown  Deflection  vs.  Fill  Height  for  Case  IIH 


type  to  a  aaxlnun  of  !.0.69  In.  occurlng  when  the  embankment  Is 
of  CL90  type.  In  this  case.  It  could  be  said  that  the  soil  type 
for  the  embankment  does  not  cause  extreme  difference  in  the 
relative  deflection  of  the  crown. 

The  third  case  is  a  system  modeled  on  SM85  soil  bedding 
with  a  3-layer  rock  cover  Installed.  In  this  case,  a  system 
model  with  the  SM90  soil  embankment  caused  the  least  amount  of 
relative  crown  deflection,  being  5.84  in.  The  Installation  of 
SC90  soil  for  the  embankment  causes  the  crown  to  deflect  by  6.J4 
in.,  half-inch  more  than  with  the  previous  condition.  In  this 
case  overall,  an  SM90  type  soil  used  for  the  embankment  material 
shall  produce  the  least  crown  deflection. 

The  fourth  and  final  case  was  found  to  be  the  most  ideal 
among  all  the  cases  studied.  The  minimum  relative  crown 
deflection  was  found  to  be  4.76  in.  and  attained  by  installing 
an  SM90  soil  type  embankment  under  a  5-layer  rock-cover  over  an 
SM85  soil  bedding.  The  maximum  is  found  to  be  4.93  in.  with  an 
SC90  soil  embankment. 

Overall,  by  examining  all  the  observations  mentioned  above, 
it  was  found  that  an  SM90  type  soil  used  for  the  embankment 
causes  the  least  amount  of  crown  relative  deflection. 

Momei.ts  developed  at  crown  are  also  highly  influenced  by 
the  soil  type  of  the  embankment.  In  this  discussion,  similar  to 
the  previous  one,  the  system  models  were  set  to  have  a  constant 
embankment  slope  of  30  degrees  for  the  same  reason  as  to  focus 
the  attention  on  the  influence  of  ‘he  embankment  soil  only. 

Also,  all  other  variables  such  as  bedding  type  and  the  number  of 
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layers  for  the  rock  cover  are  taken  Into  account  to  visualize 
the  consistency  of  the  results  (refer  to  Figure  4.22). 

Case  number  III2.a  involves  a  system  model  on  CL85  soil 
bedding  with  3  layers  of  rock  as  the  cover.  The  moment  at  the 
crown  of  the  structure  at  the  completion  of  six  load  increments 
was  governed  by  the  CL90  embankment  producing  a  minimum  moment 
magnitude  of  80039  in-lb/ln.  The  maximum  being  85160  in~lb/in 
produced  by  SC90  embankment.  This  analysis  covers  crown  moment 
after  six  load  Increments  only  because  a  system  under  these 
conditions  failed  after  eight  load  increments. 

The  next  case,  Case  III2.b  involves  a  system  modeled  in  the 
same  manner  as  the  previous  one  discussed  except  that  a  5-layer 
rock  cover  was  Installed.  The  structure  models  withstood  the 
entire  load  scheme  for  all  nine  load  increments  providing 
complete  data  for  analysis.  The  bending  moments  occurlng  at  the 
crown  had  a  standard  deviation  of  2211  in-*lb/in  with  a  maximum 
of  113096  In-lb/ln  occurlng  when  the  embankment  is  of  SC90  soli 
and  a  minimum  of  17756  In-lb/in  occurlng  when  the  embankment  is 
of  CL90  soil. 

The  third  case,  Case  1112. c  involves  a  system  modeled  on  an 
SM85  bedding  with  a  3-layer  rock  cover  installed.  In  this  case, 
the  -embankment  made  up  of  SM90  soil  produced  the  least  moment 
magnitude  among  all  three  embankment  soil  types.  At  the 
completion  of  nine  load  increments,  the  SM90  type  soil  produr-d 
a  crown  moment  magnitude  of  83773  in-lb/in,  a  magnitude  that  is 
20X  less  than  that  produced  by  SC90  type  soil,  which  in  turn 
produced  the  greatest  crown  moment. 
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Figure  ‘l.?2a  Crown  Moment  vs.  Fill  Height  for  Case  III2 
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The  fourth  and  final  case,  Case  III2.d  is  similar  to  the 
previous  one  except  that  a  5-layer  rock  cover  is  Installed.  And 
similarly,  the  SM90  soil  produced  the  least  amont  of  moment  at 
the  completion  of  nine  load  increments.  A  standard  deviation  of 
A958.071  in-lb/ln  with  a  maximum  of  90750  in-lb/in  occurlng  when 
the  embankment  is  of  SC90  type  soil,  and  a  minimum  of  78640 
in-lb/in  when  the  embankment  type  is  SM90. 

Thrust  developed  at  the  crown,  like  the  other  types  of 
structural  reactions  are  also  highly  dependent  on  the  kind  of 
soil  used  for  the  embankment.  This  discussion  shall  also  be 
limited  to  system  models  with  30-degree  slope  embankments, 
taking  into  account  all  the  other  factors  such  as  bedding 
conditions  and  the  number  of  layers  for  the  rock  cover, 
designated  as  Case  III3  (refer  to  Figure  4.23). 

The  first  case  to  be  analyzed  is  a  structure  modeled  on  CL8 
under  a  3-layer  rock  cover.  Unfortunately,  the  structural 
models  for  this  system  failed  after  eight  load  Increments. 
Therefore,  the  analysis  is  limited  to  the  findings  after  the 
completion  of  six  load  increments.  At  the  end  of  the  fourth 
load  increment,  it  was  found  that  the  CL90  embankment  soil 
produced  less  crown  thrust  magnitude  than  the  SM90  type.  But  at 
the  completion  of  six  load  increments,  it  was  found  that  the 
SM90  type  soil  produced  less  thrust  at  the  crown  than  the  CL90 
type,  which  in  turn  produced  the  highest  thrust  magnitude 
between  all  three  typts.  The  SC90  type  was  consistent 
throughout  the  entire  scheme  in  producing  the  least  amount  of 
thrust  at  the  crown,  the  final  thrust  magnitude  after  six  load 
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Figure  4.?3a  Ci'ov/n  Thrust  vs.  Fill  Height  for  Case  ITI3 
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increments  being  R87  Ib/in. 

A  different  set  of  results  appears  when  the  rock  cover  of 
the  above  case  is  replaced  by  the  5-layer  type.  At  the  end  of 
four  load  increments,  the  SM90  type  soil  produced  greater 
thrust  in  magnitude  at  the  crown  than  the  SC90  type.  But  after 
six  load  increments,  it  was  found  that  the  SM9f  type  soil 
produced  less  thrust  than  the  SC90  type.  This  condition  went  on 
until  the  completion  of  all  nine  load  increments,  making  the 
SM90  type  the  ideal  soil  type  to  use  in  this  case,  producing  a 
final  crown  thrust  of  1384  Ib/ln,  21  less  than  that  produced  by 
the  SC90  type.  The  SC90  type  was  consistent  in  producing  the 
largest  amount  of  thrust  throughout  the  entire  load  Increments 
producing  a  final  crown  thrust  of  1540  Ib/ln,  11.321  larger  than 
the  minimum. 

The  most  interesting  case  in  this  analysis  Is  the  case 
wherein  the  system  is  modeled  on  an  SM85  bedding  under  the 
3-layer  rock  cover.  The  thrust  developed  at  the  crown  Increased 
directly  proportional  to  the  fill  height  for  all  types  of 
embankment  soil  until  the  completion  of  the  sixth  load 
Increment.  The  thrust  values  then  decreased  even  though  the 
fill  height  continued  to  Increase  at  the  beginning  of  the 
seventh  load  increment.  The  final  thrust  developed  at  the  crown 
per  soil  type  was  therefore  the  least  among  all  the  cases 
studied  in  this  section.  The  SC90  type  produced  the  least 
amount  of  thrust  at  the  crown  having  a  magnitude  of  531  Ib/in. 
Followed  by  embankment  soil  type  SM90  producing  a  final  thrust 
magnitude  of  607  Ib/in  and  then  by  the  CL90  type  producing  651 
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The  final  case  studied  in  this  section  Is  a  system  modeled 
on  SM85  bedding  under  a  5-layer  rock  cover.  At  the  end  of 
fourth  load  increment,  the  SM90  soil  type  produced  the  greatest 
amount  of  thrust  at  the  crown.  But  after  the  completion  of  the 
sixth  load  Increment,  the  CL90  soil  type  produced  the  highest 
thrust  magnitude  among  all  three  embankment  types.  At  the 
completion  of  the  final  load  increment,  soil  type  CL90  was  found 
to  generate  the  largest  amount  of  thrust  at  the  crown  of  the 
structure.  Its  magnitude  being  1810  Ib/ln.  The  SC^^O  type  soil 
was  consistent  in  producing  the  least  amount  of  thrust  on  the 
crown  throughout  the  whole  loading  scheme,  generating  a  final 
thrust  magnitude  after  ninth  load  increment  of  1664  Ib/in. 
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SUMMARY,  CONCLUSIONS  AND  RECOMMENDATIONS 

A-  SUMMARY  AND  CONCLUSIONS 

This  report  presents  an  Investigation  of  the  behavior  of  a 
reinforced  concrete  arch  structure  under  various  embankment  and 
bedding  conditions.  A  finite  element  analysis  of  structural 
responses  was  performed  on  system  models  with  varying  bedding, 
embankment  slope,  embankment  soil  type,  and  rock  cover  type.  The 
findings  and  conclusions  from  this  investigation  are  as  follows: 

1.  The  results  of  the  finite  element  model  study  showed 
consistency  and  reliability  in  predicting  soli  and  structure 
responses . 

2.  The  structures  that  were  predicted  to  fall  were  the 
ones  systemed  in  the  order  of: 

a.  Having  a  CL85  (clay)  soil  type  bedding,  with  a 
3-layer  rock  cover,  and  a  30-degree  slope  embankment. 

A  structure  modeled  under  these  conditions  was 
predicted  to  fall  no  matter  what  type  of  soil  was  used 
for  the  embankment. 

b.  Having  a  CL85  (clay)  soil  type  bedding,  with 
a  5-layer  rock  cover,  and  a  45-degree  slope 
embankment.  A  structure  modeled  under  these 
conditions  was  predicted  to  fail  no  matter  what  type 
of  soil  was  used  for  the  embankment. 

3.  Between  the  two  types  of  bedding  materials  used  in  the 
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tests,  It  was  found  that  at  a  relative  compaction  of  85*,  a 
silty  sand  type  of  soil  (type  SM85)  produced  less  structural 
settlement  than  a  clayey  type  soil  (CL85). 

4.  An  embankment  constructed  with  a  AS-degree  slope  caused 
less  moments  to  develop  around  the  structure  than  an  embankment 
with  a  30-degree  slope.  In  addition  the  former  condition  also 
caused  structural  reaction  to  develop  more  in  form  of  axial 
forces  (thrust)  rather  than  bending  moments,  which  is  preferable 
in  shell  structures. 

5.  An  embankment  constructed  with  a  silty  sand  type  (SM85) 
of  soil  produces  less  moment  around  the  structure  and  encourages 
structural  reaction  to  develop  through  axial  forces  (thrust). 

6.  A  rock-cover  composed  of  a  mean  boulder  diameter  of  14 
Inches  and  set  in  five  layers  causes  less  magnitude  In 
structural  reponse  at  the  completion  of  the  embankment  than  a 
3-layer  rock  cover  with  a  mean  boulder  diamet'er  of  42  inches. 

7.  Overall,  a  structure  set  on  an  SM85  sand  soil  type 
bedding,  with  SM90  soil  embankment,  5-layer  type  rock  cover,  and 
set  on  a  30-degree  slope  proves  to  be  the  most  effective  system. 
The  finite  element  analysis  show  that  an  arch  structure 
constructed  in  this  manner  will  produce  the  least  amount  of 
bending  moment,  the  load  transfer  onto  the  structure  being 
concentrated  mainly  on  axial  thrust.  Figures  4.24  thru  4.26  show 
the  final  soil  pressure,  moment,  and  thrust  distribution  for  a 
structure  modelled  in  this  manner. 
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B.  RECOMMENDATIONS 


Recommendations  for  research  based  on  this  study  are  as 
follows : 

The  results  of  the  finite  element  model  study  presented  in 
this  report  are  based  on  the  preliminary  investigation  for  the 
feasibility  study  of  an  arch  structure  in  a  centrigfuge.  A  more 
comprehensive  (extensive)  research  including  embankments  of 
other  sizes  and  geometries,  structural  re  inf ormcement ,  different 
soil  types,  compaction  for  embankments,  and  bedding  parameters 
is  required  in  order  to  gain  a  better  understanding  of  the 
relationship  of  the  arch  structure  and  the  embankment. 

In  order  to  achieve  a  comprehensive  model  study  using  the 
centrifuge  model  technique,  an  on-board  data  acquisition  system 
is  required.  Past  experience  of  the  researchers  prove  that  this 
is  a  necessity  for  acquiring  reliable  and  qualitative  data. 

A  field  study  is  necessary  to  verify  the  results  of  the 
centrifuge  model  study  and  finite  element  analysis  once 
reasonable  information  is  obtained. 
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ABSTRACT 


Protective  military  structures  are  typically  con¬ 
structed  o£  massive,  monolithic  concrete  slabs.  This 
practice  is  considered  necesseury  in  order  to  protect 
personnel  and/or  vital  equipment  within  the  structure 
from  the  harmful  blast  effects  of  conventional  weaponry. 

A  high  intensity  compressive  longitudinal  blast  wave 
incident  upon  a  protective  shelter  propagates  through  the 
monolithic  concrete  shelter  wall,  and  reflects  as  a  tensile 
wave  at  the  Interior  face  of  the  wall.  The  reflected 
tensile  wave  can  cause  extensive  spalling  of  the  interior 
face  of  the  shelter  walls,  thus  posing  a  potentially 
dangerous  threat  to  personnel  and  equipment  in  the  shelter. 

A  conprehensive  numerical  analysis  is  conducted  to 
evaluate  the  effectiveness  of  various  "layered  structures" 
for  enhanced  protection  against  blast  effects  from  con¬ 
ventional  weaponry.  Of  primary  concern  is  the  capability 
of  layered  systems  to  mitigate  stress  wave  propagation, 
and  the  subsequent  elimination  of  interior  spalling. 

The  finite  element  me^od  (FEM)  is  employed  to 
investigate  the  wave  propagation  characteristics  of  three 
layered  systems;  1)  concrete-polystyrene-concrete  (CPC), 

2)  sand-polystyrene-concrete  (SPC) ,  and  3)  sand- 
polystyrene-concrete-polystyrene-concrete  (SPCPC) .  Both 
linear  and  nonlinear  material  behavior  are  considered. 
Stress  transmission  ratios  are  determined  for  each  layered 
system. 
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The  results  o£  the  investigation  indicate  that 
layered  media  significantly  reduce  stress  transmission, 
and  that  the  application  of  layered  media  to  protective 
shelters  can  eliminate  interior  spalling. 
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I .  INTRODUCTION 


Motivation 

Protective  military  shelters  are  typically 
constructed  of  massive,  monolithic  concrete  slabs.  This 
practice,  which  has  not  changed  significantly  since  the 
end  of  the  Second  World  War,  is  necessary  in  order  to 
protect  personnel  and/or  vital  equipment  within  the 
structure  from  the  harmful  effects  of  conventional 
weaponry.  The  drastic  devaluation  of  the  dollar  in  the 
NATO  nations  during  the  past  several  years  has  made 
construction  costs  for  these  types  of  structures 
prohibitive.  This  fact  dramatically  underscores  the 
urgency  for  the  development  of  alternative,  cost-effective 
protective  structure  designs. 

When  a  typical  monolithic  concrete  protective  shelter 
is  sxibject  to  an  impulse  loading  (blast  wave)  from 
conventional  weaponry,  a  high-intensity  stress  wave 
propagates  through  the  shelter  wall.  The  stress  wave 
reflects  from  the  interior  face  of  the  shelter  wall  as  a 
tension  wave  and  can  cause  extensive  spalling  of  the 
interior  face  of  the  wall,  spalled  concrete  fragments 

are  hurled  away  from  the  wall  at  high  velocities,  posing  a 
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potentially  dangerous  threat  to  personnel  and  equipment 
within  the  shelter. 

Recent  experimental  studies  (5,26)  and  analytical 
investigations  conducted  by  Tedesco  (20,21),  Tedesco  et  al 
(23,24),  and  Landis  (14)  have  indicated  that  "layered  ' 
structures"  may  provide  a  viable  alternative  to 
conventionally  hardened  structures.  The  intent  of  this 
investigation  is  to  evaluate  the  effectiveness  of  various 
"layered  structures"  for  enhanced  protection  against 
impulse  loading  from  conventional  weaponry. 

Previous  Research 

Extensive  information  is  available  concerning  the 
propagation  of  stress  waves  in  elastic,  isotropic  media. 
Information  is  also  available  concerning  the  propagation 
of  stress  waves  in  inelastic  media.  The  majority  of  the 
information,  however,  deals  with  stress  wave  propagation 
in  homogeneov.s  metal  wires  and  bars  of  uniform  cross- 
section. 

Little  information  is  available  on  the  propagation  of 
stress  waves  through  layered  media.  Even  less  information 
exists  on  the  use  of  layered  media  in  protective  shelters 
to  resist  high-intensity  blast  loadings.  Recent  research, 
however,  indicates  that  layered  structures  may  provide  a 
viable  alternative  to  conventionally  hardened  structures,.. 

Analytical  investigations  conducted  by  Tedesco 
(20,21)  indicate  that  layered  systems  may  be  utilized  to 
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alter  the  stress  wave  propagation  characteristics  of  a 
protective  shelter.  The  investigations  revealed  that 
certain  combinations  of  materials,  constructed  in  layers, 
significantly  reduce  the  stress  transmitted  through  the 
shelter  walls.  It  was  determined  that  the  combination  of 
materials  required  to  substantially  reduce  the  transmitted 
stress  consists  of  a  low-density,  low-strength  center 
layer  sandwiched  between  two  high-density,  high-strength 
layers. 

Subsequent  analytical  invest! gation^j  conducted  by 
Tedesco  et  al  (23,24)  examined  three  layered  systems,  each 
consisting  of  a  center  absorption  layer  surrounded  by  two 
layers  of  concrete.  The  absorption  layers  examined 
included  air,  polystyrene,  and  soil.  The  results  from 
these  preliminary  studies  revealed  that  each  layered 
system  significantly  reduced  the  stress  transmitted  to  the 
interior  of  the  structure. 

Analytical  investigations  conducted  by  Landis  (14) 
examined  layered  systems  consisting  of  an  absorption  layer 
of  polystyrene  or  soil,  sandwiched  between  two  layers  of 
concrete.  Both  layered  systems  significantly  reduced  the 
stress  transmitted  through  the  structure..  The 
polystyrene,  however,  proved  to  be  a  much  more  effective 
absorption  layer  than  the  soil,. 

The  research  conducted  in  the  aforementioned  studies 
investigated  only  the  linear-elastic  behavior  of  the 
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layered  systems.  The  effect  of  nonlinear  behavior  upon 
stress  wave  propagation  in  layered  media  has  not  yet  been 
examined. 


Oblective 

The  primary  objective  of  this  study  is  to  evaluate 
the  effectiveness  of  layered  systems  to  mitigate  stress 
wave  propagation  via  a  comprehensive  numerical  analysis. 
The  numerical  analyses  are  conducted  on  several  different 
types  of  layered  systems  subject  to  impulse  loading  in 
order  to  ascertain  the  characteristics  of  wave  propagation 
through  the  various  material  layers  comprising  the 
systems.  From  the  results  of  the  numerical  analyses, 
recommendations  for  layered  systems  as  an  alternative  to 
conventionally  hardened  protective  structures  are 
presented. 


Scope  of  Study 

The  finite  element  method  (FEM)  is  utilized  in  the 
numerical  analyses  to  investigate  the  wave  propagation 
characteristics  of  several  layered  systems.  The  layered 
systems  investigated  include  a  concrete-polystyrene- 
concrete  system,  a  sand-polystyrene-ccncrete  system,  and  a 
sand-polystyrene-concrete-polystyrene-concrete  system. 

The  dimensions  of  the  concrete  and  sand  layers  are 
constant  for  each  layered  system.  Polystyrene  layers  of 
varying  thickness  are  examined  in  order  to  properly  assess 
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the  effect  of  eOasorption  layer  thickness  upon  wave 
propagation.  Both  linear  and  nonlinear  material  behavior 
are  considered  in  the  wave  propagation  analysis.  The 
ratio  of  transmitted  stress  to  incident  stress 
(transmission  ratio)  is  determined  for  each  system 
investigated. 
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II.  STRESS  WAVE  PROPAGATION  IN  SOLIDS 


In  static  analys'es,  it  is  assumed  that  a  structure  is 
gradually  loaded  over  an  interval  of  time  such  that 
inertial  properties  have  no  effect  on  the  response  of  the 
structure.  The  structure  is  in  a  state  of  static 
equilibrium  and  time  is  not  a  variable  in  the  solution. 
Dynamic  problems,  however,  occur  when  loading  conditions 
vary  significantly  in  time  to  produce  inertial  forces  in 
the  structure.  The  structure  is  in  a  state  of  dynamic 
equilibrium,  and  the  response  of  the  structure  is  time 
dependent.  In  general,  if  the  period  of  the  forcing 
function  is  less  than  three  times  the  fundamental  period 
of  the  structure,  then  the  problem  becomes  dynamic,  and 
inertial  forces  must  be  considered  in  the  structural 
response  (8). 

When  an  elastic  body  is  subjected  to  dynamic  loading, 
stresses  and  strains  propagate  throughout  the  structure  at 
a  finite  velocity.  In  most  structural  dynamic  problems, 
the  rate  of  application  of  the  load  is  small  in  comparison 
with  the  velocity  of  stress/strain  propagation. 

Therefore,  the  entire  structure  respo’^'^'-  to  the  loading 
condition  immediately,  and  wave  propagation  is  not  a 
consideration  in  the  solution.  In  instances  where  the 


load  Is  applied  at  a  high  rate  (as  in  the  case  of  an 
impact  or  explosion)^  however,  the  propagation  of 
stress/strain  waves  must  be  considered  in  the  problem 
solution. 


Basic  Stress  Wave  Theory 

A  stress  wave  is  defined  as  a  moving  part  of  a  medium 
in  a  state  of  stress.  Since  any  stress  has  associated 
with  £t  a  corresponding  strain,  stress  waves  and  strain 
waves  are  identical.  When  a  force  is  Initially  applied  to 
a  body,  a  stress  wave  is  generated  at  the  point  of 
contact.  The  stress  wave  propagates  away  from  the  point 
of  contact  at  approximately  the  seismic  velocity  of  that 
body.  Any  point  in  the  body  remains  unstressed  and 
unaffected  by  the  disturbance  until  reached  by  the  stress 
wave.  The  stress  wave  propagates  through  the  body, 
reflects  from  the  boundaries  of  the  body,  and  interacts 
with  other  stress  waves  to  create  very  complicated  stress 
distributions  in  the  body. 

Stress  waves  are  classified  as  normal  or  tangential. 
Normal  waves  are  analagous  to  the  normal  component  of 
stress.  Mass  particles  under  the  influence  of  a  normal 
wave  move  parallel  to  the  direction  of  propagation  of  the 
wave.  A  body  subjected  to  a  purely  normal  wave  undergoes 
a  volumetric  change  due  to  lengthening  or  shortening,  but 
does  not  undergo  any  rotation.  For  this  reason,  normal 
waves  are  referred  to  as  longitudinal  waves  or  dilational 
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waves.  Tangential  waves  are  similiar  to  the  shear 
component  of  stress.  Mass  particles  under  the  influence 
of  a  tangetial  wave  move  perpendicular  to  the  direction  of 
wave  propagation.  A  body  subjected  to  a  purely  tangential 
wave  undergoes  a  rotational  displacement  but  does  not 
experience  a  change  in  volume.  For  this  reason, 
tangential  waves  are  referred  to  as  shear  waves, 
rotational  waves,  or  distortional  waves. 

The  general  equations  governing  the  propagation  of 
stress  waves  can  be  derived  by  analyzing  a  differci^cial 
volume  SxSySz,  as  shown  in  Figure  2.1.  The  residual  force 
f  in  each  direction  can  be  found  by  summing  forces  in 
each  direction.  Summing  forces  in  the  x  direction  yields 

F  •  (a  5j;  5r  -  c  5y  5z  +  (x  +  5x  5z  -  t  3;e  5z  + 

*  \  *  dx  /  *  \  dy  y  *>' 

which  simplifies  to 

(da  dx  \ 

dx  dy  dz  J 

Similiarly, 


/CT 

da  --  •' 

F  =  1 

+  ^  +  ) 

ay  Qz  / 

5.r  dy  dz 

y 

\  ox 

( 

dx  da  \ 

F,  =  1 

+  — il  +  _i  ) 

dx  oy  dz 

\  OX 

dy  dz  / 
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Flgur*  2.1.  Stress  Components  Acting  on  a  Differential 
Volume 


Applying  Newton's  2ndi  law  of  motion,  F  •  ma ,  yields 
f  ■  p  5*  5j>  or 

f  ■  p  &x  5>»  5z  ~ 

'  di^ 

f  -  p  5;e  52  ^ 

dt^ 


(2. 


Equating  equations  (2.1)  and  (2.2)  and  simplifying  yields 


n 


■“'n 

. '  i 


!j 


(’• 

L 


a^w 

^  di^  ^y 


Equation  (2.3)  is  the  governing  equation  for  stress 
wave  propagation  in  any  media/  and  is  valid  for  any 
stress-strain  relationship.  Therefore/  the  stress  wave 
propagation  in  any  medium  is  determined  by  the 
substitution  of  the  appropriate  stress-strain  relations 
into  equation  (2.3),  and  the  concurrent  application  of  the 
appropriate  boundary  conditions. 

It  has  been  proven  (12/19)  that  equation  (2,3) 
describes  the  dilatational  and  distortional  stress  wave 
components.  It  has  also  been  shown  (12/19)  that/  in  any 
isotropic  elastic  medium/  a  stress  wave  must  propagate  at 
a  velocity  dependent  upon  the  mechanical  properties  of  the 
medium.  A  dilatational  (longitudinal)  wave  must  travel  at 
a  velocity  of 


c 


I 


•  f(l-v)  f 

.(1  v)(l  -2v)p 


(2.4) 


and  a  distortional  (rotational)  wave  must  travel  at  a 
velocity  of 

Cj  »  VX  (2.5) 

'  P 


where  E  is  Young's  modulus/  v  is  Poisson's  ratio,  p  is  the 
material  mass  density,  and  G  is  the  shear  modulus  (6). 
Moreover,  for  an  elastic  isotropic  material,  equation 
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(2.3)  can  be  further  simplified  to 


^  -  cV^w 


where  v’  represents  the  operation 
^  3*^  3y^  ' 


(2.6) 


Plane  Stress  Wave  Propagation 
Stress  waves  are  propagated  in  all  directions  in  a 
media  from  the  point  of  a  disturbance.  At  a  great  enough 
distance  from  the  disturbance,  however,  the  stress  waves 
are,  for  all  practical  purposes,  plane  waves.  In  the  case 
of  a  plane  wave,  it  can  be  assumed  that  all  mass  particles 
under  the  influence  of  the  wave  are  moving  either  parallel 
to  the  direction  of  wave  propagation  (dilatational  waves) 
or  perpendicular  to  the  direction  of  wave  propagation 
(distortional  waves)  (25). 

For  propagation  in  the  x  direction  only  (refer  to 
Figure  2.2),  the  deformation  becomes  a  functi'^n  of  only. 
Equation  (2.6)  therefore  reduces  to 
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The  general  solution  to  this  equation  is 

u{x,t)  =  f^{x+c^l)  +  f^ix-c^t)  (2.8) 

where  /J  and  are  arbitrary  functions  depending  upon 
initial  conditions  and  describe  the  shape  of  the  stress 
wave.  The  function  yj  represents  a  stress  wave  traveling 
in  the  negative  x  direction  with  constant  velocity  c,.  The 
function  yj  represents  a  stress  wave  traveling  in  the 
positive  X  direction  with  constant  velocity  c,.. 

A  segment  of  a  plane  dilatational  wave  is  illustrated 
in  Figure  2.2.  The  medium  at  section  A-A  in  Figure  2.2(a) 
experiences  a  normal  stress  After  a  time  interval  ai , 
the  segment  of  the  wave  initially  corresponding  to  section 
A-A  has  traveled  a  distance  5x  to  section  B-B  (Figure 
2.2(b)).  The  stress  wave  has  imposed  an  impulse  a,  5^  on 
the  medium  between  sections  A-A  and  B-B.  Equating  this 
impulse  to  the  change  in  momentum  yields 
ot  =  p  5jr  »  p  c,  5/ 

which  reduces  to 

a  =  p  c,  K  or  K  =  (2.9) 

where  /  is  the  particle  velocity  parallel  to  the  direction 
of  wave  propagation.  Similiarly,  for  a  distortional  wave, 

T  =  p  c,  /  or  V  =  ^  (2.10) 

yt  ^  i  y  >'  pc. 
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where  K  is  the  particle  velocity  perpendicular  to  the 
direction  of  wave  propagation.  Equations  (2,9)  and  (2.10) 
Indicate  that  stress  and  particle  velocity  are  linearly 
related  by  the  material  impedance  (acoustic  resistance), 
pc 


Figure  2.2.  Section  of  Medium  Subject  to  a  Transient 
Stress  Wave 

The  stress  wave  velocity  and  the  particle  velocities 
must  not  be  construed  as  the  same.  The  stress  wave 
velocity  is  the  velocity  at  which  the  stress  wave 
propagates  from  one  point  to  another  in  the  medium.  For 
an  elastic  isotropic  medium,  this  velocity  remains 
constant  with  both  time  and  position  in  the  medium.  The 

velocity  is  the  velocity  of  the  mass  particles  at 
some  point  in  the  medium.  Since  the  wave  propagates 
through  the  .medium,  stressing  separate  sections  of  the 
medium  at  different  times,  the  particle  velocity  is  not 
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constant  throughout  the  medium.  It  should  also  be  noted 
that  the  particle  velocities  for  tension  and  compression 
waves  are  opposite  in  sense.  Mass  particles  subjected  to 
a  compressive  wave  travel  in  the  direction  of  wave 
propagation.  Mass  particles  subjected  to  a  tensile  wave, 
however,  travel  in  the  direction  opposite  that  of  wave 
propagation. 


Stress  Wave  Propagation  in  a  Uniform  Bar 
Equation  (2.3)  describes  the  motion  of  a  stress  wave 
propagating  through  an  isotropic  medium.  Theoretically, 
equation  (2.3)  can  be  used  to  describe  the  propagation  of 
stress  waves  in  any  isotropic  solid  by  solving  for  the 
appropriate  boundary  conditions.  However,  except  for  a 
few  simple  cases,  the  solution  becomes  very  complex. 
Therefore,  the  propagation  of  stress  waves  in  a  bar  of 
uniform  cross  section  is  investigated  in  order  to 
illustrate  the  phenomenon. 

A  differential  element  of  the  bar,  ojc,  is  shown  in 
Figure  2.3.  Assuming  plane  sections  remain  plane  and  the 
stress  is  uniformly  distributed  over  the  cross  section,  it 
follows  that 

a,  =  (2.11) 

*  ox 

Summing  forces  in  the  x  direction  to  obtain  the  residual 
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Equation  (2.12)  is  identical  with  equation  (2.7),  derived 
for  the  propagation  of  a  one-dimensional  plane  wave  in  a 
linear  isotropic  solid.  Therefore,  the  propagation  of  a 
plane  longitudinal  stress  wave  in  a  three-dimensional 
elastic  isotropic  medium  is  identical  to  the  propagatidn 
of  the  same  stress  wave  in  a  bar  of  uniform  cross  section. 
Equation  (2.12)  is  applicable  to  any  elastic  isotropic  bar 
of  uniform  cross  section.  The  relationship  between  c,/c 
and  Poisson's  ratio  v  is  illustrated  in  Figure  2.4.  It 
should  be  ncted  that  for  small  Poisson's  ratios,  the 
velocity  c,  in  equation  (2.4)  approaches  the  velocity  c  in 
equation  (2.13).  Therefore,  the  velocity  of  longitudinal 
stress  wave  propagation  in  a  linear  isotropic  solid  of 
small  Poisson  ratio  is  approximately  equal  to  the  velocity 
of  propagation  in  a  bar  of  uniform  cross  section.  In 
consideration  of  the  above  similarities,  stress  wave 
propagation  in  a  bar  of  uniform  cross  section  can  be  used 
to  model  the  plane  longitudinal  stress  wave  propagation  -in 
a  linear  isotropic  medium. 

The  general  solution  to  equation  (2.12)  is 

“(jf.O  =  /,{Jc  +  cO  +  fjix-ct)  (2.14) 

where  the  function  describes  a  backward  moving  wave, 
and  the  function  /j  describes  a  foward  moving  wave. 

Equating  impulse  to  the  change  in  momentum  yields 
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Figure  2.4.  Relationship  Between  Velocity  Ratio  and 
Poissons  Ratio 

<s  •  pcy  or  K  ■  — 1  (2.15) 

pc 

This  Is  consistent  with  equation  (2.9),  and  Indicates  that 
the  particle  velocity  of  the  bar  is  linearly  proportional 
to  the  stress  by  the  material  impedance  pc. 

For  a  foward  moving  stress  wave  propagating  in  the  x 
direction,  equation  (2.14)  reduces  to 

u(x,0  «/j(JC-cO.  (2.16) 

If  the  stress  wave  is  elastic  and  plane,  then  every 
segment  of  the  stress  wave  is  propagated  with  the  same 
velocity,  and  the  shape  of  the  stress  wave  does  not 
ol.dnge.  The  solution  given  by  equation  (2.16)  is  depicted 
graphically  in  Figure  2.5.  Figure  2.5(a)  illustrates  a 
foward  moving  stress  wave  ab  of  arbitrary  shape  at  ti.me  i. 
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propagating  in  the  x  direction  at  constant  velocity  c.  At 
time  +  (Figure  2.5(b))/  the  stress  wave  has  traveled 
a  distance  Bx*c5t  from  position  ab  to  position  a'b\  The 
shape  and  amplitude  of  the  stress  wave  remain  unchanged. 


Figure  2.5.  Foward  Moving  Stress  Wave  Propagating  With 
Velocity  c. 


Figure  2.5  represents  a  single  stress  wave 
propagating  through  a  uniform  bar  in  the  positive  x 
direction.  If  this  stress  wave  encounters’ a  second  stress 
wave  traveling  in  the  negative  x  direction,  the  two  waves 
will  interfere  with  one  another.  The  iwi.'ses  and 
particle  velocities  of  the  two  stress  waves  can  be 
superimposed  to  obtain  the  resultant  stresses  and  particle 
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velocities.  The  interference  of  two  compression  waves 
traveling  in  opposite  directions  is  shown  in  Figure  2.6. 
The  compression  waves  have  the  same  stress  level  ,  but 
have  particle  velocities  J',  in  opposite  directions,  as 
shown  in  Figure  2.6(a).  Figure  2.6(b)  reveals  the 
superposition  of  a  portion  of  the  two  waves.  The  stress 
level  is  additive  in  this  case,  and  the  compressive  stress 
intensity  doubles  to  2a^  The  particle  velocities,  equal  in 
magnitude  but  opposite  in  direction,  cancel  and  become 
zero. 

The  interference  of  a  compressive  wave  and  a  tensile 
wave  traveling  in  opposite  directions  is  shown  in  Figure 
2.7,  The  waves  are  of  equal  stress  level,  c, ,  but  of 
opposite  sign.  The  particle  velocities  are  identical  in 
both  magnitude  and  direction,  as  shown  in  Figure  2.7(a). 
The  superposition  of  the  two  stress  waves  is  depicted  in 
Figure  2.7(b).  In  this  case,  the  stresses  cancel  each 
other  out  and  the  particle  velocity  doubles  to  2v,. 


Boundary  Effects  on  Stress  Wave  Propagation 
A  stress  wave  will  continue  to  propagate  through  a 
linear,  isotropic  material  with  no  change  in  shape, 
intensity,  or  direction  until  it  encounters  a  boundary. 

In  general,  when  an  incident  longitudinal  stress  wave 
reaches  a  boundary,  it  is  reflected  as  two  separate  waves; 
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Figure  2.6.  Interference  and  Superposition  of  Two 
Compressive  Stress  Waves. 


1)  a  longitudinal  wave,  and  2)  a  distortional  wave. 
However,  when  the  incident  longitudinal  wave  is  normal  to 
the  boundary,  the  longitudinal  wave  is  reflected  as  a 
single  longitudinal  wave. 

Reflection  crom  a  Free  Surface 
A  longitudinal  wave  at  normal  incidence  to  a  free 
boundary  will  be  reflected  as  a  longitudinal  wave  with  a 


24-23 


(b) 

Figure  2.7.  Interference  and  Superposition  of  Compression 
and  Tension  Stress  Waves. 

change  in  phase.  The  three  boundary  conditions  governing 
the  reflection  from  a  free  surface  are:  1)  zero  stress, 

2)  nonzero  displacement,  and  3)  nonzero  particle  velocity 
at  the  free  surface.  An  incident  longitudinal  compression 
wave  will  be  reflected  at  full  stress  level  as  a 
longitudinal  tension  wave.  Superposition  of  the  incident 
and  reflected  stress  waves  yields  zero  stress  but  a 
doubling  of  the  particle  velocity  at  the  free  surface-. 
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Figure  2.8.  Reflection  of  a  Rectangular  Wave  at  a  Free 
Surface. 

The  reflection  of  a  longitudinal  stress  wave  at 
normal  incidence  to  a  free  surface  is  illustrated  in 
Figures  2.8  and  2.9  for  four  intervals  of  time.  A 
rectangular  compressive  stress  wave  reflecting  from  the 
free  end  v,r  i  bar  is  shown  in  Figure  2.8,  and  a  triangular 
compressive  stress  wave  reflecting  from  the  free  end  is 
shown  in  Figure  2.9,  The  compression  wave  is  reflected  at 
full  intensity  as  a  tensile  wave  when  it  reaches  the  free 
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surface.  It  is  noted  that,  upon  reflection  of  the 
rectangular  wave,  a  segment  of  the  bar  adjacent  to  the 
free  surface  remains  xinstressed.  However,  upon  reflection 
of  the  triangular  wave,  only  the  free  surface  remains 
unstressed.  In  either  case,  the  maximum  stress  due  to 
superposition  is  equal  to  the  maximum  incident  stress. 

a 


/ 


Figure  2.9.  Reflection  of  a  Triangular  Wave  at  a  Free 
Surface . 
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Reflection  From  a  Fixed  Boundary 

A  longitudinal  wave  at  normal  incidence  to  a  fixed 
boundary  will  be  reflected  as  a  longitudinal  wave  with  no 
change  in  phase.  The  three  boundary  conditions  governing 
the  reflection  from  a  fixed  surface  are:  1)  zero 
displacement,  2)  zero  particle  velocity,  and  3)  nonzero 
stress  at  the  fixed  surface.  An  incident  longitudinal 
compressior.  wave  will  be  reflected  at  full  stress  level  as 
a  compression  wave.  Superposition  of  the  incident  and 
reflected  stress  waves  yields  zero  particle  velocity,  but 
a  doubling  of  the  stress  intensity  at  the  fixed  surface. 

The  reflection  of  a  longitudinal  stress  wave  at 
normal  incidence  to  a  fixed  surface  is  illustrated  in 
Figures  2.10  and  2.11  for  four  intervals  of  time.  A 
rectangular  compressive  stress  wav«  reflecting  from  the 
fixed  end  of  a  bar  is  shown  in  Figure  2.10,  and  a 
triangular  compressive  stress  wave  reflecting  from  the 
fixed  end  is  shown  in  Figure  2.11.  The  compression  wave 
is  reflected  at  full  intensity  without  a  change  of  phase 
when  it  reaches  the  fixed  surface.  It  is  noted  that,  in 
the  case  of  an  incident  triangular  wave,  a  rectangular 
wave  of  varying  length  and  intensity  is  formed  as  the 
incident  and  reflected  waves  are  superimposed.  In  either 
case,  however,  the  maximum  stress  due  to  superposition  is 
twice  that  of  the  maximum  stress  of  the  incident  wave. 
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Figure  2.10.  Reflection  of  a  Rectangular  Wave  at  a  Fixed 
Surface. 


Stress  Wave  Propagation  in  Layered  Media 
The  relative  strength  and  density  of  two  dissimiliar 
materials  forming  a  boundary  affects  stress  wave 
propagation  in  the  media.  In  general,  when  an  incident 
longitudinal  stress  wave  reaches  a  boundary  between  two 
dissimiliar  materials,  the  wave  is  reflected  and  refracted 
(transmitted)  as  four  separate  waves:  1)  a  reflected 
longitudinal  wave,  2)  a  reflected  d: "tortional  wave,  3)  a 
transmitted  longitudinal  wave,  and  4)  a  transmitted 
distortional  wave.  When  the  incident  longitudinal  wave  is 
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Figure  2.11.  Reflection  of  a  Triangular  Wave  at  a  Fixed 
Surface. 

normal  to  the  boundary,  however,  the  incident  longitudinal 
wave  is  reflected  and  transmitted  as  a  longitudinal  wave. 

The  boundary  conditions  governing  the  reflection  and 
refraction  at  the  boundary  between  two  dissimiliar 
materials  are  defined  as  the  equality  of  stress  and 
particle  velocities  in  each  material  at  the  boundary, 
(assuming  the  materials  remain  in  constant  contact). 

These  conditions  are  described  by  the  expressions 

aj(x,t)  +  a^{x,i)  =  (2.17) 

and 
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(2.18) 


where  <s,,  a,,  Qj.  are  the  incident,  reflected,  and 
transmitted  stresses,  respectively;  V,,  Vj.  are  the 
respective  incident,  reflected,  and  transmitted  particle 
velocities,  where 


P. 


-11. 

P,c, 


p.c 


2 '•2 


as  given  by  equation  (2.15).  The  subscripts  1  and  2  refer 
to  the  first  and  second  materials,  respectively. 
Substituting  equation  (2.15)  for  V  in  equation  (2.18) 
yields 


_ 


Pi  <=i 


P2=2 


(2.19) 


Solving  equations  (2.17)  and 

cr  *  —  ,  ■■  ■  <j 
p,c, +  P2C2 

^  P,c,  +  PjCj  / 


(2.19)  simultaneously  yields 

(2.20) 

(2.21) 


Rearranging  equations  (2.20)  and  (2.21)  yields 

,  2p^c^ 

”“piCi  +  PjC, 


a,  P,c,  +  p,Cj 


(2.22) 

(2.23) 


and  dividing  equation  (2.22)  by  equation  (2.23)  results  in 
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(2.24) 


Of  ^  ^  ^2  ^2 

PjCj-p.c, 

Equations  (2.20)  through  (2.24)  also  describe  the 
fixed  and  free  boundary  conditions  previously 
investigated.  WhenpjCj»0,  i.e./  a  free  boundary 

(7 

condition,  then  C;.  ■  0  and— «  -1,  denoting  that  an  incident 

compression  wave  will  be  reflected  at  full  stress  level  as 
a  tension  wave.  When  p,c,»oo,  describing  a  fixed  boundary 

Cr 

condition,  —  *  2  and  —  «  I,  indicating  that  a  compression 
a,  a, 

wave  reflects  at  full  intensity  without  a  change  in  phase, 
and  the  boundary  will  experience  a  stress  equal  to  twice 
that  of  the  incident  stress. 

The  transmission  ratio  reflection  ratio 

are  plotted  in  Figure  2.12  for  a  stress  wave  at 

normal  incidence.  When  p,  c,  »  Pj  Cj ,  as  in  a  continuous 
media,  there  will  be  no  reflection  and  the  stress  wave 
will  be  transmitted  at  full  intensity.  When  p,  c,  <  p^c,  ,  the 
reflection  ratio  is  positive,  and  an  incident  stress  wave 
will  be  reflected  with  no  change  in  phase.  When  p,  c,  >  PjC,  , 
the  reflection  ratio  is  negative,  and  an  incident  stress 
wave  will  be  reflected  with  a  change  in  phase.  It  should 
be  noted,  however,  that  if  the  boundary  between  the 
materials  cannot  support  tension,  then  an  incident  tensile 
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wave  will  be  reflected  as  if  at  a  free  boundary.  The 
lower  and  upper  limits  of  the  transmission  ratio  are  0  and 
2,  respectively,  representing  the  extreme  cases  of  free 
and  fixed  boundaries.  The  upper  and  lower  limits  of  the 
reflection  ratio  are  -1  and  1,  representing  the  limiting 
cases  of  free  and  fixed  boundaries,  respectively. 


Figure  2.12.  Transmission  and  Reflection  Ratios  at  Normal 
Incidence  vs.  the  Ratio  Between  Material 
Impedances . 


The  reflection  and  transmission  of  a  triangular 
stress  wave  is  illustrated  in  Figure  2.13  for  p,  c,  <  p, c,.  The 
stress  wave  is  reflected  at  less  than  full  intensity  with 
no  change  in  phase,  and  is  transmitted  at  an  increased 
stress  level.  The  length  of  the  stress  wave  increases  in 
the  second  material  due  to  the  increased  wave  velocity  c,. 
of  the  second  .material. 
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Figure  2.13.  Transmission  and  reflection  of  a  stress  wave 
for  P,  c,  <  Pj  Cj  . 

The  reflection  and  transmission  of  a  triangular 
stress  wave  is  illustrated  in  Figure  2.14  for  p,  c,  >  Pj  c,.  The 
stress  wave  is  reflected  at  less  than  full  intensity  with 
a  change  in  phase,  and  is  transmitted  at  a  reduced  stress 
level.  The  length  of  the  stress  wave  decreases  in  the 
second  material  due  to  the  reduced  wave  velocity  c,  of  the 
second  material. 


Plastic  Stress  Wave  Propagaticu 


The  propagation  of  stress  waves  in  a  nonlinear 


material  is  a  very  complex  phenomenon.  This  study  will 
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Figur®  2.14.  Transmission  and  reflection  of  a  stress  wave 
for  p,c,  >PjCj. 

address  the  basic  concepts  of  the  subject.  For  a  more 
detailed  investigation,  the  reader  is  referred  to  (10), 
(12),  (17),  and  (27). 

The  velocity  of  stress  wave  propagation  in  an  elastic 
medium  has  been  previously  defined  by  equation  (2.13)  as 


The  velocity  of  a  plastic  stress  wave  is  a  function  of  the 
slope  of  the  stress-strain  curve  and  is  defined  (6)  by 


I 


If  the  material  is  stressed  in  the  elastic  region  only, 
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equation  (2.25)  simplifies  to.  equation  (2.13)  by 

substituting  E  for  in  equation  (2.25).  A  stress- 

dc 

strain  curve  for  a  typical  elasto-plastic  material  is 
shown  in  Figure  2.15  along  with  the  corresponding  stress 
wave  velocity  obtained  from  equation  (2.25).  It  is 
indicated  in  Figure  2.15  that,  for  a  typical  elasto- 
plastic  material,  the  plastic  stress  wave  travels  at  a 
much  lower  velocity  than  the  elastic  portion  of  the  stress 
wave.  The  plastic  wave  front  will  therefore  lag  behind 
the  elastic  wave  front. 


Figure  2.15.  Stress-Strain  Curve  for  a  Typical  Elasto- 
Plastic  Material  and  Corresponding  Wave 
Velocity. 


A  representation  of  a  typical  stress-strain  curve  for 
an  elasto-plastic  material  in  compression  is  shown  in 
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Figure  2.16.  If  a  high  intensity  shock  wave,  such  as  that 
shown  in  Figure  2.17,  is  applied  to  a  material  having  a 
constitutive  relation  similiar  to  that  shown  in  Figure 
2.16,  the  stress  wave  will  propagate  through  the  medium  as 
illustrated  in  Figure  2.18  (10)-.  The  plastic  portion  of 
the  wave  propagates  at  a  velocity  less  than  that  of  the 
elastic  portion  and  lags  behind  the  elastic  wave  front. 

The  wave  can  be  separated  into  an  unstable  shock  wave,  a 
plastic  wave,  and  an  elastic  wave.  After  the  wave  has 
traveled  a  sufficient  distance,  the  unstable  shock  wave  is 
reduced  to  a  plastic  wave.  The  plastic  wave  continues  to 
lag  behind  the  elastic  wave  front  until  it  has  been 
completely  reduced  to  an  elastic  wave. 


Figure  2.16.  Stress-Strain  Curve  for  an  Elasto-Plastic 
Material  in  Compression. 


In  general,  any  high  intensity  elasto-plastic  stress 
wave  will  continually  change  shape  until  it  has  been 
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Figure  2.17.  Applied  High  Intensity  Shock  Wave. 


Figure  2.18.  Propagation  of  a  High  Intensity  Stress  Wave 
in  an  Elasto-Plastic  Medium. 

completely  reduced  to  an  elastic  wave.  The  stress  wave  is 
not  only  reduced  in  intensity,  but  is  also  lengthened,  due 
to  the  reduced  velocity  of  the  plastic  .portion  of  the 
wave . 
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Spalling* is  defined  as  fracturing  caused  when  a  high 
intensity  transient  stress  wave  reflects  from  a  free 
surface  (19).  Spalling  is  a  result  of  the  interference  of 
an  incident  compression  wave  with  its  reflected  tensile 
wave  at  a  free  surface.  The  incident  triangular 
compression  wave  shown  in  Figure  2.19(a)  is  reflected  as  a 
tensile  wave  from  the  free  surface  as  shown  in  Figure 
2.19(b).  The  superpostion  of  the  incident  and  reflected 
stress  waves  results  in  a  net  tensile  stress  near  the  free 
surface.  The  tensile  stress  increases  from  zero  at  the 
free  surface  to  a  maximum  value  of  c,  some  distance  from 
the  free  surface.  The  distance  from  the  free  surface  at 
which  the  net  tensile  stress  reaches  the  maximum  value  a, 
is  a  function  of  the  wavelength  of  the  incident  stress 
wave.  The  maximum  possible  value  of  tensile  stress  due  to 
the  superposition  of  the  incident  and  reflected  triangular 
waves  is  shown  in  Figure  2.20.  The  maximum  tensile  stress 
9,  is  reached  at  a  distance  of  one  half  the  wave  length 
from  the  free  surface. 

For  a  rectangular  incident  compression  wave,  as  shown 
in  Figure  2.21(a),  the  reflection  from  a  free  surface  is 
illustrated  in  Figure  2.21(b).  The  superposition  of  the 
incident  and  reflected  stress  waves  results  in  an 
unstressed  section  of  the  medium  adjacent  to  the  free 
surface.  However,  a  tensile  stress  of  intensity 


instantaneously  develops  at  some  distance  from  the  free 
surface.  The  distance  from  the  free  surface  at  which  the 
net  tensile  stress  forms  is  a  function  of  the  wavelength 
of  the  incident  stress  wave.  The  maximum  possible  value 
of  tensile  stress  due  to  the  superposition  of  the  incident 
and  reflected  rectangular  waves  is  shown  in  Figure  2.22. 
The  maximum  tensile  stress  a,  is  attained  at  a  distance  of 
one  half  the  wave  length  from  the  free  surface.  For  both 
the  rectangular  and  triangular  wave  shapes,  the  stress  at 
the  free  surface  is  always  zero. 


Figure  2.19.  Superposition  of  Stresses  During  the 

Reflection  of  of  a  Triangular  Wave  from  a 
Free  Surface. 


For  a  medium  capable  of  supporting  the  tensile  stress 
inflicted  by  wave  reflection,  the  reflected  tensile  wave 
will  continue  to  propagate  with  no  change  in  shape  or 
intensity,  and  no  spalling  will  occur.  If,  however,  the 
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Figure  2.20.  Maximum  Attainable  Tensile  Stress  During  the 
Reflection  of  a  Triangular  Wave  from  a  Free 
Surface. 


(a) 


compression 

tension 


(b) 


Figure  2,21.  Superposition  of  Stresses  During  the 

Reflection  of  of  a  Rectangular  Wave  from  a 
Free  Surface. 


medium  is  unable  to  support  the  tensile  stress,  spalling 
will  occur  and  propagation  of  the  reflected  stress  wave 
will  be  altered.  The  occurence  of  sralling  and  the 
location  of  the  fractures  are  dependent  upon  three 
factors:  1)  the  resistance  of  the  material  to  fracture; 
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Flgurt  2.22.  Maximum  Attainable  Tensile  Stress  During  the 
Reflection  of  a  Rectangular  Wave  from  a  Free 
Surface. 

2)  the  Intensity  of  the  tensile  stress;  and  3)  the  shape 
of  the  stress  wave  (19). 

The  stress  level  at  which  fracturing  occurs  in  a 
material  is  called  the  critical  normal  fracture  strength. 
In  general,  the  critical  normal  fracture  strength  is 
considerably  higher  than  the  static  tensile  strength  of 
the  material.  However,  it  is  difficult  to  experimentally 
determine  the  exact  state  of  stress  at  the  instant  of 
fracture  because  of  the  instantaneous  stress  reversal  at 
the  point  of  fracture  (19). 

Spalling  occurs  when  the  critical  normal  fracture 
strength  of  the  material  is  reached.  Multiple  spalling  is 
the  occurrence  of  spall  in  layers  parallel  to  the  free 
surface.  Multiple  spalling  occurs  when  the  stress  level 
is  greater  than  twice  the  critical  normal  fracture 
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strength.  The  occurence  of  a  multiple  spall  caused  by  an 
Incident  compression  wave  is  shown  in  Figure  2.23.  The 
tensile  stress  resulting  from  the  superposition  of  the 
incident  and  reflected  waves  increases  until  the  fracture 
strength  is  attained.  Once  the  f];acture  strength  is 
exceeded,  the  material  spalls,  creating  a  second  free 
surface.  If  the  remaining  portion  of  the  incident  wave  is 
greater,  than  the  fracture  strength  of  the  material,  then  a 
second  spall  will  occur  due  to  the  reflection  from  the 
newly  created  free  surface.  This  process  conti'^vea  until 
the  stress  wave  no  longer  exceeds  the  critical  normal 
fracture  strength  of  the  material. 

The  thickness  of  the  spalled  layer  is  affected  by  the 
stress  wave  profile,  as  indicated  in  Figure  2.23.  A 
rectangular  stress  wave  causes  thicker  spall  than  a 
triangular  stress  wave.  This  is  a  result  of  the  shape  of 
the  maximum  net  tensile  stress  distribution  formed  by  the 
superposition  of  the  incident  and  reflected  stress  waves 
(refer  to  Figures  2.20  and  2.22).  The  net  tensile  stress 
resulting  from  a  triangular  stress  wave  builds 
immediately,  thus  restricting  spall  close  to  the  vicinity 
of  the  free  surface.  However,  the  net  tensile  stress 
resulting  from  a  rectangular  stress  wave  remains  zero  for 
a  distance  equal  to  one  half  the  length  of  the  incident 
wave  before  reaching  its  maximum  value,  thus  causing 
thicker  spall. 
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Figure  2.23.  Multiple  Spalling  Caused  by  a  Single 
Incident  Compression  Wave. 


The  particle  velocity  is  linearly  related  to  stress 
by  equation  (2.15).  Also,  the  superposition  of  particle 
velocities  for  compression  and  tension  waves  traveling  in 
opposite  directions  is  additive.  Therefore,  the  greater 
the  stress  level  of  the  incident  wave,  the  greater  the 
velocity  at  which  the  spall  is  propagated  away  from  the 
free  surface.  In  multiple  spalls,  caused  by  a  stress  wave 
such  as  that  shown  in  Figure  2.23,  the  first  spall  will 
have  the  greatest  velocity,  and  each  subsequent  spall  will 
have  a  decreased  velocity  due  to  the  reduced  intensity  of 
the  incident  stress  wave. 


III.  LAYERED  MEDIA  FOR  PROTECTIVE  MILITARY  SHELTERS 


Protective  military  shelters  are  designed  to  house 
personnel,  vital  functions,  and  equipment  of  value.  As  a 
result,  the  survivability  of  the  structure  is  critical  and 
takes  precedence  over  its  appearance.  The  structures  are 
invariably  massive,  constructed  primarily  of  concrete  and 
soil.  Damage  to  protective  shelters  from  conventional 
(non-nuclear)  weapons  occurs  as  a  result  of  one  or  more  of 
the  following  effects;  penetration,  fragmentation,  ground 
shock,  and  blast  (9).  Penetration  is  a  consideration  only 
in  the  case  of  a  direct  hit.  Ground  shock  is  a 
consideration  only  in  the  case  of  a  near  miss. 
Fragmentation  and  blast  are  considerations  in  both  cases. 
This  study  considers  only  the  effects  of  the  blast. 


Typical  Military  Protective  Shelters 
Typical  aboveground  and  underground  protective 
shelters,  designed  to  survive  a  near  miss  of  a 
conventional  weapon,  are  shown  in  Figures  3.1  and  3.2, 
respectively.  The  protective  shelter  is  typically 
consx-ructed  of  thick,  monolithic  concrete  walls  and  roof. 
This  practice  is  considered  necessaiy  in  order  to  resist 
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Monolithic 


Figure  3.2*.  Typical  Underground  Protective  Shelter. 

the  blast  pressure  and  prevent  interior  spalling  of  the 
concrete  walls. 

A  cross  section  of  a  typical  protectiVv-  shelter  wall 
is  shown  in  Figure  3.3.  When  the  shelter  wall  is  subject 
to  a  blast,  the  compressive  longitudinal  blast  wave 
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contacts  the  exterior  surface  of  the  wall  and  propagates 
through  the  monolithic  wall  with  little  or  no  change  in 
shape  or  intensity.  When  the  compressive  stress  wave 
reaches  the  interior  face  of  the  wall,  it  is  reflected  at 
full  stress  level  as  a  tension  wave.  The  concrete  will 
spall  at  the  interior  face  of  the  wall  if  the  intensity  of 
the  reflected  tension  wave  is  greater  than  the  critical 
normal  fracture  strength  of  the  concrete.  When  the  stress 
wave  is  intense,  the  spalled  concrete  fragments  are  hurled 
away  from  the  wal’  at  high  velocities,  posing  a 
potentially  dangerous  threat  to  personnel  and  equipment 
within  the  shelter. 


Figure  3.3.  Blast  Wave  Propagation  in  Monolithic  Wall  of 
a  Typical  Protective  Shelter. 
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Current  procedures  to  reduce  or  eliminate  spall  in 
protective  shelters  include  the  construction  of  earth 
berms  against  the  exterior  face  of  the  shelter  walls 
(Figure  3.4(a)),  the  installation  of  steel  spall  plates  on 
the  interior  face  of  the  shelter  walls  (Figure  3.4(b)), 
and  the  use  of  thicker  concrete  walls  (Figure  3.4(c))  (7). 
Of  these  methods,  earth  berms  have  proven  to  be  the  most 
successful,  while  the  use  of  thicker  walls  has  been  the 
least  successful. 


(a)  BERMS  (b)  SPALL  PLATES  (c)  THICKER  WALLS 


Figure  3.4.  Current  Procedures  to  Reduce  Spall  in 
Protective  Shelters. 

Layered  Media  Applied  to  Protective 
Military  Shelters 

It  was  shown  in  Chapter  II  that  the  propagation  of 
stress  waves  in  solid  media  is  altered  at  the  boundaries 
between  dissimiliar  materials.  Therefore,  it  is  logical 
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to  assxime  that  layered  media  can  be  utilized  to  alter  the 
stress  wave  propagation  in  the  walls  of  a  protective 
structure.  The  use  of  layered  media  to  mitigate  the 
incident  blast  wave,  and  therefore  eliminate  spall,  in 
protective  shelters  is  investigated  in  this  study.  The 
layered  systems  considered  consist  of  an  exterior  layer  of 
concrete  or  sand,  an  absorption  layer,  and  an  interior 
layer  of  concrete. 

The  transmission  ratio  of  a  three-layered  syste.m  is 
determined  by  two  consecutive  applications  of  equation 

(2.22) .  Equation  (2.22)  is  first  applied  to  determine  the 
stress  wave  transmitted  from  the  first  layer  to  the  second 
layer.  The  stress  transmitted  to  the  second  layer  is  the 
incident  stress  for  the  third  layer.  Equation  (2.22)  is 
applied  a  second  time  to  obtain  the  stress  transmitted  to 
the  third  layer.  Two  consecutive  applications  of  equation 

(2.22)  yields 


1  ^ 

11 

2  p,  c. 

2  pj  c. 

_P,c,  +  p,c,_ 

p,  C,  +  pj  Cj 

where  the  subscripts  1,  2,  and  3  refer  to  the  first, 
second,  and  third  layers,  respectively.  For  a  layered 
system  in  which  the  first  and  third  layers  consist  of  the 
same  material  (concrete)  and  the  second  layer  is  an 
absorption  layer,  equation  (3.1)  becomes 
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where  the  si:ibscripts  c  and  a  represent  the  concrete  and 
absorption  layers,  respectively.  A  BASIC  computer  program 
was  developed  to  calculate  the  transmission  ratio  for 
layered  systems.  The  program  listing  and  sample  output  is 
presented  in  Appendix  A. 

Equation  (3.2)  is  depicted  graphically  in  Figure  3.5. 

The  transmission  ratio  for  the  system,  is  plotted 

^1 

against  the  ratio  of  impedances  of  the  concrete  and 


absorption  layers. 


P  c 

-i-i.  If  the  absorption  layer  has  the 

P.c. 


same  impedance  as  the  concrete  layer,  then  the 
transmission  ratio  is  equal  to  one,  and  no  reflection 
occurs  at  the  boundary  between  layers.  As  the  impedance 
ratio  diverges  from  one  in  either  direction,  the 
transmission  ratio  of  the  system  decreases.  An  impedance 
ratio  less  than  one  is  necessary  if  the  blast  wave  is  to 
be  significantly  reduced,  as  is  illustrated  in  Figure  3.5. 

The  blast  wave  propagation  through  a  layered  system 
with  an  impedance  ratio  less  than  one  is  illustrated  in 
Figure  3.6.  The  incident  compressive  longitudinal  blast 
wave  propagates  through  the  exterior  layer  of  concrete 
until  the  first  material  boundary  is  encountered.  Upon 
reaching  the  boundary  between  the  exterior  concrete  and 
the  absorp  .ayer,  a  portion  of  the  blast  wave  is 

transmitted  to  the  absorption  layer.  The  remainder  of  the 
blast  wave  is  reflected  back  through  the  exterior  layer  of 
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Figure  3.5.  Transmission  Ratio  vs.  the  Impedance  Ratio 
for  the  Three  Layered  System. 


concrete  as  a  tensile  wave.  The  transmitted  compression 
wave  propagates  through  the  absorption  layer  until  it 
encounters  the  interior  layer  of  concrete.  At  the 
boundary  between  the  absorption  layer  and  the  interior 
layer  of  concrete,  the  compression  wave  is  reflected  back 
through  the  absorption  layer  as  a  compression  wave  of 
diminished  intensity.  Because  the  compression  wave;  is 
reflected  with  no  change  in  phase,  the  stress  wave 
transmitted  to  the  interior  layer  of  concrete  is  greater 
than  the  reflected  wave  in  the  absorption  layer.  H'^wo^er, 
the  stress  wave  transmitted  to  the  interior  layer  of 
concrete  is  significantly  less  than  the  intial  incident 
blast  wave  which  impinged  upon  the  e.xterior  layer  of  the 
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system.  If  the  stress  wave  transmitted  to  the  interior 
layer  of  concrete  is  less  than  the  critical  normal 
fracture  strength  of  the  concrete,  then  no  spalling 
occurs. 


Figure  3.6.  Blast  Wave  Propagation  in  a  Layered  System 

Consisting  of  an  Exterior  Concrete  Layer,  an 
Absorption  Layer,  and  an  Interior  Concrete 
Layer  for  (p,  c,)  <  (p^cj. 


Previous  research  (14,23,24)  has  shown  polystyrene  to 
be  extremely  effective  as  an  absorption  layer.  It 
exhibits  very  low  impedance  characteristics  and  is 
available  in  large  sheets  which  afford  expedient 
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construction.  For  these  reasons,  polystyrene  was  selected 
for  the  absorption  layer  in  this  study. 

The  layered  systems  investigated  include  a  concrete- 
polystyrene-concrete  (CPC)  system,  a  sand-polystyrene- 
concrete  (SPC)  system,  and  a  sand-polystyrene-concrete- 
polystyrene-concrete  (SPCPC)  system.  The  static,  linear 
properties  of  the  materials  are  given  in  Table  3.1.  The 
elastic  wave  speeds  and  material  impedances  for  these  same 
materials  are  listed  in  Table  3.2. 


Table  3.1.  Static  Linear  Material  Properties. 


Strength 

£  ipsi) 

Mass 

Density 

p  (10-^^ 

tn 

Poisson' s 
Ratio 

V 

4000  psi 
Concrete 

3605000 

217.164 

.18 

110  pcf 

Sand 

80000 

164.745 

.25 

5.744  pcf 
Polystyrene 

2100 

1 

8.604 

_ d!_ 

The  intensity  of  the  blast  wave  is  capable  of 
stressing  the  shelter  into  the  inelastic  range.  This 
nonlinearity  is  addressed  to  some  extent.  It  is  assumed 
that  the  concrete  behaves  linearly  at  all  times.  Although 
this  IS  not  the  case  in  the  actual  structure,  the 
assumption  of  linearity  does  not  significantly  affect  the 
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Table  3.2.  Elastic  Wave  Speeds  and  Impedances. 


Elastic 

Wave  Speed 

c  inisec 

Impedance 

lb-see 
pc  -T7" 
m 

4000  psi 
Concrete 

128842 . 4 

28.0 

110  pcf 

Sand 

22036.3 

3.63 

5.744  pcf 
Polystyrene 

15622.8 

.134 

overall  behavior  of  the  layered  system.  Both  linear  and 
nonlinear  material  models  are  employed  for  the  sand  and 
polystyrene.  A  linear  elastic-perfectly  plastic 
constitutive  law,  similiar  to  that  shown  in  Figure  3.7,  is 
used  to  model  the  nonlinear  behavior  of  the  sand  and 
polystyrene.  The  parameters  for  the  nonlinear  sand  and 
polystyrene  models  are  presented  in  Table  3.3. 

The  material  properties  given  in  Tables  3.1  and  3.3 
represent  the  material  behavior  under  static  loading 
conditions.  A  protective  shelter  subjected  to  a  blast 
loading  experiences  a  high  intensity,  short  duration  load 
at  a  very  rapid  rate  of  application.  Materials  are  often 
able  to  tolerate  increased  stress  levels  when  loaded  at  a 
high  rate.  Laboratory  tests  performed  on  concrete 
specimens  under  high  rates  of  loading  have  demonstrated 
that  the  dynamic  secant  modulus,  ultimate  strength,  and 
tensile  strength  greatly  exceed  the  corresponding 


24-53 


Figure  3.7.  Linear  Elastic-Perfectly  Plastic  Constitutive 
Law  Used  to  Model  Nonlinear  Behavior. 


Table  3.3.  Static  Nonlinear  Material  Properties. 


Strength 

E  (jjst) 

Yield 

Stress 

(pit) 

Mass 

Density 

m 

Poisson*  s 
Ratio 

V 

Unconfined 

Sand 

80000 

20 

164.745 

.25 

5.744  pcf 
Polystyrene 

2100 

55 

8.604 

.  15 

strengths  determined  from  static  tests  (9,15,16,18). 
Under  dynamic  loading,  the  compressive  strength  of 
concrete  can  increase  by  as  much  as  a  factor  of  1.8,  and 
the  tensile  strength  can  increase  by  as  much  as  a  factor 
of  2.2. 

Tests  on  polystyrene  at  high  rates  of  Loading  have 
indicated  that  the  dynamic  ultimate  strength  of 
polystyrene  is  much  greater  than  the  static  strength. 
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Polystyrene  subjected  to  a  dynamic  loading  can  withstand  a 
stress  as  much  as  1.6  times  that  of  the  static  ultimate 
strength  (11).  Soils  also  exhibit  a  tendancy  to  withstand 
higher  loads  under  dynamic  loading  conditions.  The 
increase  in  strength  is  a  function  of  moisture  content, 
porosity,  and  pre-existing  static  stress.  However,  due  to 
the  uncertainties  involved  in  accurately  assessing  the 
high  strain  rate  effects  on  material  strength,  the  use  of 
the  dynamic  material  properties  in  this  study  is  limited. 

Concrete-Polystyrene-Concrete  System 

The  concrete-polystyrene-concrete  (CPC)  system  is 
intended  to  be  used  for  an  aboveground  protective  shelter. 
A  typical  cross  section  of  the  CPC  system  is  shown  in 
Figure  3.8.  The  exterior  concrete  layer  is  intended  to 
partially  withstand  blast  effects  and  to  stop  any  incident 
fragments.  Thn  polystyrene  absorption  layer  causes  a 
significant  portion  of  the  blast  wave  to  be  reflected  back 
through  the  exterior  concrete  layer.  The  polystyrene 
layer  also  traps  any  concrete  fragments  spalled  from  the 
exterior  layer.  The  interior  layer  of  concrete  resists 
that  portion  of  the  stress  wave  which  propagates  through 
the  polystyrene  layer. 

The  aboveground  CPC  system  may  be  constructed  by 
forming  the  concrete  walls  against  sheets  of  polystyrene. 

A  possible  construction  procedure  consists  of  casting  the 
interior  concrete  wall  against  polystyrene  sheets  backed 
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EXTERIOR 

CONCRETE 


Figure  3.8.  Cross  Section  of  the  Concrete-Polystyrene- 

Concrete  (CPC)  Layered  System. 

by  plywood  or  metal  forma.  After  the  forms  are  removed, 
the  outside  surface  of  the  exterior  wall  is  formed,  and 
the  exterior  layer  of  concrete  is  placed. 

Sand-Polystyrene-Concrete  System 
The  sand-polystyrene-concrete  (SPC)  system  is  intended  to 
be  used  for  both  aboveground  and  underground  protective 
shelters.  A  typical  cross  section  of  the  SPC  system  is 
shown  in  Figure  j.9.  The  layer  of  sand  is  intended  to 
partially  withstand  the  blast  -effects  and  to  stop  any 
incident  fragments.  The  polystyrene  absorption  layer 
refler^*-  portion  of  the  blast  wave  back  through  the  sand 
and  reduces  the  intensicy  of  the  blast  which  wave  reaches 
the  interior  layer  of  concrete.  The  interior  layer  of 
concrete  resists  the  remaining  portion  of  the  blast 
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Figure  3.9.  Cross  Section  of  the  Sand-Polystyrene- 
Concrete  (SPC)  Layered  System. 

Construction  of  the  aboveground  SPC  system  consists 
of  placing  sheets  of  polystyrene  against  the  walls  of  an 
existing  monolithic  concrete  shelter.  A  sand  berm  is  then 
constructed  against  the  polystyrene  sheets  to  complete  the 
layered  system.  A  typical  aboveground  SPC  system  is 
depicted  in  Figure  3.10. 

Construction  of  the  underground  SPC  system  consists 
of  placing  sheets  cf  polystyrene  against  the  walls  and 
roof  of  a  monolithic  concrete  shelter.  The  layered  system 
is  completed  by  backfilling  against  the  polystyrene  sheets 
with  sand.  A  typical  underground  SPC  system  is 
illustrated  in  Figure  3.11. 
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^ORtGJNAL  GROUND 

Figure  3.10.  Typical  Aboveground  SPC  Layered  System. 


Figure  3.11.  Typical  Underground  SPC  Layered  Syste.m. 


Sand-PoLy Styrene-Concrete- Polystyrene- 
Concrete  System 

The  sand-polystyrene-concrete-polystyrene-concrete  ( 3PCPC ) 
layered  system  is  intended  for  use  as  an  aboveground 
shelter  when  very  high  intensity  blasts  are  expected. 
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A  typical  cross  section  of  the  SPCPC  system  is  shown  in 
Figure  3.12.  The  sand  layer  partially  withstands  the 
effects  of  the  blast  wave  and  stops  any  incident 
fragments.  The  first  polystyrene  layer  reflects  a  portion 
of  the  blast  wave  back  through  the  sand,  and  reduces  the 
blast  wave  which  reaches  the  exterior  layer  of  concrete. 
The  remaining  concrete-polystyrene-concrete  layers  resist 
the  transmitted  blast  wave  as  described  in  the  section  on 
the  CPC  layered  system. 


Figure  3.12.  Cross  Section  of  the  SPCPC  Layered  System. 

Construction  of  the  aboveground  SPCPC  system  consists 
of  adding  the  exterior  polystyrene  and  soil  layers  to  a 
CPC  system.  After  the  CPC  system  has  been  constructed  and 
the  exterior  forms  removed,  sheets  of  polystyrene  are 
placed  against  the  exterior  concrete  walls.  A  soil  berm 
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is  then  constructed  against  the  polystyrene  sheets, 
completing  the  layered  system. 

Protective  Shelter  Loading  Condition 
An  explosion  is  a  violent  release  of  energy  from 
within  a  confined  space.  The  explosion  causes  a  sudden 
pressure  increase  in  the  surrounding  medium,  resulting  in 
the  formation  of  a  blast  wave.  The  blast  wave  travels 
away  from  the  point  of  detonation  at  a  velocity  in  excess 
of  t'*y.  sonic  velocity  of  the  surrounding  medium.  The 
velocity  diminishes  with  distance,  but  remains  in  excess 
of  the  sonic  velocity.  The  variation  in  velocity  leads  to 
the  rapid  formation  of  a  shock  front,  as  illustrated  in 
Figure  3.13  (9).  The  shock  front  decays  with  distance 
from  the  point  of  detonation,  as  illustrated  in  Figure 
3.14,  and  eventually  drops  below  atmospheric  pressure  and 
becomes  negative  (9). 

The  free-field  time  history  of  the  blast  wave  at  a 
specified  distance  from  the  point  of  detonation  is 
depicted  in  Figure  3. IS.  The  shock  front  arrives  at  time 
with  a  peak  pressure  of  The  pressure  then  decays  to 
the  ambient  pressure  at  time  r, ,  thus  ending  the  positive 
phase  of  the  blast.  The  pressure  then  becomes  negative, 
beginning  the  negative  phase  of  the  blast.  The  negative 
phase  lasts  longer  than  the  positive  phase,  but  its 


24-50 


Figure  3.13.  Formation  of  a  shock  front 


Figure  3.14.  Variation  of  a  shock  wave  with  distance  at 
selected  times  (13). 

intensity  is  small  in  comparison  with  the.  positive  phase 
and  is  usually  neglected. 

When  the  blast  wave  comes  in  contact  with  a  rigid 
surface,  such  as  the  wall  of  a  protective  shelter,  a 
reflected  pressure  in  excess  of  the  incident  pressure  is 
immediately  formed.  A  typical  reflected  pressure- time 
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Figure  3.15.  Fre«-field  pressure-time  history  of  a  given 
distance  from  the  point  of  detonation  (9). 

history  is  shown  in  Figure  3.16.  The  incident  pressure 
causes  the  reflected  pressure  /*,  to  develop.  The 
magnitude  of  the  reflected  pressure  is  a  function  of  the 
incident  pressure  and  the  angle  formed  between  the  wave 
front  and  the  rigid  surface. 

The  magnitude  and  distribution  of  the  blast  pressure 
on  a  structure  are  dependent  on  the  type  of  explosive,  the 
weight  of  the  explosive,  the  location  of  the  explosion 
with  respect  to  the  structure,  the  interaction  of  the 
shock  front  with  the  ground  and  surrounding  obstructions, 
and  the  interaction  of  the  shock  front  with  the  structure 
itself  (9).  A  typical  reflected  blast  pressure-time 
history  measured  in  half-scale  field  tests  (7)  is  shown  in 
Figure  3.17.  The  pressure- time  history  can  be  represented 
by  the  simple  triangular  pulse  shown  in  Figure  3,18,  thus 
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Figure  3.16.  Typical  reflected  pressure-time  history  (9). 

neglecting  the  short  rise  time  and  the  negative  phase  of 
the  blast.  The  pressure  time  histories  used  in  this 
investigation  are  illustrated  in  Figures  3.18  through 
3.20.  A  triangular  pulse  of  duration  =  50  n  sec  is 
depicted  in  Figure  3.18  and  is  used  to  investigate  the 
response  of  the  layered  systems  to  a  short  duration  blast 
wave.  A  rectangular  pulse  of  the  same  duration,  shown  in 
Figure  3.19,  is  used  to  model  the  blast  condition  in  order 
that  the  results  can  be  more  readily,  interpreted.  A 
longer  triangular  pulse  of  =  850  n  sec,  illustrated  in 
r-i -ure  3.20,  is  used  to  approximate  the  actual  blast 
durations  measured  in  reference  (7). 

The  magnitude  of  the  blast  loading  is  varied.  A  peak 
pressure  of  1,000  psi  is  used  in  the  linear  models  in 


24-63 


Figure  3.17.  Typical  reflected  pressure-time  history 
measured  in  tests  in  reference  (7). 


Figure  3.18.  Simplified  triangluar  pressure- time  history 
used  in  the  present  analysis. 

order  to  facilitate  interpretation  of  the  results.  In  the 
linear  analyses,  the  resultant  stresses  can  be  scaled  to 
any  blast  load  intensity.  A  peak  pressure  of  10,000  psi 
is  used  in  the  nonlinear  analysis  of  the  layered  systems. 
The  10,000  psi  pressure  is  representative  of  that  induced 
by  currently  e.mployed  conventional  weaponry. 
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Figure  3.19.  Rectangular  pressure- time  history  used  in 
the  present  analysis. 


TIME 


Figure  3.20.  Approximation  to  measured  pressure- time 
history  used  in  the  present  analysis. 


IV.  ANALYSIS  OF  STRESS  WAVE  PROPAGATION  BY  THE  FINITE 

ELEMENT  METHOD 

A  comprehansive  finite  element  method  (FEM) 
investigation  of  stress  wave  propagation  in  the  layered 
systems  described  in  Section  III  was  conducted.  The  FEM 
study  was  performed  through  implementation  of  the  ADINA 
(1,2, 3, 4)  computer  prograjns.  The  computer  analyses  ’’eje 
executed  on  the  Auburn  University  IBM  3083  academic 
mainframe  and  the  Alabama  Supercomputer  Network  Cray 
X-MP/24. 

General  Description  of  the  Finite  Element  Models 

Any  FEM  model  used  in  a  wave  propagation  study  must 
consist  of  a  very  fine  mesh  in  order  to  produce  accurate 
results.  For  this  reason,  only  a  very  small  cross  section 
of  each  layered  system  is  analyzed.  Each  FEM  model 
consists  of  a  single  column  of  two-dimensional,  nine-node 
isoparametric  finite  elements.  A  typical  mesh  is 
illustrated  in  Figure  4.1(a)  and  a  typical 
nine-node  isoparametric  finite  element  is  depicted  in 
Figure  4.1(b) . 

The  major  difference  between  a  wave  propagation 
problem  and  a  structural  dynamics  problem  is  the  number  of 
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LAYER  2 


LAYER  3 


(a)  Typical  FEM  Mesh 


(b)  Typical  FEM  Element 


Figure  4.1.  Typical  Finite  Element  Mesh  Consisting  of 
Nine-Node  Isoparametric  Elements. 


modes  that  significantly  contribute  to  the  response  of  the 
structure.  In  most  structural  dynamics  problems,  only  a 
few  lower  frequencies  are  excited.  Therefore,  only  these 
lower  modes  contribute  significantly  to  the  response  of 
the  system.  In  a  wave  propagation  problem,  a  large  number 
of  frequencies  are  excited.  Therefore,  a  large  number  of 
modes  significantly  contributes  to  the  structural 
response.  For  this  reason,  modal  analyses  generally  do 
not  yield  cost  effective,  accurate  results,  and  a  direct 
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numerical  integration  procedure  must  be  utilized  in  the 
wave  propagation  analysis. 

In  the  present  study,  the  Newmark  method  of  implicit 
time  integration  and  a  consistent  mass  matrix  formulation 
are  employed.  The  nonlinear  equations  are  solved  with  the 
BFGS  matrix  update  equilibrium  iteration  method  with  line 
searches.  A  modified  Newton  equilibrium  iteration  method 
is  used  in  several  instances  to  corroborate  the  results 
obtained  with  BFGS  method. 

The  time  step  selected  for  temporal  integration  in  a 
wave  propagation  problem  is  critical  to  the  accuracy  and 
stability  of  the  solution.  Since  the  Newmark  method  is 
unconditionally  stable,  selection  of  the  time  step  can  be 
baaed  entirely  upon  accuracy  considerations.  In  a  wave 
propagation  problem,  the  maximum  time  step  is  related  to 
wave  speed  in  the  material  and  element  size.  The  maximum 
time  step. is  selected  so  that  the  stress  wave  propagates 
the  distance  between  element  integration  points  within 
that  time  increment.  For  the  nine-noded  element  depicted 
in  Figure  4.1(b),  the  maximum  time  step  is  defined  by 


i^i) 


max 


(4.1) 


where  is  the  length  of  an  element  in  the  direction  of 
wave  propagation,  and  c  is  the  velocity  of  wave 
propagation.  It  has  been  determined  from  previous 
experience  (24)  that  a  time  step  of 
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(4.2) 


yields  accurate  results. 

The  investigation  of  stress  wave  propagation  in  a 
nonlinear  material  by  the  finite  element  method  requires  a 
very  fine  discretization.  Due  to  the  prohibitive 
computational  cost  required  for  the  analysis  of  a 
completely  nonlinear  layered  system,  models  consisting 
only  of  a  single,  nonlinear  material  layer  are 
investigated.  However,  the  results  of  the  nonlinear 
single  layer  models  are  used  to  predict  the  transmission 
ratios  for  the  completely  nonlinear  layered  systems,  as 
described  in  Chapter  V. 

Linear  Concrete  Model 

A  single  layer  concrete  model  is  investigated  to 
determine  the  discretization  necessary  to  produce  accurate 
results.  The  discretization  is  considered  adequate  when 
the  incident  wave  shape  is  retained  as  the  stress  wave 
propagates  through  the  FEM  model.  This  model  is  also  used 
as  a  control  to  verify  that  the  FEM  analysis  correctly 
predicts  the  .stress  wave  propagation  and  reflection 
characteristics  of  an  elastic  medium.  A  typical  mesh  for 
the  single  layer  concrete  mod^.’  xs  illustrated  in  Figure 
4.2.  The  material  properties  for  the  concrete  were 
previously  presented  in  Table  3.1. 
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24' CONCRETE  LAYER 
7S  ELEMENTS 


Figure  4.2.  Typical  Finite  Element  Mesh  for  the  Linear 
Concrete  Model. 

Nonlinear  Concrete  Model 
A  nonlinear,  single  layer  concrete  FEM  model  is 
analyzed  to  verify  that  the  wave  propagation  properties  of 
the  layered  systems  are  not  significantly  altered  by  the 
assumption  of  linear  behavior  of  the  concrete.  A  linearly 
elastic-perfectly  plastic  constitutive  law,  such  as  that 
illustrated  in  Figure  3.7,  is  used  to  approximate  the 
nonlinearity.  The  nonlinear  model  is  investigated  usi.ng 
both  the  static  and  dynamic  material  properties  of 
concrete.  The  dynamic  strength  is  taken  as  1.6  times  the 
static  strength  (15,16,18).  The  static  ar.d  dynamic 
nonlinear  concrete  parameters  are  presented  in  Table  4.1. 
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Table  4.1.  Nonlinear  Concrete  Parameters. 


".I 

'  n 


Yield 

Mass 

Poisson' s 

Strength 

Stress 

Density 

Ratio 

1  1 

E  (pst) 

(Pit) 

P  (10^) --r- 
//? 

V 

4000  psi 

¥ 

■  J 

Concrete 

3605000 

4000 

217.164 

.18 

(static) 

V 

t 

,  ( 

4000  psi 
Concrete 

4560000 

6400 

217.164 

.18 

•  < 

(dynamic) 

The  FEM  model  considered  is  a  16"  long  column 
comprised  of  200  elements  (Figure  4.3).  Several  different 
wave  types  are  examined  in  order  to  ascertain  the  effects 

'  of  the  nonlinearity  upon  wave  shape  and  wave  duration. 

! 

'  The  propagation  of  a  50  sec  short  duration  triangular 

stress  wave,  an  850  n  sec  long  duration  triangular  stress 
wave,  and  a  50  n  sec  rectangular  stress  wave  are 

■  investigated.  Stress  wave  intensities  of  7,500  and  10,000 

psi  are  considered. 

Nonlinear  Polystyrene  Model 

,  A  nonlinear  polystyrene  model  is  investigated  in 

order  to  examine  the  effects  of  the  nonlinear  absorption 
layer  upon  wave  propagation.  A  linearly  elastic-perfectly 
plastic  constit”^^/e  law,  such  as  that  illustrated  in 
Figure  3.7,  is  used  to  approximate  the  nonlinearity.  The 
nonlinear  polystyrene  parameters  were  previously  presented 
in  Table  3.3. 
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16'  C0NC3?eTE  LAYER 
200  ELEMENTS 


Figure  4.3,  Typical  Finite  Element  Mesh  for  the  Nonlinear 
Concrete  Model. 

The  FEM  model  considered  is  a  16"  long  column 
comprised  of  2C0  elements  (Figure  4.4).  Several  different 
wave  types  are  investigated  in  order  to  determine  the 
effects  of  the  nonlinearity  upon  wave  shape  and  wave 
duration.  The  propagation  of  a  50  ^  sec  short  duration 
triangular  stress  wave,  an  ,850  ^  sec  long  duration 
triangular  stress  wave,  and  a  50  sec  rectangular  stress 
wave  are  investigated.  Stress  wave  intensities  of  6  'w  and 
1,000  psi  are  considered. 
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16' POLYS-m?ENe  LAYER 
200  ELEMENTS 


Figure  4.4.  Typical  Finite  Element  Mesh  for  the  Nonlinear 
Polystyrene  Model. 

Nonlinear  Sand  Model 

A  nonlinear  FEM  sand  model  is  analyzed  in  order  to 
examine  the  effect  of  nonlinearity  in  the  exterior  sand 
layer  upon  wave  propagation.  A  linearly  elastic-perfectly 
plastic  constitutive  law,  such  as  that  illustrated  in 
Figure  3.7,  is  used  to  approximate  the  nonlinearity.  This 
constitutive  law  is  a  very  crude  representation  for  the 
’•^lilinear  response  of  sand,  but  does  provide  some  insight 
to  assessing  the  effects  of  the  nonlinearity  upon  wave 
propagation.  The  nonlinear  sand  parameters  were 
previously  presented  in  Table  3.3. 
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Th'e  FEM  model  considered  is  a  16"  long  column 
comprised  of  200  elements  (Figure  4.5).  Several  different 
wave  types  are  investigated  in  order  to  determine  the 
effects  of  nonlinearity  upon  wave  shape  and  wave  duration. 
The  propagation  of  a  50  n  sec  short  duration  triangular 
stress  wave,  an  850  n  sec  long  duration  triangular  stress 
wave,  and  a  50  n  sec  rectangular  stress  wave  are 
investigated.  Stress  wave  intensties  of  5,000  and  10,000 
psi  are  considered. 


— 


16"  SAND  LAYER 
aOO  ELEMENTS 


Figure  4.5,  Typical  Finite  Element  Mesh  for  the  Nonlinear 
Sand  Model. 
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Concrete  -  Polystyrene  -  Concrete  System 


Three  concrete-polystyrene-concrete  (CPC)  layered 
systems  are  investigated.  The  length  of  the  polystyrene 
absorption  layer  is  varied  from  2"  to  8”.  The  length  of 
the  exterior  layer  of  concrete  is  10”,  and  the  length  of 
the  interior  layer  of  concrete  is  8”.  A  typical  mesh  for 
the  CPC  system  is  illustrated  in  Figure  4.6.  The  length 
of  the  polystyrene  layer  and  the  number  of  elements  in 
each  modal  are  given  in  Table  4.2.  Typical  ADINA-IN  and 
ADINA-PLOT  data  files  for  the  CPC-4  FEM  models  are 
presented  in  Appendix  B. 

The  CPC  systems  are  investigated  using  both  linear 
and  nonlinear  constitutive  laws  for  the  polystyrene 
absorption  layer.  The  properties  for  the  linear  and 
nonlinear  materials  were  previously  presented  in  Tables 
3.1  and  3.3,  respectively.  An  incident  stress  wave  having 
an  intensity  of  1,000  psi  is  applied  to  the  linear 
systems.  A  10,000  psi  incident  stress  wave  is  applied  to 
the  CPC  layered  systems  with  the  nonlinear  polystyrene 
absorption  layer. 

Additionally,  the  nonlinear,  single  layer  analyses  of 
the  concrete  and  polystyrene  are  utilized  to  esti;rate  the 
transmission  ratios  of  the  completely  nonlinear  CPC 
layered  systems.  An  incident  stress  wave  of  10,000  psi  is 
used  in  the  nonlinear  single  layer  analyses. 
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30  ELEMENTS 


Figure  4.6.  Typical  Finite  Ele.nent  Mesh  for  the  CPC 
Layered  Systems. 


Table  4.2.  CPC  Model  Dimensions  and  Number  of  Elements. 


Total 

Length 

Exterior 

Polystyrene 

Interior 

and 

No. 

Concrete 

Absorption 

Concrete 

Elements 

•  Layer 

Layer 

Layer 

no. 

no. 

no. 

no. 

Model 

L 

el. 

el. 

^2 

el. 

^3 

el. 

CPC- 2 

20” 

90 

10” 

30 

2” 

30 

8” 

30 

CPC-4 

22” 

90 

10” 

30 

4” 

30 

8” 

30 

CPC- 8 

26” 

110 

10” 

30 

8” 

50 

8” 

30 

Sand  -  Polystyrene  -  Concrete  System 

Three  sand-polystyrene-concrete  (SPC)  layered  systems 
are  investigated.  The  length  of  the  polystyrene 
absorption  layer  is  varied  from  2”  to  8”.  The  length  of 
the  sand  layer  is  10”,  and  the  length  of  the  concrete 
layer  is  12”.  A  typical  mesh  for  the  SPC  system  is 
illustrated  in  Figure  4.7.  The  length  of  the  polystyrene 
layer  and  the  number  of  elements  in  each  model  are  given 
in  Table  4.3 . 

The  finite  element  method  is  utilized  only  in  the 
investigation  of  the  linear  SPC  systems.  The  properties 
for  the  linear  material  models  were  previously  presented 
in  Table  3.1.  An  incident  stress  wave  of  1,000  psi 
intensity  is  used  in  the  linear  investigation. 

The  transmission  ratios  for  the  nonlinear  SPC  systems 
are  estimated  using  the  results  of  the  nonlinear  single 
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Figure  4.7.  Typical  Finite  Element  Mesh  for  the  S?C 
Layered  Systems. 
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Table  4.3.  SPC  Model  Dimensions  and  Number  of  Elements. 
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SPe-2 

24" 

90 

10" 

30 

2" 

30 

12" 

30 

SPC-4 

26" 

90 

10" 

30 

4" 

30 

■>2" 

30 

SPC-8 

30" 

110 

10" 

30 

! 

8" 

50 

12" 

30 

layer  analyses.  An  incident  stress  wave  of  10,000  psi 
intensity  is  used  in  the  nonlinear  single  layer 
investigation. 

Sand-Polvstvrene-Concrete-Polvstyrene- 
Concrete  System 

Three  sand-polystyrene-concrete-polystyrene-concrete 
(SPCPC)  layered  systems  are  investigated.  However,  an  FEM 
analysis  is  not  conducted  for  this  model  due  to  the 
prohibitive  computational  costs.  The  length  of  the 
polystyrene  absorption  layer  is  varied  from  2”  to  8".  The 
length  of  the  sand  layer  is  10",  and  the  length  of  each 
concrete  layer  is  12".  A  typical  SPCPC  system  is 
illustrated  in  Figure  4.8.  The  length  of  the  l=ivers  for 
each  of  the  SPCPC  models  is  given  in  Table  4.4.  The 
properties  for  the  linear  and  nonlinear  materials  were 
previously  presented  in  Tables  3.1  and  3.3. 
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Figure  4.3.  Typical  SPCPC  Layered  System. 

The  transmission  ratios  of  the  linear  systems  is 
determined  by  the  successive  application  of  egviation 
(2.24).  The  transmission  ratios  of  the  completely 
nonlinear  SPCPC  layered  systems  is  estimated  from  the 
nonlinear  single  layer  analyses.  A  10,000  psi  intensity 
incident  stress  wave  is  assumed  for  the  nonlinear 
investigation. 
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Table  4.4.  SFCFC  Layered  System  Dimensions. 


Model' 

Total 

Length 

Ext. 

Sand 

Layer 

Foly. 

Abs. 

Layer 

Ext. 

Cone. 

Layer 

Foly. 

Abs. 

Layer 

Int. 

Cone. 

Layer 

L 

L, 

SPCFC-2 

34" 

10" 

2" 

12" 

2" 

.....  1 

8" 

1  --  . 

SFCFC-4 

38" 

10" 

IQI 

12" 

4" 

8"  i 

! 

SFCFC-8 

46" 

10" 

8" 

12" 

8" 

8"  1 

V.  RESULTS  OF  ANALYSIS 


Linear  Concrete  Model 

A  typical  stress- time  history  for  a  point  in  the 
linear  concrete  model,  described  in  Chapter  IV,  is 
illustrated  in  Figure  5.1.  The  figure  depicts  an  incident 
rectangular  compressive  stress  wave  of  intensity  1,000  psi 
and  the  corresponding  tension  wave  reflected  from  the  free 
end.  The  stress  wave  is  reflected  at  full  intensity.  The 
high  frequency  oscillations  defining  the  shape  of  the 
stress  wave  are  caused  by  the  excitation  of  high  frequency 
modes  in  the  structure.  The  theoretical  wave  propagation 
and  reflection  are  also  shown  in  Figure  5.1  as  a  dashed 
line. 

A  typical  stress  trace  from  wave  propagation  tests 
conducted  on  the  Split  Hopkinson  Pressure  Bar  (SHPB)  (22) 
is  depicted  in  Figure  5.2.  The  bar  is  loaded  axially  with 
a  high  intensity,  short  duration  pressure  wave  at  one  end. 
Each  end  of  the  bar  is  a  free  surface.  An  incident 
compressive  wave  and  two  reflected  waves  are  depicted  in 
the  figure.  The  first  reflection  is  caused  by  the 
compressive  wave  reflecting  from  the  unloaded  end  of  the 
bar  as  a  tension  wave.  The  second  reflection  is  caused 
when  the  tensile  wave  propagates  back  to  the  loaded  end  of 
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stress 


Figure  5.1.  Typical  Stress-Time  History  for  a  Point  in 
the  Linear  Concrete  Model. 


Figure  5.2.  Typical  Stress-Time  History  for  a  Point  in 
the  SHPB. 
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the  bar  and  is  reflected  as  a  compression  wave.  The 
stress  wave  is  reflected  from  each  free  surface  at  full 
intensity.  It  is  noted  that  the  stress  trace  obtained 
from  the  SHPB  contains  the  same  high  frequency 
oscillations  exhibited  by  the  FEM  analysis. 

Nonlinear  Concrete  Model 
The  mitigation  of  stress  waves  in  the  nonlinear 
concrete  model,  described  in  Chapter  IV,  is  illustrated  in 
Figure  5.3.  The  maximum  stress  resulting  from  three  types 
of  incident  stress  waves  (  1.  a  short  duration  triangular 
wave,  2.  a  short  duration  rectangular  wave,  and  3.  a  long 
duration  triangular  wave)  is  shown.  It  is  noted  that  the 
rate  at  which  the  stress  wave  degenerates  is  a  function  of 
the  shape  and  duration  of  the  incident  stress  wave.  The 
short  duration  triangular  wave  decays  very  rapidly. 

However,  both  the  short  duration  rectangular  wave  and  the 
long  duration  triangular  wave  decay  very  gradually. 

It  is  observed  that  the  inelastic  behavior  of  the 
concrete  (as  modeled  by  the  bilinear  constitutive  law 
illustrated  in  Figure  3.7)  does  not  significantly  affect 
the  propagation  of  either  the  rectangular  or  the  long 
duration  triangular  stress  waves.  The  intensity  of  the 
short  duration  triangular  stress  wave,  however,  is 
drastically  reduced  in  a  short  distance. 

The  attenuation  of  stress  waves  in  the  nonlinear 
concrete  model,  condsidering  the  dynamic  stre.ngth 
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Figure  5.3.  Maximum  Stress  Wave  Intensity  Versus  Distance 
in  Nonlinear  Concrete. 

characteristics  specified  in  Table  4.1,  is  represented  in 
Figure  5.4.  The  nonlinearity  of  the  concrete  has  a 
significant  effect  on  the  propagation  of  the  short 
duration  triangular  wave.  However,  the  rate  of  decay  of 
the  stress  wave  is  lessened.  The  intensities  of  the 
rectangular  and  long  duration  triangular  waves  are  not 
significantly  affected. 

Nonlinear  Polystyrene  Model 
The  propagation  of  a  plastic  stress  wave  in  a 
nonlinear  polystyrene  medium  is  illustrated  in  Figures  5.5 
through  5.7.  The  results  are  very  similiar  to  the 
theoretical  predictions  of  plastic  wave  propagation 
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dynamic  C0NCR6TE 


Figure  5,4.  Maximum  Stress  Wave  Intensity  Versus  Distance 
in  Nonlinear  Concrete  When  Dynamic  Strength  ■ 
is  Considered. 

presented  in  Figure  2.18.  In  Figure  5.5,  the  propagation 
of  a  short  duration  triangular  stress  wave  is  represented. 
The  plastic  portion  of  the  wave  propagates  at  a  lower 
velocity  than  the  elastic  portion  and  lags  behind  it.  As 
the  stress  wave  propagates  through  the  medium,  the  wave 
length  increases  and  the  intensity  decreases.  The  entire 
stress  wave  becomes  elastic  after  the  wave  has  propagated 
a  sufficient  distance. 

The  propagation  of  an  incident  shor-^  duration 
rectangular  stress  wave  is  depicted  in  Figure  5.6.  The 
rectangular  stress  wave  disperses  and  attenuates  with 
distance.  The  incident  rectangular  stress  wave  becomes 
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Figure  5.5.  Propagation  of  a  Plastic  Short  Duration 
Triangular  Stress  Wave  in  Polystyrene. 
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Figure  5.6.  Propagation  of  a  Plastic  Short  Duration 
Rectangular  Stress  W;ive  in  Polystyrene. 
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Figure  5.7. 


Propagation  of  a  Plastic  Long  Duration 
Triangular  Stress  Wave  in  Polystyrene. 


entirely  elastic  in  a  greater  distance  than  required  for 
the  short  duration  triangular  wave. 

The  propagation  of  a  long  duration  triangular  stress 
wave  is  illustrated  in  Figure  5.7.  The  plastic  portion  of 
the  wave  lags  behind  the  elastic  portion  but  exhibits 
little  dispersion.  The  long  duration  triangular  wave  must 
travel  a  great  distance  through  the  medium  before  it  is 
fully  mitigated  to  an  elastic  wave. 

The  attenuation  of  several  plastic  stress  waves  in 
nonlinear  polystyrene  is  illustrated  in  Figure  5.8.  It  is 
evident  from  Figure  5.8  that  the  degree  of  attenuation  of 
a  plastic  stress  wave,  as  it  propagates  through  the 
polystyrene,  is  a  function  of  the  incident  wave  shape  and 
initia'.  intensity. 

The  short  duration  triangular  wave  exhibits  a  high 
attenuation  rate  as  it  propagates  through  the  polystyrene.. 
After  a  short  distance,  the  intensity  of  the  short 
duration  triangular  stress  wave  is  reduced  to  the  yield 
stress.  The  rate  of  attenuation  exhibited  by  the  Long 
duration  triangular  wave,  however,  is  significantly  less. 

A  much  greater  distance  is  required  for  the  intensity  of 
the  long  duration  triangular  stress  wave  to  be  reduced  to 
the  yield  -tr^ss  of  th'-i  polystyrene.  The  attenuation  rate 
of  the  short  duration  rectangular  wave  is  less  than  that 
of  the  short  duration  triangular  wave,  but  is  much  greater 
than  that  of  the  long  duration  triangular  wave. 
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Figure  5.11.  Effect  of  Strength  and  Yield  Stress  on 
Maximum  Stress  Wave  Intensity  in  Sand. 


stress  axis  for  the  polystyrene  absorption  layer  and  the 
interior  concrete  layer  ranges  from  -45  to  45  psi. 

The  propagation  of  a  short  duration  rectangular 
stress  wave  is  illustrated  in  Figures  5.12  through  5.14. 
The  stress  wave  profile  at  a  time  before  the  wave 
encounters  the  boundary  between  the  exterior  concrete 
layer  and  the  polystyrene  absorption  layer  is  depicted  in 
Figure  5.12.  The  stress  profile  at  a  time  after  the  wave 
encounters  the  boundary  between  the  exterior  concrete 
layer  and  the  polystyrene  absortion  layer  is  presented  in 
Figure  5.13.  It  is  observed  that  a  very  small  percentage 
(0.96%)  of  the  stress  wave  is  transmitted  to  the 

polystyrene  absorption  layer,  and  that  the  major  portion 
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Figure  5.12.  Stress  Profile  of  the  CPC-4  Layered  System 
Sub j  ect  to  a  1000  psi  Short  Duration 
Rectangular  Stress  Wave.  Time  40.E-6  sec 


5,13*  Stress  Profile  of  the  CPC-4  Layerpd  Syste.m 
Subject  to  a  1000  psi  Short  Duration 
Rectangular  Stress  Wave.  Time  =  140. E-6 


sec 


Figure  5.14.  Stress  Profile  of  the  CPC-4  Layered  System 
Subject  to  a  1000  p'si  Short  Duration 
Rectangular  Stress  Wave.  Time  »  380. E-6  sac. 

of  the  wave  is  reflected  back  through  the  exterior 
concrete  layer  as  a  tensile  wave. 

The  stress  condition  in  the  model  at  a  time  after  the 
wave  encounters  the  interior  layer  of  concrete  is 
illustrated  in  Figure  5.14.  The  stress  wave  is 
transmitted  from  the  polystyrene  layer  to  the  interior 
layer  of  concrete  at  an  increased  intensity.  The  wave 
length  in  the  polystyrene  is  much  shorter  than  in  the 
concrete  layers  because  the  wave  speed,  .  in  the 
polystyrene  in  much  less  than  the  wave  speed,  c. ,  in  the 
concrete . 

A  second  stress  wave  propagating  in  the  polystyrene 
Layer  is  also  shown  in  Figure  5,,  14.  This  wave  is  caused 
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by  reflections  generated  in  the  exterior  layer  of 
concrete.  These  reflections  will  not  occur  in  the  actual 
structure  because  the  concrete  will  not  withstand  the 
reflected  tensile  stress.  The  exterior  layer  of  concrete 
will  spall,  thereby  diminishing  the  stress  wave  until  the 
intensity  no  longer  exceeds  the  spall  threshold  of  the 
concrete.  Additional  stress  profiles  for  the  CPC  layered 
systems  are  presented  in  Appendix  D. 

The  transmission  ratios  for  the  CPC  layered  system 
obtained  using  linear  wave  transmission  and  reflection 
theory  are  presented  in  Table  5.1.  The  theoretical 
transmission  ratios  were  obtained  using  the  BASIC  computer 
program  WAVEPROP.BAS  (Appendix  A).  The  program  results 
for  the  CPC  layered  system  are  given  in  Appendix  C.  The 
transmission  ratios  obtained  from  the  linear  FEM  analysis 
are  also  listed  in  Table  5.1.  The  results  of  the  linear 
FEM  analysis  exhibit  excellent  correlation  with  the 
predictions  of  linear  stress  wave  propagation  theory. 

Transmission  ratios  for  the  CPC  layered  systems, 
based  upon  the  single  layer  nonlinear  FEM  analyses,  are 
given  in  Table  5.2  for  a  short  duration  rectangular 
incident  stress  wave.  Plastic  wave  attenuation  in  each 
lay<^r  is  determined  directly  from  the  single  layer 
nonlinear  analyses  decribed  in  Chapter  IV.  The  stress 
transmitted  across  material  boundaries  is  estimated  from 
linear  wave  propagation  theory. 
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Table  5.1. 


Linear  Transmission  Ratios  for  the  CPC  Layered 
Systems . 


Transmission  Ratio 

Model 

Linear 

Theory 

Linear 

FEM 

Analysis 

CPC-2 

0.019 

0.018 

CPC-4 

0.019 

0.020 

CPC- 8 

0.019 

0.019 

Table  5.2.  Estimated  Nonlinear  Transmission  Ratios  for 
the  CPC  Layered  Systems  Subject  to  a  Short 
Duration  Rectangular  Incident  Stress  Wave. 


CPC-2 

CPC-4 

'  "■  . .  1 

CPC- 8  1 

Incident  Upon 
Ext.  Concrete 

10,000 

10,000 

j 

10,000  1 

Intensity 
at  Boundary 

9,900 

9,900 

9 , 900  ; 

Transmitted  to 
Polystyrene 

95 

95 

i 

95  i 

Intensity 
at  Boundary 

55 

55 

1 

55  ; 

Transmitted  to 
Int.  Concrete 

110 

110 

110  I 

j Nonlinear 

1  Transmi ssion 
j  Ratio 

.0110 

i  i 

.0110  j  .0110 

1 

t 
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Figure  5.15. 


Stress  Profile  of  the  SPC-4  Layered  System 
Subject  to  a  1000  psi  Short  Duration 
Rectangular  Stress  Wave.  Time  =  350. E-6  sec 
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Figure  5.16.  Stress  Profile  of  the  SPC-4  Layered  System 
Subject  to  a  1000  psi  Short  Duration 
Rectangular  Stress  Wave.  Time  =  550. E-6  sec 


Figure  5.17.  Stress  Profile  of  the  SPC-4  Layered  System 
Subject  to  a  1000  psi  Short  Duration 
Rectangular  Stress  Wave.  Time  *  750. E-6  sec. 

The  propagation  of  a  short  duration  rectangular 
stress  wave  is  illustrated  in  Figures  5.15  through  5.17. 
The  stress  wave  profile  at  a  time  before  it  encounters  the 
boundary  between  the  sand  and  polystyrene  layers  is 
depicted  in  Figure  5.15.  The  stress  profile  at  a  time 
after  the  boundary  between  the  exterior  sand  layer  and  the 
polystyrene  absortion  layer  is  encountered  is  represented 
in  Figure  5.16.  It  is  observed  that  a  small  percentage 
(7.14%)  of  the  stress  wave  is  transmitted  to  the 
polystyrene  -br-^rption  layer,  and  that  the  major  portion 
of  the  wave  is  reflected  back  through  the  exterior  sand 
layer  as  a  tensile  wave. 
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The  stress  condition  in  the  model  at  a  time  after  the 
wave  encounters  the  interior  layer  of  concrete  is 
illustrated  in  Figure  5.17.  The  stress  wave  is 
transmitted  from  the  polystyrene  layer  to  the  interior 
layer  of  concrete  at  an  increased  intensity.  Additional 
stress  profiles  for  the  SPC  layered  systems'  are  presented 
in  Appendix  E. 

The  transmission  ratios  for  the  SPC  layered  systems 
obtained  using  linear  wave  transmission  and  reflection 
theory  are  listed  in  Table  5.5.  The  theoretical 
transmission  ratios  were  obtained  using  the  BASIC  computer 
program  WAVEPROP.BAS  (Appendix  A).  The  program  results 
for  the  SPC  layered  system  are  given  in  Appendix  C.  The 
transmission  ratios  obtained  from  the  linear  FEM  analysis 
are  also  listed  in  Table  5.5.  The  results  of  the  linear 
FEM  analysis  exhibit  excellent  correlation  with  the 
predictions  of  linear  wave  propagation  theory. 

Transmission  ratios  for  the  SPC  layered  systems, 
based  upon  the  single  layer  nonlinear  FEM  analyses,  are 
given  in  Table  5.6- for  a  short  duration  triangular  stress 
wave.  Plastic  wave  attenuation  in  each  layer  is 
determined  directly  from  the  single  layer  nonlinear 
analyses  described  in  Chapter  IV.  The  stress  trar. -..v  -ted 
across  material  boundaries  is  estimated  from  linear  wave 
propagation  theory. 
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Table  5.5.  Linear  Transmission  Ratios  for  the  SPC  Layered 
Systems . 


Transmission  Ratio 

Model 

Linear 

Theory 

Linear 

FEM 

Analysis 

SPC-2 

0.142 

0.130 

SPC-4 

0.142 

0.150  . 

SPC-8 

0.142 

0.150 

Table  5.6.  Estimated  Nonlinear  Transmission  Ratios  for 
the  SPC  Layered  Systems  Subject  to  a  Short 
Duration  Triangular  Incident  Stress  Wave. 


SPC-2 

SPC-4 

SPC-8 

Incident  Upon 
Exterior  Sand 

10,000 

10,000 

j 

10,000  1 

Intensity 
at  Boundary 

8,900 

8,900 

1 

8,900  ! 

Transmitted  to 
Polystyrene 

635 

635 

j 

635 

Intensity 
at  Boundary 

274 

.  152 

1 

55  ! 

Transmitted  to 
Int.  Concrete 

545 

303 

1 

1 

110  1 
1 

Nonlinear 

Transmission 

Ratio 

1 

.0545 

.0303 

1 

.0110  j 

24-104 


The  transmission  ratios  estimated  from  the  single 
layer  nonlinear  FEM  analyses  are  summarized  in  Table  5.7. 
The  transmission  ratio  decreases  as  the  thickness  of  the 
polystyrene  layer  is  increased.  A  plot  of  transmission 
ratio  versus  the  thickness  of  the  polystyrene  absorption 
layer,  for  both  the  linear  and  nonlinear  materials,  is 
shown  in  Figure  5.18.  It  is  indicated  that  a  thicker 
polystyrene  absorption  layer  increases  the  capacity  of  the 
layered  system  to  mitigate  an  incident  stress  wave, 
thereby  increasing  the  structures  resistance  to  interior 
spalling. 


Table  5.7.  Estimated  Nonlinear  Transmission  Ratios  for 
the  SPC  Layered  Systems  Subject  to  a  10,000 
psi  Incident  Stress  Wave. 


Estimated  Nonlinear  Transmission  Ratio 

SPC-2 

SPC-4 

SPC-8 

Short  Duration 
Triangular 

Stress  Wave 

.0545 

.0303 

.0110 

Short  Duration 
Rectangular 
Stress  Wave 

.1260 

.0967 

.0496 

Long  Duration 
Triangular 

Stress  Wave 

.1304 

.  1226 

.1107 
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SPC  LAYERED  SYSTEM 


Figure  S.18.  Transmission  Ratio  Versus  the  Thickness  of 
the  Polystyrene  Absorption  Layer  of  the  SPC 
Layered  System. 

Sand«Poly3tvrene-Concrete~Polvstvrene- 
Concrete  System 

The  transmission  ratios  for  the  SPCPC  layered 
systems,  based  upon  the  single  layer  nonlinear  FEM 
analyses,  are  given  in  Table  5.8  for  a  long  duration 
triangular  incident  stress  wave.  Plastic  wave  attenuation 
in  each  layer  is  determined  directly  from  the  single  layer 
nonlinear  analyses  described  in  Chapter  IV.  The  stress 
transmitted  across  material  boundaries  is  estimated  from 
linear  wave  propagation  theory. 

The  transmission  ratios  estimated  from  the  single 
layer  nonlinear  FEM  analyses  are  summarized  in  Table  5.9. 
It  IS  noted  that  the  transmission  ratio  decreases  as  the 
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Table  5.8.  Estimated  Nonlinear  Transmission  Ratios  for 

the  SPCPC  Layered  Systems  Subject  to  a  10,000 
psi  Long  Duration  Triangular  Incident  Stress 
Wave. 


SPCPC-2 

SPCPC-4 

SPCPC-8 

Incident  Upon 
Exterior  Sand 

10,000 

10,000 

10,000 

Intensity 
at  Boundary 

9,900 

9,900 

9,900 

Transmitted  to 
Polystyrene 

707 

707 

707 

Intensity 

at  Boundary 

-  -  . . . 1 

655 

616 

556 

Transmitted  to 
Ext.  Concrete 

1304 

1226 

1107 

- . -  I 

Transmitted  to 
Polystyrene 

13 

12 

Transmitted  to 
Int.  Concrete 

26 

1 

24  1 

1 

22 

Nonlinear 

Transmission 

Ratio 

.0026 

.0024 

.0022 

thickness  of  the  polystyrene  layer  is  increased.  A  plot 
of  transmission  ratio  versus  the  thickness  of  the 
polystyrene  absorption  layer  is  shown  in  Figure  5.19.  It 
is  indicated  that  a  thicker  polystyi-ene  absorption  layer 
increases  the  capacity  '"■*'  the  layered  system  to  mitigate 
an  incident  stress  wave,  thereby  increasing  the  structures 
resistance  to  interior  spalling. 
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Table  5.9. 


Estimated  Nonlinear  Transmission  Ratios 
the  SPCPC  Layered  Systems  Subject  to  a  10,000 
psi  Incident  Stress  Wave. 


Estimated  Nonlinear  Transmission  Ratio 


Short  Duration 
Triangular 
Stress  Wave 

Short  Duration 
Rectangular 
Stress  Wave 


SPCPC-2 


.0010 


.0024 


SPCPC- 4 


.0006 


SPCPC-8 


.0002 


.0018 


.0010 


Long  Duration 
Triangular 
Stress  Wave 


.0026 


.0024 


.0022 


The  thickness  of  the  polystyrene  absorption  layer  has 
no  effect  on  the  transmission  ratio  of  the  linear  system 
(see  Figure  5.19).  The  linear  transmission  ratio  for  the 
SPCPC  layered  systems,  obtained  from  the  BASIC  computer 
program  WAVEPROP.BAS  (Appendix  A),  is  0.0027.  The  program 
results  for  the  SPCPC  layered  system  are  given  in 
Appendix  C. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


Conclusions 

The  finite  element  method  (FEM)  has  been  successfully 
employed  to  predict  the  propagation,  reflection,  and 
transmission  of  stress  waves  in  both  homogeneous  and 
layered  media.  However,  very  fine  discretizations  and 
small  integration  time  steps  are  required  in  order  for  the 
FEM  analyses  to  yield  accurate  results.  This  has  resulted 
in  significant  computational  costs,  even  for  simple  FEM 
models. 

It  has  been  shown  that  layered  media  significantly 
alter  stress  wave  propagation  characteristics.  The  proper 
combination  of  high  impedance  and  low  impedance  materials 
has  demonstrated  substantial  enhancement  of  stress  wave 
attenuation.  This  finding  suggests  that  the 
implementation  of  layered  media  in  the  construction  of 
protective  shelters  can  eliminate  blast  induced  spalling 
at  the  interior  surface  of  the  shelter  walls. 

The  concrete-polystyrene-concrete  (CPC),  sand- 
polystyrene-concrete  (SPC),  and  sand-polystyrene-concrete- 
polystyrene-concrete  (SPCPC)  layered  systems  have  been 
shown  to  be  extremely  effective  in  mitigating  an  incident 
stress  wave  as  it  propagates  through  the  syscem.  The 
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transmission  ratios  for  these  layered  systems,  calculated 
on  the  assumption  of  linear  material  behavior,  indicate 
that  their  deployment  in  the  construction  of  protective 
shelters  will  eliminate  blast  induced  spalling.  The 
transmission  ratios  for  these  systems,  calculated  on  the 
assumption  of  inelastic  behavior  of  the  media,  indicate  an 
even  higher  degree  of  attenuation  for  propagating  stress 
waves,  thus  lending  additional  credibility  to  the  layered 
system  concept. 

The  SPCPC  layered  system  exhibits  the  highest  degree 
of  stress  wave  attenuation.  The  linear  transmission  ratio 
for  this  system  is  calculated  to  be  0.0027.  The  nonlinear 
transmission  ratio  of  the  SPCPC-8  layered  system  is 
estimated  to  be  as  low  as  0.0002  for  a  short  duration 
triangular  stress  wave,  and  0.0022  for  a  long  duration 
triangular  stress  wave. 

The  CPC  layered  system  also  exhibits  a  significant 
degree  of  stress  wave  attenuation.  The  linear 
transmission  ratio  for  the  CPC  systems  is  0.019,  The 
nonlinear  transmission  ratio  of  the  CPC  layered  system  is 
estimated  to  be  as  low  as  0.0098  for  a  short  duration 
triangular  stress  wave,  and  0.0110  for  a  long  duration 
triangular  stress  wave. 

The  SPC  layered  system  is  the  least  effective  in 
terms  of  stress  wave  mitigation.  The  linear  transmission 
ratio  for  the  SPC  layered  system  is  0.142.  The  nonlinear 
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transmission  ratio  of  the  SPC-8  system  is  estimated  to  be 
0.011  for  a  short  duration  triangular  stress  wave  and 
0.1107  for  a  long  duration  stress  wave.  The  SPG  systems 
considered  ware  comprised  of  a  10"  thick  layer  of  sand.  A 
structure  surrounded  by  a  sand  be  m,  or  a  buried 
structure,  would  be  tantamount  to  a  sand  layer  much 
thicker  than  10".  Therefore,  the  actual  nonlinear 
transmission  ratios  for  the  SPG  systems  would  be  much  less 
than  those  determined  for  the  SPG  models  investigated  in 
this  study. 

In  all  the  layered  systems;  examined,  the  thickness  of 
the  absorption  layer  significantly  effects  the  nonlinear 
transmission  ratio.  A  thicker  polystyrene  absorption 
layer  results  in  a  reduced  nonlinear  transmission  ratio. 

Recommendations 

The  FEM  analyses  conducted  in  this  study  indicate 
that  the  GPG,  SPG,  and  SPGPG  layered  systems  are  highly 
effective  in  reducing  transmitted  stresses.  However,  the 
transmission  ratios  were  obtained  without  consideration  of 
the  dynamic  properties  of  the  concrete,  sand,  and 
polystyrene.  Also,  the  nonlinear  soil  and  concrete 
approximations  of  linear  elastic-perfectly  plastic 
constitutive  relations  are  very  crude. 

The  present  study  has  examined  only  the  wave 

propagation  and  attenuation  characteristics  of  the  CPC, 

SPC,  and  SPCPC  layered  systems.  It  is  recommended  that 
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additional  numerical  studies  be  conducted  to  investiqate 
the  overall  structural  response  of  protective  shelters 
employing  the  layered  system  concept. 

It  is  recommended  that  field  tests  of  the  CPC,  SPC, 
and  SPCPC  layered  systems  be  conducted  to  determine  their 
effectiveness  in  reducing  spall.  It  is  recommended  that 
layered  systems  consisting  of  polystyrene  layers  of  4", 

8”,  and  12”  thicknesses  be  tested. 

The  SPC  layered  system  is  the  most  cost-effective 
alternative  for  retrofitting  existing  protective  shelters. 
It  is  recommended  that  extensive  field  testing  be 
conducted  on  the  SPC  layered  systems  to  determine  the 
optimum  thickness  of  the  polystyrene  absorption  layer. 
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